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Abstract 


Photorefractive  materials  are  an  important  class  of  electrooptic  materials  used  in 
a  wide  variety  of  nonlinear  optical  devices.  Although  recent  theories  explaining  the 
photorefractive  charge  transport  properties  of  these  materials  have  been  relatively 
successful,  the  origin  of  the  charge  centers  is  unknown.  Also,  there  are  several 
anomalous  effects  observed  in  BaTi03  which  indicate  that  the  charge  transport 
models  are  incomplete. 

In  this  investigation,  photo-induced  absorption,  and  isotropic  and  anisotropic 
photorefractive  gratings  were  used  to  study  the  photo-induced  charge  transport 
properties  of  BaTi03  in  an  attempt  to  identify  the  photorefractive  species.  Previous 
models  and  experiments  on  BaTi03  have  not  taken  into  account  the  effects  of  photo- 
induced  absorption  in  studying  photorefractive  gratings.  Also,  most  of  this  early 
work  concentrated  on  the  isotropic  nature  of  these  gratings  to  determine  the  material 
properties.  The  experiments  and  analyses  presented  here  extend  these  techniques, 
both  theoretically  and  experimentally,  to  include  photo-induced  absorption  and 
anisotropic  photorefractive  gratings. 

The  photo-induced  absorption  experiments  presented  here  represent  evidence  for 
the  existence  of  a  distribution  of  localized  trapping  states  in  the  band  gap  of  BaTi03. 
The  intensity  dependence  of  the  photo-induced  absorption  correlates  well  with  the 
intensity  dependences  of  the  trap  density  and  recombination  time,  quantities  which 
until  now  were  assumed  to  be  constants. 

Several  novel  photorefractive  scattering  processes  were  also  studied.  Isotropic 
scattering  of  a  single  intense  pump  beam  was  observed  and  modeled.  Numeri¬ 
cal  calculations,  performed  using  material  parameters  consistent  with  the  previous 
experiments,  yielded  scattering  cross  sections  in  excellent  agreement  with  the  ob¬ 
servations.  Anisotropic  scattering  of  a  single  extraordinary  polarized  beam  into 
conical  rings  of  light  was  also  observed.  Using  a  phase  matching  condition  and  the 
electrooptic  tensor  symmetry  properties  of  BaTi03,  the  scattering  process  of  the 
inner  ring  was  uniquely  determined.  A  model  which  describes  the  formation  of  the 
outer  ring  was  also  proposed. 

Anisotropic  scattering  into  higher  orders  were  observed.  To  derive  approximate 
solutions  for  the  higher  order  gratings  a  perturbation  technique  was  used  to  modify 
the  transport  equations.  The  measurements  were  found  to  be  well  described  by 
the  model.  Finally,  a  new  anisotropic  four-wave  mixing  configuration  was  demon¬ 
strated.  It  was  shown  that  amplified  reflection  could  be  performed  without  phase 
conjugation. 
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1  Introduction 

1.1  The  Photorefractive  Effect 


The  wide  spread  proliferation  of  laser  technology  and  the  growing  interest  in  non¬ 
linear  optics  during  the  early  sixties,  led  to  an  intensive  search  for  new  optical  ma¬ 
terials.  This  early  research  was  primarily  concerned  with  parametric  wave-mixing, 
such  as  frequency  doubling.1'2  and  parametric  amplification.3’4  Since  these  processes 
required  a  variety  of  phase  matching  conditions,  many  birefringent  crystals,  such 
as  ammonium  dihydrogen  phosphate  (ADP),5  potassium  dihydrogen  phosphate 
(KDP).  6  and  barium  sodium  ninbafe  Ba2NaNb50i5,7  were  examined.  Unfortu¬ 
nately,  the  high  optical  intensities  required  for  optical  parametric  processes  were 
also  a  source  of  optical  damage.6 

In  LiNb03  and  LiTa03,  this  damage  was  manifest  as  distortions  in  the  refractive 
index  along  the  path  of  the  beam.8  which  could  be  reversed  by  heating  the  crystal, 
or  by  flooding  it  uniformly  with  light.  This  damage  was  also  found  to  be  a  function 
of  total  energy  incident  on  the  crystal.  That  is,  short  exposure  times  and  high 
intensities  gave  the  same  effect  as  long  exposure  times  and  low  intensities.  The 
ability  of  the  crystal  to  integrate  the  optical  power  indicated  that  some  type  of 
microscopic  transport  process  could  be  involved.  Since  the  refractive  index  damage 
was  photo-induced,  the  term  photorefractive  was  coined. 

In  a  similar  phenomenon,  known  as  the  photochromic  effect,  the  optical  absorp¬ 
tion  could  be  changed  by  the  application  of  light.  This  change  in  the  absorption 
was  shown  to  be  the  result  of  photo-induced  charge  transfer  between  impurity  ions 
in  the  crystal.9,10  Since  different  valence  states  of  the  ions  had  different  absorption 
spectra,  spatially  varying  intensity  patterns  generated  spatially  varying  absorption 
centers  in  the  crystal.  An  extension  of  this  model  was  used  by  Chen11  to  explain 
the  optical  damage  seen  in  LiNb03. 

In  Chen’s  model,  Fig.  1,  charge  trapped  in  defects  or  impurities  in  the  band 
gap,  is  photoexcited  in  the  high  intensity  regions  of  the  crystal.  The  liberated  car¬ 
riers  then  migrate  to  the  dark  regions  of  the  crystal,  resulting  in  a  redistribution 
of  charge  with  the  same  spatial  variation  as  the  incident  light.  This  redistributed 
charge  has  an  associated  space  charge  field  which,  in  noncentrosymmetric  crystals, 
induces  a  change  in  the  refractive  index  through  the  linear  electrooptic  effect.  Since 
the  symmetry  requirement  for  phase  matching  is  the  same  as  that  for  the  electroop¬ 
tic  effect,  no  inversion  symmetry,  many  materials  of  possible  use  in  parametric- 
wave  mixing  were  found  to  be  subject  to  photorefractive  damage.  Some  of  the  ma¬ 
terials  found  to  be  photorefractive  include  ferroelectric  BaTi03,12  SrxBa1_J.Nb206 
(SBN),13  Ba2NaNb50i5,14  and  KNb03,15  non-ferroelectric  Bii2Si02o  (BSO),16  and 
GaAs,1,  and  organic  polymers.18  Also  certain  PLZT  ceramics,19,20  which  possess 
only  a  second  order  (Kerr)  electrooptic  effect,  have  been  found  to  be  photorefrac¬ 
tive. 
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Figure  1:  Basic  mechanism  of  the  photorefractive  effect.  Photoexcited  electrons  are 
redistributed,  resulting  in  a  space-charge  field  which  modulates  the  refractive  index 
through  the  electrooptic  effect. 

Chen  was  the  first  to  devise  the  standard  experiment  for  studying  the  photore¬ 
fractive  effect.  He  realized  that  the  refractive  index  variations  could  be  used  to 
store  phase  holograms  in  the  crystal.21  In  his  experiment  ,  the  interference  fringes  of 
two  coherent  plane  waves  intersecting  in  the  crystal,  produced  a  pure  phase  grating 
through  the  photorefractive  effect.  This  enables  characterization  of  the  effect  in  a 
crystal  by  measuring  the  diffraction  efficiency  of  a  light  induced  grating.  In  one 
experiment,  more  than  100  gratings  were  recorded  in  a  LiNb03  crystal.22 

Since  this  early  work,  the  photorefractive  effect  has  been  used  in  a  large  number 
of  nonlinear  optical  applications.  The  most  widely  known  of  these  is  optical  phase 
conjugation.23,24,25  This  will  be  discussed  in  Chapter  6.  Other  applications  include; 
image  subtraction,26  associative  holographic  memories,2'  phase-locking  of  lasers,28 
and  optical  bistability.29  While  many  of  these  applications  have  shown  great  promise 
in  the  optical  signal  processing  field,  it  is  a  promise  which  has  so  far  been  fruitless, 
since  these  applications  are  confined  to  relatively  few  research  laboratories.  Some  of 
the  reasons  photorefractive  materials  have  not  been  successfully  integrated  into  use¬ 
ful  devices  are  (1 )  a  lack  of  a  clear  understanding  of  the  photorefractive  centers  and 
their  role  in  the  photorefractive  process,  and  (2)  several  anomalous  photorefractive 
scattering  processes  that  tend  to  degrade  device  performance.  Other  considerations 
are  more  in  the  nature  of  engineering  problems  such  as  the  requirements  of  stability 
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•  ABSORPTION  COEFFICIENT 
■  EMPTY  TRAPS  (BEAM  COUPLING) 
A  F*3  *  CONCENTRATION  (EPR) 


10 


4 

# 


R2 


i 

G  L 


♦  * 
R  1 

1334 


40  80  120  1  BO 

TOTAL  IRON  CONCENTRATION,  ppm 


E  E 


*“2 


I  £ 


z  - 

3  « 


Figure  2:  Correlation  between  iron  concentration  and  the  photorefractive  trap  den¬ 
sity  in  BaTi03.  (After  Ref.  34) 

and  coherence  needed  for  holography. 

Identification  of  the  photorefractive  species  in  LiNb03  was  attempted  in  the 
early  seventies  by  Peterson  et  al.  30  Their  method  was  to  combine  a  number  of 
techniques  such  as  electron  spin  resonance,  oxidation  and  reduction,  and  optical 
absorption  to  correlate  the  affects  of  selected  doping  on  the  photorefractive  prop¬ 
erties.  They  measured  the  diffraction  efficiency  of  undoped  crystals  and  crystals 
doped  with  450  ppm  Fe.  The  diffraction  efficiency  of  the  doped  crystals  was  found 
to  be  two  orders  of  magnitude  larger  than  the  undoped  crystals.  After  annealing 
in  a  reduced  oxygen  partial  pressure  the  diffraction  efficiency  of  all  the  crystals  was 
smaller  than  that  of  the  as-grown  or  air- annealed  crystals. 

Their  ESR  measurements  indicated  that  the  dominant  valence  state  of  the  iron 
in  the  as  grown  crystal  was  Fe+3.  Since  the  reduced  crystals  showed  both  a  decrease 
in  the  diffraction  efficiency  and  a  decrease  in  the  Fe+3  ESR  signal,  they  concluded 
that  the  predominant  photorefractive  species  in  LiNb03  was  Fe+3.  Unfortunately, 
most  of  the  early  photorefractive  work  did  not  include  light-induced  grating  erasure 
and  although  there  has  been  some  work  on  the  erasure  rate  kinetics,  31  ■  32,  33  the 
intensity  dependence  of  the  erasure  rate  has  not  been  reported.  Therefore,  it  is  not 
known  if  the  decay  rate  of  LiNb03  exhibits  a  sublinear  intensity  response.  Also, 
photo-induced  absorption  effects  were  not  reported. 

The  first  attempts  to  determine  the  photorefractive  species  in  BaTi03  were 
only  recently  attempted  by  Klein  and  Schwartz.  34  Their  research  was  patterned 
after  that  of  Peterson  et  al,  30  with  the  exception  that  they  used  beam  coupling 
to  characterize  the  photorefractive  effect.  Chemical  analyses  of  the  crystals  showed 
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that  iron  was  the  predominant  impurity  in  their  commercial  samples,  although  other 
transition  metal  impurities  were  also  present .  These  iron  concent  rations  also  seemed 
to  correlate  to  the  Fe+3  ESR  signals  from  the  samples.  Since  the  beam  coupling 
results  showed  a  linear  correlation  to  the  iron  concentration,  they  concluded  that 
iron  was  also  the  dominant  photorefractive  species  in  BaTi03.  Their  data  showed 
a  four-fold  increase  in  the  trap  density  for  only  a  two  fold  increase  in  the  iron 
concentration  Fig.  2. 

This  result  however,  is  in  direct  conflict  with  the  results  of  Schunemann  et  al.  3' 
and  the  work  presented  in  this  thesis.  The  origin  of  this  conflict  is  not  completely 
understood,  however  the  crystals  used  in  this  work  were  grown  from  purified  mate¬ 
rials  and  systematically  doped  with  varying  amounts  of  iron,  constituting  a  better 
control  group. 

Ducharme  and  Feinberg  3,1  conducted  oxidation  reduction  experiments  on  similar 
commercially  available  samples  of  BaTi03.  They  found  that  the  dominant  carrier, 
which  was  holes  in  as  grown  and  oxidized  crystals,  could  be  converted  to  electrons 
after  reduction  at  650‘c  at  a  partial  pressure  of  10-<>  atm.  The  sublinear  intensity 
dependence  of  the  decay  rate  usually  observed  for  BaTiO.j,  3l  was  observed  for  their 
as  grown  and  oxidized  crystals.  However,  the  crystals  which  had  been  converted  to 
n  type  exhibited  a  linear  intensity  dependence  for  the  decay  rate. 
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2  Barium  Titanate:Physical  Properties  and  Crystal 
Growth 

The  electrical  and  optica!  properties  of  Barium  titanate  have  been  studied  inten¬ 
sively  over  the  past  forty  years  with  much  of  this  work  performed  at  the  MIT  Lab¬ 
oratory  for  Insulation  Research.  7'his  section  primarily  serves  to  give  background 
on  the  material  properties  of  BaTi03  relevant  to  the  photorefractive  effect. 

2.1  Physical  Properties 

Above  130  °C,  barium  titanate  exists  in  the  cubic  perovskite  structure,  AB03, 
shown  in  Fig.  2.1.  1  In  this  structure,  the  Ti4+  ions  (B),  are  at  the  center  of  six  02~ 
ions  forming  regular  octahedrons.  The  octahedrons  arc  connected  at  the  corners 
with  Ba2+  ions.  (A),  occupying  holes  in  this  framework.  Each  Ba2+  ion  is  surrounded 
by  twelve  equidistant  02~  ions. 

As  the  temperature  is  lowered  below  130  °C,  BaTi03  experiences  three  first  order 
ferroelectric  phase  transitions.  At  130  °C  the  unit  cell  becomes  distorted  along  a 
[100]  direction  and  the  crystal  symmetry  reduces  to  tetragonal  4mm.  It  is  believed 
that  the  Ti  and  0  ions  move  to  new  equilibrium  positions  in  opposite  directions  [sec 
Fig.  4]  resulting  in  a  net  polarization.  This  model  was  shown  to  fit  neutron  and 
x-ray  diffraction  data  with  the  barium  ions  not  contributing  to  the  polarization.  1 
Since  this  shift  can  occur  in  any  one  of  six  directions,  cooling  through  the  Curie 
temperature  causes  simultaneous  nucleation  of  multiple  domain  structures  in  the 
crystal. 

The  first  order  clectrooptic  effect  has  it’s  origin  in  the  quadratic  electrooptic 
effect.  To  understand  this  consider  the  cubic  structure  of  BaTi03  above  130  C'C.  In 
this  phase  the  refractive  index  is  isotropic.  However,  if  we  apply  an  electric  field  to 
the  crystal  we  can  induce  birefringence  through  the  quadratic  electrooptic  effect  2 

4(1),!-  (!) 


In  this  case  the  polarization  P  is  the  result  of  the  applied  field  E. 

Now  consider  the  same  BaTi03  crystal  below  130  °C.  In  this  tetragional  phase 
there  is  a  spontaneous  polarization  P3  and  from  Eq.  1  we  expect  that  the  crystal 
should  be  naturally  birefringent.  When  an  electric  field  E  is  applied  to  the  crystal 
the  net  polarization  P  is  composed  of  both  the  spontaneous  polarization  Ps  and  the 
polarization  due  to  applied  fields  P/,-.  From  Eq.  1  we  see  that  the  change  in  the 
refractive  index  now  has  several  terms  with 


A 


oc  P *  +  Pf:  +  Pa  PE  . 


(2) 


Figure  3:  Perfect  perovskite  structure  showing  (a)  the  unit  cell  centered  around  the 
Ti'M  ion,  and  (b)  the  TiOe  octahedra  surrounding  the  Ba+2  ion.  (After  Ref.  1). 
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Figure  4:  Ionic  displacement  in  tetragonal  BaTi03. 


As  stated,  the  natural  birefringence  is  due  to  the  spontaneous  polarization  P,,  with 
n0-ne  ocPj.  For  small  applied  fields  we  can  neglect  the  P\  term.  This  leaves  the 
cross  term  which  is  known  as  the  first  order  electrooptic  effect.  Taking  into  account 
the  tensoral  properties  of  the  crystal,  the  electrooptic  effect  can  be  written  as  3 

An*  =  -y  rijk  Ek  (3) 

where  r is  known  as  the  electrooptic  tensor.  Using  the  symmetry  properties  of 
the  tetragonal  phase  4mm,  the  allowed  tensor  elements  are  ri31=r3U=r232=r322, 
rn3=r223>  and  r333.  3 

The  two  remaining  phase  transitions  are  found  to  be  very  destructive.  Below 
5  °C,  the  unit  cell  distorts  to  mm  orthorombic  symmetry  with  the  spontaneous 
polarization  along  the  [110]  directions.  Finally,  below  -90  °C,  the  spontaneous 
polarization  is  along  the  [111]  direction  with  3m  trigonal  symmetry. 

The  room  temperature  tetragonal  phase  is  the  one  of  interest  to  most  optics 
researchers  since  its  both  ferroelectric  and  electrooptic  with  only  a  only  a  slight 
elongation  of  the  unit  cell  in  the  polar  direction. 

2.2  Crystal  Growth  and  Preparation 

At  present,  the  best  technique  of  growing  single  crystalline  BaTi03  is  the  top-seeded 
solution  growth  technique  developed  at  the  MIT  Crystal  Physics  and  Optical  Elec¬ 
tronics  Laboratory.  In  this  method,  a  solution  of  BaO  and  Ti02  is  prepared  using 
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the  phase  diagram  shown  in  Fig.  6(a).  A  composition  of  66%Ti02  and  34%BaO  is 
found  to  give  the  widest  growth  range.  The  furnace,  Fig.  2.2,  is  heated  by  silicon 
carbide  heating  elements  and  the  thermal  gradients  are  kept  low  by  insulating  fire 
bricks.  The  temperature  is  raised  to  approximately  1400  °C  and  a  [100]  seed  of 
BaTiOa  is  introduced  into  the  furnace  above  the  melt.  After  equilibrium  is  reached, 
the  seed  is  introduced  into  the  melt.  To  initiate  growth,  the  temperature  is  slowly 
lowered  at  approximately  0.4  °C/hr.  After  5  to  10  hours,  the  seed  is  pulled  from 
the  melt  at  a  rate  of  0.2  to  1.0  mm/hr.  After  about  60  hours  of  pulling,  the  crystal 
is  removed  from  the  melt  and  the  furnace  slowly  cooled  to  room  temperature. 

Because  of  the  large  differences  in  the  size  of  Ba  and  Ti  ions,  doping  is  expected 
to  be  highly  site  selective.  For  instance,  the  transition  metals  are  all  similar  in 
radius  to  the  Ti  ion,  while  the  rare  earth  ions  are  closer  in  radius  to  the  Ba  ion. 
Therefore  Fe+3  would  be  expected  to  be  incorporated  in  the  lattice  at  the  Ti+4  site, 
while  Ce  would  probably  reside  in  the  Ba  site. 

The  first  attempt  at  doping  was  made  on  crystals  grown  from  a  melt  containing 
reagent  grade  Ti02  and  BaC03.  These  crystals  were  light  green  in  color  before  dop¬ 
ing.  The  addition  of  20  ppm  Cr  to  the  melt,  {0.66[(l-x)  Ti02+x  C'r203]  +0.34  BaO], 
resulted  in  crystals  which  were  very  dark  red  in  color.  Although  these  were  not  use¬ 
ful  for  the  normal  photorefractive  studies,  the  strong  beam  fanning  discussed  in 
Chapter  6  was  first  seen  in  these  crystals. 

The  second  doping  attempt  was  made  on  crystals  known  as  nominally  pure,  that 
is,  they  were  grown  from  a  melt  containing  relatively  pure  Ti02  and  reagent  grade 
BaC03.  Crystals  of  this  type  are  the  same  as  those  reported  in  the  literature.  These 
crystals  were  found  to  have  a  light  yellow-orange  color.  The  addition  of  20  ppm  Ce 
to  the  melt,  [0.66TiO2+0.34(l-x)BaO+0.34(x)CeO2],  resulted  in  moderately  dark 
red  samples.  Strong  beam  fanning  was  observed  in  the  Ce  doped  samples  and 
photorefractive  characterization  was  again  hindered. 

The  third  attempt  at  doping  was  performed  by  Peter  Schunemann  on  a  melt 
containing  Ti02  which  he  purified  through  a  distillation  process  and  BaC03  purified 
by  David  Gabbe  through  a  solvent  extraction  processes.  The  undoped  crystals 
showed  very  little  absorption  in  the  visible.  The  doped  crystals  were  prepared  by 
substituting  Fe  for  Ti  in  the  melt  [0.66(l-x)Ti02+0.66(x/2)Fe203+0.34BaO]  where 
x  is  the  doping  level.  The  melt  doping  levels  were  50  ppm,  500  ppm,  and  1000  ppm. 
The  color  of  the  as-grown  doped  samples  varied  from  light  yellow  for  the  50  ppm  to 
brown  for  the  1000  ppm.  The  visible  absorption  in  50  ppm  and  500  ppm  crystals 
was  small  enough  to  allow  photorefractive  characterization. 

Crystal  growth  from  a  BaTi03  Ba204  melt  [see  Fig.  6(b)]  was  also  attempted. 
This  has  the  potential  of  increasing  the  growth  range  and  lowering  the  growth  tem¬ 
perature.  In  this  case  only  one  growth  run  was  attempted.  The  composition  used 
was  [0.4B2O3+0.6TiO2  +  BaO]  and  the  resulting  seeding  temperature  was  1287°. 
This  melt  seemed  to  be  extremely  viscous  and  the  resulting  crystal  was  simply  a 
mass  of  dendritic  growth.  However,  recently  crystals  have  been  successfully  grown 


Figure  5:  Schematic  diagram  of  a  typical  furnace  used  for  top- seeded  solution  growth 
of  oxide  crystals.  (After  Ref.  4). 
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(a)  (b) 

Figure  6:  (a)  The  Ba0-Ti02  phase  diagram  and  (b)  the  BaB204-BaTi03  phase 
diagram.  (After  Refs.  6  and  7) 

from  a  similar  melt  by  David  Gabbe.  5 

As  a  final  note,  a  weighing  error  by  the  author  lead  to  a  melt  with  approximately 
a  50/50  Ba0/Ti02  ratio.  This  resulted  in  a  crystal  which  was  extremely  cracked 
and  very  unsymmetric.  X-ray  diffraction  data  taken  by  David  Gabbe  showed  this 
crystal  to  have  a  hexagonal  structure  as  expected  by  the  phase  diagram  Fig.  6(a). 

Since  the  Curie  point  is  around  130  °C,  crystals  cut  from  the  boules  are  found 
to  have  twinning  of  90°  domains  as  shown  in  Fig.  7.  These  crystals  are  cut  along 
[100]  directions  and  polished  first  with  20  micron  grit  and  then  3  micron  diamond 
paste.  The  domains  are  eliminated  by  two  methods.  First  the  samples  are  subjected 
to  hydrostatic  pressure  alternately  along  two  perpendicular  axes  to  eliminate  the 
twinning.  After  each  cycle  the  crystal  must  be  repolished  due  to  surface  changes 
which  can  result  in  cracking  if  the  poling  were  continued.  This  mechanical  poling 
removes  the  90°  domains  leaving  only  the  domains  which  are  180°  apart.  These 
are  aligned  by  heating  the  crystal  near  or  above  the  Curie  temperature,  applying 
an  electric  field  of  approximately  1000  V/cm,  and  then  slowly  cooling  the  sample 
back  to  room  temperature.  Figure  8(a)  shows  the  typical  domain  structure  of  an 
unpoled  crystal  placed  between  cross  polarizers.  Figure  8(b)  shows  the  same  crystal 
after  mechanical  and  electrical  poling  with  the  c-axis  parallel  to  the  page  and  the 
long  axis  of  the  crystal.  The  size  difference  is  due  to  the  successive  repolishing  of 
the  crystal  during  the  poling  cycles. 
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Figure  7:  Dtvnain  twinning  in  BaTiOa-  The  arrows  denote  the  polar  axis  (After 
Ref.  1). 


Figure  8:  Photographs  of  a  BaTi03  crystal  (a)  unpoled  and  (b)  poled.  The  crystal 
was  placed  between  crossed  polarizers  and  in  (b)  the  c-axis  is  parallel  to  the  page. 
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Three  samples  doped  with  50  ppm,  500  ppm,  and  1000  ppm  Fe  were  electroded 
on  the  c  faces,  placed  in  a  quartz  tube  and  annealed  at  800  °C  in  an  argon  at¬ 
mosphere  with  100  ppm  oxygen  (lO-4  atm)  for  36  hours.  The  samples  were  then 
quenched  to  175  °C  by  moving  the  quartz  tube  to  another  furnace  at  175  °C.  A  pol¬ 
ing  field  was  then  applied  and  the  temperature  slowly  lowered  to  room  temperature 
at  a  rate  of  approximately  5  °C/hour.  9<  8 
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3  Role  of  Iron  Centers  in  the  Photorefractive 
Effect  in  Barium  Titanate 
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Variably  valent  iron  impurities  have  been  suggested  as  the  source  of  photorefractive  charge  carriers  in  BaTiO 
High-purity  BaTiO,  crystals  were  grown  with  transition-metal  impuritv  levels  below  the  He'  cm  baseline 
estimated  for  photorefractivitv  Iron  doped  trvstals  were  grown  with  concentrate  of  5.  .50.  1150.  500.  750.  and 
1000  parts  in  10*'  of  iron  Changes  in  iron  \alence  monitored  by  optical  absorption  in  the  Pi,  range  1  10  1  atm  were 
found  to  be  consistent  with  a  defect  chemical  model  indicating  Fe"  and  Be'*  to  be  the  stable  valence  states  in  this 
range.  Photorefractive  characterization  of  high-purity  BaTiO,  suggests  that  variably  valent  iron  ions  are  not  the 
dominant  photorefractive  species  in  these  crystals,  whereas  the  role  of  iron  centers  in  doped  crystals  is  complicated 
by  large  absorption  losses. 


INTRODUCTION 

BaTiOi  is  of  particular  interest  among  photorefractive  ma¬ 
terials  because  it  has  demonstrated  the  highest  four-wave 
mixing  reflectivities  (greater  than  a  factor  of  201  by  virtue  of 
its  large  electro-optic  coefficients.1  Despite  the  promise  of 
barium  titanate  for  applications  involving  phase  conjugation 
and  optical  information  processing,  the  centers  responsible 
for  the  photorefractive  effect  in  this  material  have  not  been 
conclusively  identified,  and  this  lack  of  understanding 
hinders  the  control  and  optimization  of  its  nonlinear-optical 
and  photorefractive  properties. 

The  basic  physics  of  the  photorefractive  effect  is  currently 
well  understood1:  (1)  When  the  material  is  nonuniformly 
illuminated,  carriers  (electrons  or  holesi  are  optically  excit¬ 
ed  from  energy  levels  (donors  or  acceptors)  in  the  band  gap; 
(2)  the  excited  carriers  migrate  by  diffusion  and  drift  and  are 
retrapped  and  reexcited  until  they  reach  the  less  bright 
regions  of  the  material;  (3)  this  redistribution  of  charge  cre¬ 
ates  internal  electric  fields  that  modulate  the  refractive  in¬ 
dex  of  the  material  through  the  electro-optic  effect.  The 
exact  identity  of  the  photorefractive  sources  and  traps,  how¬ 
ever,  remains  unclear  for  BaTiOi. 

Previous  efforts  to  isolate  the  photorefractive  centers  in 
BaTiO  i  have  focused  mainly  on  nominally  undoped  crystals 
from  a  single  commercial  source.4  Klein  and  Schwartz  ’  de¬ 
termined  the  impurity  concentrations  in  a  series  of  such 
crystals  by  chemical  analysis  and  noted  a  general  increase  in 
the  photorefractive  trap  density,  Fe1*  electron  paramagnet¬ 
ic  resonance  (EPR)  signal,  and  absorption  coefficient  with 
increasing  iron  content  over  a  narrow  concentration  range 
[--50- 150  parts  in  10fi  (ppmlj.  On  the  basis  of  their  observa¬ 
tions  and  similar  results  reported  for  LiNbO,."  "  Klein  and 
Schwartz  concluded  that  Fe1*  and  Fe:*  are  the  sources  and 
'raps,  respectively,  of  photorefractive  carriers  when  the  car¬ 
riers  are  holes  (and  vice  versa  for  electrons).  Reduction  of 
the  same  crystals  (Po;  =  10“"  atm)  yielded  increases 

in  photorefractive  speed  and  a  change  in  carrier  sign  from 
positive  to  negative.  These  were  explained  bv  a  model 


based  on  variations  in  the  Fe-'* 'TV1*  ratio.1''  Ducharmeand 
Feinberg11  reported  the  effect-  of  oxidation  and  reduction  at 
650°C  on  a  nominally  undoper!  crystal  and  also  noted  a 
change  in  photorefractive  carrier  Ggn.  which,  in  their  crys¬ 
tal,  corresponded  to  a  minimum  in  trap  density  near  0.5  atm. 
These  results  led  them  to  conclude  that  oxygen  vacancies  (or 
impurity-oxygen  vacancy  complexes)  were  the  photorefrac¬ 
tive  donors  ot  negative  charge 

The  only  reported  study  of  the  photorefractive  properties 
of  intentionally  doped  BaTiOi  is  the  recent  work  of  Gode- 
froy  et  al.1-  They  measured  the  maximum  diffraction  effi¬ 
ciency,  writing  speed,  and  electro-optic  coefficients  in  a  se¬ 
ries  of  crystals  doped  with  a  wide  range  of  Fe  concentrations 
(undoped  to  2500  ppm).  They  reported  that  all  measured 
parameters  were  highest  in  the  750-ppm  Fe-doped  crystal, 
with  the  diffraction  efficiency  increasing  by  a  factor  of  1000 
over  that  of  the  undoped  crystal.  A  major  shortcoming  of 
this  work,  however,  is  that  the  measured  crystals  were  nei¬ 
ther  electricallv  nor  mechanically  poled,  and  the  measure¬ 
ments  that  were  performed  are  particularly  sensitive  to  the 
presence  of  90°  and  180°  domains.  In  fact,  their  reported 
investigations  of  the  crystal-domain  structure  (performed 
using  a  unique  optical  topography  approach)  indicated  a  low- 
degree  of  poling  and  notable  variation  among  crystals. 

The  goal  of  this  research  was  to  gain  more  insight  into  the 
role  of  Fe  centers,  if  any.  in  the  photorefractive  process  in 
BaTiOi  by  controlling  the  concentration  and  valence  slates 
of  the  Fe  in  a  series  of  crystals  and  by  correlating  the  results 
with  measured  optical  and  photorefractive  properties.  Our 
approach  is  based  on  the  popular  theory'’  that  variably  va¬ 
lent  transition-metal  ions,  in  particular  divalent  and  triva- 
lent  Fe  impurities,  are  the  species  responsible  for  photore¬ 
fractive  charge  carriers  in  BaTiO ,. 

Our  strategy  for  evaluating  the  effects  of  such  centers  on 
the  photorefractive  behavior  of  BaTiO,  was  fourfold  Be 
cause  it  has  been  estimated  that  the  density  of  charge  earn 
ers  involved  in  the  photorefractive  effect  is  of  the  order  of 
10;,;/cm 1  the  first  goal  of  this  research  was  to  prepare  high- 
purity  single  crvstals  with  transition-metal  and  aliovalent 
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impurity  levels  below  this  level.  The  second  goal  of  this 
work  was  to  grow  a  series  of  Fe-doped  oaTiO single  crystals 
with  low  background  impurity  levels  by  systematically  dop¬ 
ing  the  melt  from  which  the  above-mentioned  high-purity 
crystals  were  grown.  Third,  we  employed  a  defect-chemical 
approach  to  control  the  Fe  valence  through  thermal  oxida¬ 
tion  and  reduction  treatments.  Finally,  these  variations  in 
Fe  concentration  and  valence  were  correlated  with  general 
trends  in  the  photorefractive  behavior  of  the  crystals  as 
characterized  by  saturation  diffraction  efficiency,  light-in¬ 
duced  grating  erasure  rate,  and  two-beam  coupling  mea¬ 
surements.  Preliminary  results  of  this  study  were  recently- 
presented. 14  1,1 

HIGH-PURITY  BARIUM  TITANATE  SYNTHESIS 
Purification  of  Starting  Materials 

A  major  effort  in  this  research  was  devoted  to  the  synthesis 
of  high-purity  starting  materials  because  titanium  dioxide 
and  barium  carbonate  feed  materials  are  not  commerciallv 
available  with  impurity  levels  below  the  nominal  concenfa 
tion  of  lO^/cm1,  as  is  apparently  required  tor  photoretrac 
tivity. 

The  TiOi  feed  material  was  prepared  by  the  hydrolysis  of 
titanium  isopropoxide,  TicOCiH  tj,  according  to  the  reac¬ 
tion 

C  H  OH 

Til'OC  ,H:)4  +  :iH  .0  *■  TiO  -  H  O  +  4(  .H  OH 

ill 

The  titanium  alkoxide  precursor  (Stauffer  Chemical  Com¬ 
pany)  was  first  purified  by  a  reduced-pressure  fractional 
distillation.  This  purification  technique  was  based  on  the 
observation  that  titanium  isopropoxide  has  a  much  lower 
boiling  point  than  the  isopropoxides  of  other  transition- 
metal  impurities.17  A  600  mm  x  12.5  mm  vacuum  jacketed 
distilling  column  packed  with  1/8-in.  (0.32-cm)  glass  helices 
provided  -'15-20  theoretical  plates. 1M  Purification  was 
qualitatively  evidenced  by  the  transformation  of  the  alkox¬ 
ide  from  deep  yellow  to  water  clear  in  color,  leaving  a  Uark- 
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(a)  PURIFICATION  Ot  TITANIUM  DIOXIDE 


!b)  PURIFICATION  or  BARIUM  CARBONATE 

Fig  I  Flow  charts  for  'a1  preparation  ot  high  purtn  I  it)  anti  'hi 
purification  and  ntheM-  of  Hal  (l 


Fig  i.  Typical  single  crystal  boule  grown  by  TSSG. 


Fig  3.  Photograph  of  90°  domain  walls  (diagonal  lines!  in  unpoled 
3aTiO  i  single  crystal. 


brown  liquid  and  metal  particulates  in  the  still  pot  The 
collected  distillate  was  removed  under  a  dry-nitrogen  atmo¬ 
sphere  and  was  reacted  with  a  dilute  solution  of  deionized 
water  and  semiconductor-grade  isopropyl  alcohol  (Nlal- 
linckrodt)  to  form  hydrated  TiO.  powder,  which  was  then 
vacuum  dried  in  the  same  vessel.  A  flow  chart  for  the  pro¬ 
cess  is  shown  in  Fig  Hal. 

Barium  carbonate  was  formed  by  reaction  of  a  purified 
aqueous  solution  of  BaCl.  and  urea.  ( NHi.CO.  with  COj  in  a 
pressurized  vessel  Reagent  grade  BaCl.  and  urea  were  dis¬ 
solved  in  water,  filtered,  and  purified  hv  solvent  extraction. 
In  the  solvent -extraction  process  the  aqueous  solution  was 
intimately  mixed  with  a  dense  immiscible  liquid  phase 
(chloroform)  that  contained  an  extracting  or  chelating  agent 
(diethylammonium  diethyidithiocarbonate).  which  readily 
forms  complexes  with  a  large  number  of  metal  ions.  These 
metal  complexes  segregated  to  the  dense  chloroform  phase, 
where  their  presence  was  indicated  by  a  bright  yellow  or 
green  color,  and  were  subsequently  removed  by  using  a  se 
paratory  funnel  The  extraction  process  was  repeated  for 
several  solution  pH  values  until  the  extracting  phase  w?s 
colorless.  After  extraction,  the  purified  solution  was  refil- 
tered  and  reacted  under  CO.  pressure  (1.03  x  10'  Fa, 
~8<)°C)  to  form  crystalline  BaCO,.  which  was  centrifuged, 
washed,  and  dried.  A  flow  chart  for  the  entire  process, 
carried  out  in  a  class  lot'  clean  room,  is  shown  in  Fig  1  <bi. 
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Single-Crystal  Growth 

High-purity  BaTiCU  single  crystals  were  grown  by  the  top- 
seeded  solution  growth  (TSSG)  technique.19-0  Extreme 
care  was  exercised  to  minimize  furnace  contamination  of  the 
purified  feed  materials.  The  growth  furnace  was  rebrcked 
with  new  lining  materials,  new  SiC  heating  elements  were 
used,  and  the  platinum  crucible  was  kept  cover' d  while  the 
closing  bricks  were  lowered  in  place.  Crystal  growth  was 
first  initiated  on  a  platinum  wire  to  avoid  contamination 
from  an  impure  single-crystal  seed.  Later,  high-purity  sin¬ 
gle-crystal  seeds  cut  in  the  100  ■  direction  were  used  with 
dimensions  of  approximately  10  mm  X  2.5  mm  X  2.5  mm. 

Single  crystals  were  grown  from  a  titanium-rich  melt  con¬ 
taining  66  mol  °i  TiCF.  corresponding  to  a  liquidus  tempera¬ 
ture  of  approximately  1380°C.  (Experience  has  indicated 
that  growth  from  melts  at  temperatures  >1400°C  yields 
poor  results.)  The  melt  was  heated  slowly  at  first  to  drive 
out  moisture  and  was  then  raised  from  200  to  1450°C  in  20  h 
and  soaked  for  24-40  h  to  ensure  complete  dissolution  of  the 
ceramic  powders.  Seeding  temperatures,  indicated  by  equi¬ 
librium  between  the  melt  and  the  air-cooled  seed,  varied 
between  1382  and  1392°C  for  successful  growth  runs  because 
of  variations  in  the  moisture  content  of  th.  hydrated  TiOj. 
After  seeding,  the  melt  was  cooled  at  0.2°C/h  for  25-35  h 
until  pulling  began  at  0.15-0.2  mm/h.  During  pulling,  the 
melt  was  cooled  at  0.5°C/h  to  temperatures  slightly  above 
the  eutectic  (1320°C),  after  which  the  crystal  was  removed 
from  the  melt  and  slowly  cooled  to  room  temperature.  Most 
growth  runs  yielded  single-crystal  boules  of  good  optical 
quality,  as  shown  in  Fig  2. 
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Optical  Sample  Preparation 

Crystals  were  cut  along  |100|  planes,  as  determined  by  tht- 
backreflection  Laue  method,  and  faces  were  ground  to  near 
parallel  with  20-Mm  SiC.  Cut  and  ground  cubes  were  then 
polished  by  hand  with  3-gm  diamond  paste  on  a  lightly  oiled 
paper  lap.  The  presence  of  90°  domains  in  the  polished 
samples  was  clearly  indicated  by  sharp  planar  boundaries 
within  the  polished  samples.  These  boundaries  ran  diago- 
nally  across  the  •  redominantly  a  faces,  as  shown  in  Fig.  3. 
These  domain  walls,  in  fact,  served  to  identify  the  a  faces, 
whereas  the  c  direction  was  indicated  by  a  characteristic 
color  pattern  when  viewed  between  crossed  polarizers. 

Single-domain  samples  were  achieved  by  a  combination  of 
mechanical  and  eleci  rical  poling.  Uniaxial  stresses  of  10.3  x 
106  to  20.7  X  10H  Pa  applied  along  an  a  axis  were  generally 
sufficient  to  remove  all  visible  90°  domain  walls  During 
application  of  the  stress,  performed  by  using  a  pressurized 
piston,  the  90°  domain  walls  visible  in  the  unstressed  u 
direction  would  migrate  toward  the  corners  of  the  crystal. 
This  migration  would  create  surface  steps  on  faces  along 
which  the  ends  of  the  domain-wall  boundaries  moved.  The 
moving  domain  walls  would  eventually  get  hung  up  on  these 
surface  steps  so  that  the  crystal  needed  to  be  repolished 
before  poling  could  continue.  After  sufficient  iterations,  all 
directly  visible  90°  domain  walls  couH  be  removed,  although 
a  cloudiness  visible  in  crossed  polarizers  often  revealed  the 
presence  of  residual  90°  domains.  These  were  removed  by- 
applying  an  electric  field  of  1000  V/cm  in  an  oil  bath  near  or 
above  the  Curie  temperature  (133°C).  Nickel  electrodes 
were  polished  to  optical  flatness  to  ensure  good  electrical 


Table  1.  Spark-Source  Mass  Spectrographic  Analysis  of  Reagent  Precursors,  Purified  Feed  Materials,  and 

Resulting  BaTiOj  Single  Crystals  ' 


Element1’ 

TiO. 

Cnpurified 

Purified 

BaCI. 

(Cnpurified) 

Crea 

(Cnpurified) 

BaCO, 

(Purified) 

BaTiO, 

Li* 

0.5 

2 

1 

B  * 

2 

3 

2 

l>  2 

1 

5 

Na* 

<2 

<  2 

2 

o.r> 

2 

<2 

Mg  * 

2 

2 

i 

<  i 

5 

2 

Al" 

; 

•) 

2 

2 

3 

0.5 

Si1* 

5 

7) 

4 

i" 

4 

4 

P  * 

0  1 

0.03 

0.2 

n  i 

0.1 

0.05 

S-  - 

<2 

<2 

<2 

Do 

<2 

<2 

Cl- 

10 

3 

Major 

2 

150 

10 

K* 

2 

2 

<0  2 

o  i 

2 

<0  1 

Ca" 

3 

5 

200 

250 

2 

Cr1-* 

006 

<0.02 

<do  2 

0.05 

<0  02 

Mn- 

<0.02 

<-0.02 

Fe  " 

O 

().-> 

0.3 

0  2 

1 

0.3 

NT  '* 

0,3 

0  1 

<0.04 

<0.04 

0  04 

Cu:  * 

0  1 

<0  02 

<n  02 

0.06 

Zn  * 

<0  1 

0.1 

As'* 

0.5 

0.05 

0  1 

Sr  * 

30 

20 

r>  o 

Zr" 

0  1 

<0  04 

0  1 

NV 

003 

La1* 

0.U6 

0  06 

Ce" 

<0  02 

<0  02 

Pr  " 

<0  02 

Nd1* 

<0  08 

Pb-‘ 

<0.04 

<0  04 

7  Quantities  are  given  in  part*  per  million  atomic 
h  Si1*,  ("«-*  Sr-  *.  Zr1*,  and  Ph  ‘  are  isovaient  impurities 
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Fig  4  Measured  TV  spectra  of  reagent  BaCI  and  purified  B.i<  O 
dissolved  in  H<’1  Beak  at  J  in  nm  corresponds  to  Fe  impurities 
removed  by  solvent  extraction 


Fig  5  Measured  optical  absorption  spectra  comparing  relative 
purity  of  BaTiO  single  crystals  grown  from  commercial  and  puri¬ 
fied  feed  materials. 


contact,  and  voltage  was  applied  slowly  ( 10  V/min)  to  reduce 
the  tendency  toward  surface  cracking  at  the  positive  elec¬ 
trode.  which  sometimes  occurs  during  rapid  domain  reorien¬ 
tation. 

Critical  factors  in  achieving  poled  BaTiO  samples  proved 
to  be  (1)  having  a  good  polish  so  that  moving  domain  walls 
would  not  hang  up  on  surface  flaws,  (2)  using  multiple  itera¬ 
tions  during  mechanical  poling  to  avoid  stressing  and  crack¬ 
ing,  and  13)  having  stress-free  single  crystals.  Residual 
stresses  from  cooling  through  the  phase  transition  during 
crystal  growth  or  stresses  induced  during  mechanical  polish¬ 
ing  and  poling  could  sometimes  prevent  complete  mechani¬ 
cal  poling  and  would  often  cause  severe  cracking  during 
electrical  poling 

Evaluation  of  Sample  Purity 

A  combination  of  spark  source  mass  spectrometry  (North¬ 
ern  Analytical  Laboratory,  Amherst,  New  Hampshire).  UV 
spectrophotometry,  and  optical-absorption  spectroscopy 
were  used  to  confirm  the  impurity  concentration  levels  in 
the  starting  materials,  the  purified  feed,  and  the  resulting 
high-purity  BaTiO  single  crystals.  The  results  of  spark- 
source  mass  spectrometry  are  listed  in  Table  1,  and  the  TV 
spectrophotometric  data  and  visible  absorption  spectra  are 
given  in  Figs  4  and  5.  respectively  The  results  of  these 
analyses  can  be  summarized  as  billows: 


1.  The  spark -source  data  indicated  that  the  TiCb  frac¬ 
tional  distillation  procedure  reduced  the  Fe  content  hv  a 
factor  of  4  and  removed  all  other  transition  metals  to  below 

detection  limits. 

2.  L  \  spectrophotometric  data  of  purified  BaCO  ,.nd 
the  reagent  BaCL  precursor  dissolved  in  HC1  (Fig.  4)  indi¬ 
cated  that  the  absorbance  due  to  transition-metal  chloro- 
complexes  in  the  BaCO:l  is  significantly  lower  than  in  the 
BaCf,.  The  peak  at  230  nm  in  the  BaCL  spectra  corre¬ 
sponds  well  to  the  Fe  peak  reported  in  the  literature.-1  In 
contrast,  then,  to  the  mass  spectrometric  data,  which  gave 
Fe  concentrations  of  0.3  and  1  ppm  in  Bad.,  and  BaCO,. 
respectively,  these  results  show,  that  the  BaCO  purification 
process  significantly  reduced  the  Fe  levels  of  the  Bad  pre¬ 
cursor. 

3.  Optical-absorption  spectra  of  crystals  grown  from  (at 
reagent  BaCO,  and  Baker  Ultrex  TiOj,  (b)  reagent  grade 
BaCO,  and  purified  TiOj,  and  (c)  purified  BaCO:)  and  puri- 
fiec  TiO.  synthesized  in  this  laboratory  (Fig.  5)  reveal  that 
the  band  edge  of  the  high-purity  crystal  is  shifted  to  a  slight¬ 
ly  shorter  wavelength  and  drops  off  more  steeply,  both  of 
which  are  indications  of  decreased  impurity  levels. 

4.  As  indicated  by  mass  spectrometry,  transition-metal 
impurity  levels  in  the  high-purity  BaTiO.3  single  crystal  are 
near  or  below  the  nominal  concentration  of  10"7cm3  (1  ppm 
atomic  =  1.56  X  lO'Vcm3)  predicted  for  observing  the  pho- 
torefractive  effect.  Fe  and  Ni  were  measured  at  concentra¬ 
tions  of  0.3  and  0.04  ppm  (atomic),  respectively,  whereas  all 
other  transition  metals  were  below  detection  limits  (0.01 
ppm).  The  only  aliovalent  impurities  reported  at  levels 
greater  than  this  were  B  (5  ppm),  Al  (0.5  ppm),  and  Cl  (10 
ppm),  which  is  a  remnant  of  the  BaCL  precursor.  These 
results  reveal  that  these  high-purity  BaTiOn  crystals  are,  to 
our  knowledge,  the  purest  yet  achieved  in  this  laboratory  or 
ever  used  for  photorefractive  studies. 

SYSTEMATIC  IRON  DOPING 

To  investigate  systematically  the  effects  of  Fe  concentration 
on  the  photorefractive  properties  of  BaTiCL,  crystals  were 
grown  from  melt  compositions  containing  5.  50,  250,  500, 
750,  and  1000  ppm  of  Fe.  Doping  was  accomplished  by 
substituting  FeO  for  TiOj  in  the  above-mentioned  high- 
purity  melt  according  to  the  proportion  34%  BaCO:l  +  66^ 
[( 1  -  x)Ti02  +  (t/2)Fe.>Oi],  where  x  is  the  doping  level  in 
parts  in  10h  (atomic).  The  Fe  concentration,  therefore,  re¬ 
fers  to  the  ratio  Fe/Ti  in  the  melt.  The  actual  Fe  concentra¬ 
tions  in  the  resulting  crystals  were  confirmed  by  chemical 
analysis.  Atomic  absorption  analysis  (Northern  Analytical 
Laboratories)  was  conducted  on  the  doped  BaTiOi  samples, 
which  were  dissolved  by  fusing  with  lithium  tetraborate. 
The  results  are  listed  in  Table  2  and  indicate  a  segregation 
coefficient  near  unity.  This  value  for  the  segregation  coeffi¬ 
cient  is  somewhat  higher  than  that  reported  by  Godefroy  et 
al.2-  for  TSSG  BaTiO  1.  Their  Fe  concentrations,  however, 
were  substantially  higher  than  those  used  in  this  study,  and 
the  few  data  points  at  lower  concentrations  do  suggest  values 
approaching  unity. 

Single-domain  optical  samples  of  the  Fe-doped  crystals 
were  prepared  using  the  procedures  described  above  for  the 
case  of  high-purity  crystals.  Two  notable  differences  in  be¬ 
havior,  however,  were  observed  during  the  electrical-poling 
procedure.  Fir-t.  crystals  containing  500.  750.  and  1000 
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Table  2.  Chemical  Analysis  of  Fe-Doped  BaTiO^ 
Single  Crystals" 


Sample  Number 

Iron  Doping  Level 
in  Melt  (ppm) 

Iron  Level  in 

Crystal.  Analyzed  (ppmC 

1434 

0 

0.3' 

1442 

5 

<10 

1446 

50 

49 

1453 

500 

530 

1457 

750 

731 

1462 

10OU 

980 

J  Results  indicate  a  segregation  coefficient  near  unity. 

*  Atomic  absorption  i  Northern  Analytical* 

Spark -source  mass  spectrometry  \  Northern  Analytical! 


ppm  of  Fe  were  found  to  develop  an  unusual  bright  yellow 
color  band  during  application  of  the  applied  field  (1000  \7 
cm)  at  temperatures  near  the  Curie  temperature.  T  This 
color  band  was  defined  by  a  sharp  boundary  parallel  to  the 
positive  electrode,  and  with  time  this  boundary  propagated 
toward  the  negative  electrode.  A  similar  coloration  phe¬ 
nomenon  was  apparently  observed  by  Godefroy  et  a/.1-’  and 
prevented  them  from  electrically  poling  their  samples.  For¬ 
tunately,  in  these  experiments  the  color  front  advanced  only 
1  or  2  mm  into  the  crystal  in  the  times  required  for  electrical 
poling  (1.5-2. 5  h),  thereby  leaving  a  region  large  enough  for 
optical  evaluation. 

Electrocoloration  phenomena  were  previously  reported  in 
both  BaTiO!  (Refs.  23  and  24)  and  SrTiO:t.-'s  Blanc  and 
Staebler-"  noted  that  the  application  of  dc  electric  fields  to 
transition-metal-doped  SrTiO.i  in  the  temperature  range  of 
~100-325°C  resulted  in  the  appearance  of  colored  regions 
characteristic  of  oxidized  material  at  the  positive  electrode 
and  reduced  material  at  the  negative  electrode.  Their  re¬ 
sults  were  consistent  with  a  simple  model  based  on  the  drift 
of  doubly  ionized  oxygen  vacancies  under  the  influence  of 
the  applied  field,  which  caused  oxidation  and  reduction  of 
stationary  transition-metal  ions  near  the  positive  and  nega¬ 
tive  electrodes,  respectively.  This  effect  was  observed  for 
all  transition  metals  studied:  Fe,  Ni,  Co.  and  Mo.  Assum¬ 
ing  that  the  same  phenomenon  occurs  in  our  heavily  doped 
BaTiOcFe  samples,  we  applied  their  analysis  to  our  results 
to  obtain  an  approximate  value  for  the  mobility  of  oxygen 
vacancies,  uv,  in  BaTiO  t  given  by 

dar/df  =  mv  V’/x.  (2) 

where  dx/dt  is  the  rate  of  the  color  boundary  movement,  x  is 
the  width  of  the  untransformed  region,  and  V  is  the  applied 
voltage.  This  yielded  a  n\  value  of  ~3  X  10'"  cm-A’  sec  at 
~128°C,  which  compares  favorably  with  the  1.5  X  10"'’  cm -7 
V  sec  obtained  by  Blanc  and  Staebler-"  for  SrTiO  <  at  200°C. 

The  second  effect  of  added  Fe.  which  was  noted  during 
electrical  poling,  was  a  slight  decrease  in  the  Curie  tempera¬ 
ture  evidenced  by  visual  observation  of  the  tetragonal-to- 
cubic  phase  transition.  More-accurate  Curie  temperature 
measurements  were  performed  using  capacitance  tech¬ 
niques  on  thin  (100)  slabs  cut  from  high  purity,  50-,  500 -. 
and  1000-opm  Fe-doped  samples.  These  samples  were 
equilibrated  at  800°C  in  oxygen  partial  pressures  of  1.  10"-’, 
and  10"'  atm.  and  in  each  case  they  were  subsequently 
quenched  to  ~1T5°C  and  slowly  cooled  to^l25°C  while  the 
ac  capacitance  (100  kHz)  was  monitored.  In  each  case  the 
phase  change  was  marked  by  a  sudden  plunge  in  the  relative 


dielectric  constant  from  ~  12.000  in  the  cubic  phase  to 
~2,000  after  the  transition.  The  results  indicated  that  the 
Curie  temperature  dropped  approximately  2°C  per  loon 
ppm  of  Fe.  in  close  agreement  with  the  value  of  —  21°C’ 
mol  ^ )  reported  by  Hagemann  and  lhrig  Annealing  in  the 
oxygen  partial  pressure  range  indicated  above  had  no  signifi¬ 
cant  effect  on  T, ,  however,  which  is  in  agreement  with  the 
results  of  previous  workers-'*’  as  well  as  those  of  Weschler  et 
a/.1"  but  is  in  sharp  contrast  with  the  results  of  Ducharme 
and  Feinberg,11  who  noted  a  6°C  drop  in  T  for  annealing  in 
argon  as  compared  with  oxygen. 

CONTROL  OF  IRON  VALENCE 
Defect  Chemistry  of  BaTiOyFe 

The  understanding  and  control  of  the  various  valence  states 
of  Fe  and  their  relative  concentrations  is  clearly  critical  to 
evaluating  the  possible  role  of  these  centers  in  the  photore- 
fractive  process  in  BaTiO;.  The  defect  chemistry  of 
Ba-TiO.i,-’  ’-  and  in  particular  transition-metal-doped 
BaTiO "  IS  has  been  extensively  studied  by  high-tem- 
perature  electrical  conductivity,  thermogravimetry.  Moss- 
bauer  spectroscopy,  and  EPR.  The  work  of  Hagemann  and 
co- workers-8  11  ,r'  on  Fe-doped  BaTiO.;  yielded  the  thermo¬ 
dynamic  defect  model  presented  below,  which  permits  the 
concentrations  of  Fe->+,  Fe:!+,  and  Fe’"  to  be  determined  as  a 
function  of  oxygen  partial  pressure,  temperature,  and  total 
Fe  doping  level. 

When  BaTiOt  is  treated  at  low  oxygen  partial  pressures 
and  elevated  temperatures,  oxygen  vacancies  are  created 
and  balanced  by  the  generation  of  mobile  electrons  accord¬ 
ing  to  the  reaction 

<V  =  V.O,  +  IV  *  +  2e\  (3) 

with  an  equilibrium  constant  for  reduction  given  by 

Kn  =  [VV  ’ ]n-P„  1  ’  =  K„-  exp<-Atf„7 kT),  (4) 

where  n  is  the  electron  concentration  and  A Hn“  is  the  stan¬ 
dard  enthalpy  of  reduction  associated  with  the  formation  of 
doubly  ionized  oxygen  vacancies.  The  notation  is  that  of 
Kroger  and  Vink.'8  where  the  cross,  dot,  and  slash  super¬ 
scripts  correspond  to  neutral,  positive,  and  negative  effec¬ 
tive  charges,  respectively.  The  electroneutrality  condition 
in  this  regime  is  simply 

2[V(;*]  =  n.  (5) 

where  the  atmosphere-controlled  oxygen -vacancy  concen¬ 
tration  is  much  greater  than  the  total  Fe  concentration. 

At  higher  oxygen  partial  pressures,  however,  the  oxygen- 
vacancy  concentration  becomes  fixed  by  the  total  Fe  concen¬ 
tration.  depending  on  the  stable  valence  state(s)  present. 
For  every  two  Fe1'  ions  reduced  to  Fe  <+,  one  oxygen  vacancy 
is  formed,  whereas  for  every  Fe-4  ion  formed,  one  oxygen 
vacancy  is  generated.  Consequently  the  condition  of  charge 
balance  requires  that 

(FeTl'j  =  2(1’,,*  '|  (6) 

when  trivalent  and  tetravalent  iron  are  present  or 

I  Fe ,  ’]  +2|Fe-,r|  =  2|V0”I  •') 

at  lower  f',,  where  both  Fe"  and  Fe  '  are  present  and  are 
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represented  by  Fe-r/  and  Fe-n ",  respectively.  The  relevant 
ionization  reactions  are  given  by 

FeT,"  =  FeTl'  +  e'  (8) 

and 

FeT,'  =  Fe,/  +  e'.  (9) 

These  reactions  have  the  corresponding  mass-action  rela¬ 
tions14  ts 

[FeTl']n/(FeXl")  =  K„”  =  1/2.Y,  exp[-t£,  -  Ey/~)ik T). 

(10) 

(FeTl*l«/|FeTl']  =  K„'  =  2.Vrexp[-(£,  -  E¥e'))kT\.  (Ill 

where  .NY  is  the  density  of  states  in  the  conduction  band 
( 1.55  X  10-’-  cm'3),  £c  is  the  band  gap  of  BaTiO,  (3.1  eV),  and 
£ f.-~  and  E¥t~  are  the  energy  levels  of  Fe!*  and  Fe4*.  re¬ 
spectively,  as  measured  from  the  top  of  the  valence  band. 

The  concentration  of  Fe  ions  in  a  given  valence  state  is 
related  to  the  total  Fe  concentration  through  the  mass-bal¬ 
ance  equation 

[FeTi|,l>t  -  [FeTl']  +  [FeTl*]  (12) 

and 

[FeTl]„„  =  [FeTl1  +  [FeTt']  (13) 

for  oxidizing  and  reducing  conditions,  respectively. 

The  equilibrium  defect  model  represented  by  Eqs.  (3)— 
(13)  can  be  used  to  predict  the  relative  concentration  levels 
of  Fe2*,  Fe3*,  and  Fe4”'  quantitatively  for  a  given  set  of 
annealing  conditions  and  doping  levels  if  values  for  the  equi¬ 
librium  constants  Kn",  KD',  and  KD"  are  known.  The  equi¬ 
librium  constants  were  determined  by  thermogravimetric 
measurements,34-33  and  high-temperature-electrical  conduc¬ 
tivity  measurements,29-32  which  give  K„"  3  X  1073  cm-9 

Pa1  2  and  A/f„"  ~~  9-6  X  10*19  J.  The  ionization  energies 
were  determined  using  thermogravimetry  and  Mossbauer 
spectroscopy  by  Hagemann.13  These  results  yielded  values 
for  £f>~  and  £f>-’~  of  0.8  and  2.4  eV,  respectively.  These 
values  inserted  into  Eqs.  (4),  (10),  and  (11)  and  combined 
with  Eqs.  (6)  and  (7)  yield  values  for  (Fe-’*j,  (Fe'*|.  and 
(Fe4*)  as  a  function  of  oxygen  partial  pressure  and  tempera¬ 
ture.  Results  of  such  calculations  for  an  equilibrium  tem¬ 
perature  of  800°C  and  the  total  Fe  concentrations  relevant 
to  this  study  (0.3,  50,  500,  and  1000  ppm)  are  plotted  in  Fig 
6. 

A  number  of  important  trends  are  illustrated  in  Fig.  6. 
For  each  Fe  concentration  labeled  at  the  left  of  the  figure, 
the  three  plateaus  going  from  left  to  right  correspond  to 
saturation  of  the  Fe  ions  in  the  divalent,  trivalent,  and  tetra- 
vaient  states,  respectively,  with  increasing  £<»,.  Note  that 
Fe1*  and  Fe3*  are  the  predominant  species  under  the  rela¬ 
tively  oxidizing  atmospheres  encountered  during  crystal 
growth,  whereas  the  Fe2*  concentration  is  extremely  low 
(CIO13  cm'1)  even  for  the  highest  Fe  concentration  and 
lowest  P<>,  encountered  in  this  study  (1000  ppm  of  Fe.  10~4 
atm).  The  specific  Pu  -dependent  trends  illustrated  in  Fig. 
6  can  be  derived  by  combining  Eqs.  (4),  (7),  and  ( 10)  or  Eqs. 
<41.  <6i,  and  'll)  and  differentiating  with  respect  to  log  P<>. 
Performing  analysis  yields 

d  log[Fe  *|/d  log  P(l  = -d  log[Fe- ‘ |  d  log  P„  =14  (14) 


Fig.  6.  Calculated  variations  in  the  concentrations  of  Fe2*,  Fe1*. 
and  Fe4*  as  a  function  of  oxygen  partial  pressure  at  various  Fe 
concentrations  for  an  equilibration  temperature  of  800°C. 


for  low  values  of  Po„  whereas  at  higher  values  of  Po, 

d  log[Fe4*]/d  log  P0  =1/4  (15) 

and 

d  log(Fe  '*]/d  log  P0  =  -1/6.  (16) 

Similar  manipulation  of  these  equations  with  respect  to  the 
total  Fe  concentration,  |Fe],„,.  at  a  fixed  Pq,  yields 


d  iog(Fe2*]/d  log[Fe]1(„  =  1/2 

(17) 

d  log[Fe4*]/d  log(Fe]„„  =  3/2. 

(18) 

[Equations  ( 14)-l  16)  apply  to  the  regions  in  Fig.  6  where  the 
concentration  of  a  given  valence  state  is  changing  and  break 
down  as  saturation  is  approached  ]  Note  from  Eqs.  (17)  and 
(18)  that  the  concentration  of  a  given  valence  state  does  not 
increase  linearly  with  Fe  content  as  one  might  intuitively 
expect.  It  is  also  notable  that  the  Fe4*  concentration  in¬ 
creases  rather  rapidly  with  Fe  content,  as  compared  with 
Fe2*. 

The  above  defect  model  governs  the  changes  in  oxygen- 
vacancy  concentration  in  addition  to  Fe  valence.  These  are 
plotted  in  Fig.  7  as  a  function  of  oxygen  partial  pressure  for 
the  same  conditions  as  those  in  Fig.  6.  Note  that  the  oxy¬ 
gen-vacancy  concentration  is  fixed  over  a  wide  Po  range  by 
the  Fe1*  concentration  according  to  Eq.  (6).  At  lower  Po 
values,  however,  Eq.  (5)  dominates  Eq.  (7),  and  additional 
oxygen  vacancies  are  charge  compensated  by  conduction 
electrons,  resulting  in  high  electrical  conductivity.  Com¬ 
parison  with  Fig.  6  reveals  that  this  intrinsic  compensation 
mechanism  takes  over  before  substantial  Fe2*  can  be 
formed  At  high  oxygen  partial  pressures,  the  oxygen-va¬ 
cancy  concentration  reaches  a  plateau  determined  by  the 
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level  of  background  impurities,  which  in  Fig.  7  is  determined 
from  the  chemical  analysis  results  in  Table  1.  An  important 
consequence  of  the  trends  illustrated  in  Fig.  7  is  that  oxygen 
vacancies  can  be  ruled  out  as  the  donors  responsible  for  Pu  - 
dependent  changes  in  photorefractive  properties  because 
their  concentration  is  fixed  over  a  wide  P(),  range  by  the 
presence  of  acceptor  impurities  or  dopants. 

Finally,  although  the  above  trends  apply  at  room  tempera¬ 
ture,  some  care  should  be  taken  in  applying  the  actual  defect 
concentrations  shown  in  Fig.  6  to  crystals  at  25°C.  In  the 
case  of  oxygen  vacancies,  rapidly  quenching  the  samples 
from  800°C  effectively  freezes  in  the  high-temperature  ionic 
defect  structure  so  that  the  concentrations  in  Fig.  7  should 
still  hold.  The  relative  concentrations  of  Fe-’*,  Fe1*,  and 
Fe4*,  however,  depend  on  the  degree  of  ionization  of  defects 
in  the  system,  which  is  highly  temperature  dependent.  On 
cooling,  some  deionization  will  occur,  which,  in  the  case  of 
holes,  will  shift  the  curves  in  Fig.  6  to  lower  oxygen  partial 
pressures,  thus  favoring  the  formation  of  higher  oxidation 
states. 

Optical  Absorption  Spectroscopy 

The  variations  in  Fe  valence  governed  by  the  above  model 
were  monitored  by  optical-absorption  measurements  based 
on  the  characteristic  spectra  of  Fe-*,  Fe1*,  and  Fe4*  in 
BaTiOi.  Because  no  published  absorption  spectra  were 
available  for  melt-grown  BaTiO:1:Fe,  results  for  Fe-doped 
SrTi03  (Ref.  35)  served  as  a  guide  to  the  interpretation  of 
our  measured  spectra.  In  SrTiO  :2300  ppm  Fe,  Fe4* 
charge-transfer  bands  have  been  associated  with  an  absorp¬ 
tion  peak  and  shoulder  at  440  and  590  nm,  respectively,  in 
oxidized  samples  (P0,  =  1-10*4  atm),  giving  them  a  reddish- 
brownish  color.35  Under  more  reducing  conditions  (Pq.  = 
10'4-10-17  atm),  Fe  exists  almost  entirely  in  the  trivalent 
state,  leaving  the  same  crystals  nearly  colorless.  Significant 
amounts  of  Fe2*  are  formed  only  in  heavily  reduced  crystals 
(P0,  <  10*1')  and  have  been  associated  with  a  double  band 
with  peaks  at  1030  and  825  nm,  which  cause  the  crystals  to 
appear  greenish  in  color.35 

Optical-absorbance  measurements  of  our  BaTiOcFe  sam¬ 
ples  yielded  similar  results.  These  were  performed  using  a 


Fig.  7  Calculated  variations  in  oxygen-vacancv  concentration  \er 
sus  Po  for  different  Fe  doping  levels  and  background  acceptor 

levels  in  Table  1. 


Fig.  8.  Measured  optical-absorption  spectra  as  a  function  of  Fe 
doping  for  as-grown  crystals  after  mechanical  poling:  (a)  E  parallel 
to  c,  tb)  E  perpendicular  to  c. 

Perkin-Elmer  Lambda  9  double-beam  spectrophotometer. 
Visible  absorption  spectra  were  measured  between  380  and 
860  nm  on  as-grown  and  annealed  samples  for  light  polarized 
parallel  and  perpendicular  to  the  c  axis.  Automatic  back¬ 
ground  correction  compensated  for  the  effects  of  matched 
polarizers  and  a  4  mm  circular  beam  mask.  The  absorption 
spectra  of  the  high-purity  sample  and  of  those  doped  with  5, 
50,  250,  500,  750,  and  1000  ppm  of  Fe  were  measured  in  the 
as-grown  condition  after  mechanical  poling,  and  the  results 
are  shown  in  Fig.  8.  In  each  case  the  absorption  coefficient  a 
was  calculated  from  the  measured  transmission  T  according 
to 

u  =  |2  ln(l  —  /?)  —  In  T\/d,  (19) 

where  d  is  the  sample  thickness  and  R  is  the  reflectance, 
which  is  related  to  the  refractive  index  by 

R  =  in  -  1  )Vin  +  l)2  (20) 

for  the  case  of  normal  incidence.3.  Corrections  for  reflec¬ 
tance  were  made  using  Eq.  (20)  and  the  refractive-index 
data  of  Wemple  et  al.  '*  for  BaTiO  i  as  a  function  of  wave¬ 
length  and  polarization. 

In  Fig.  8  the  effect  of  Fe  doping  on  the  absorption  coeffi¬ 
cient  of  as-grown  crystals  is  clear.  The  data  for  £|lc  and 
£lc  are  given  in  Figs.  8(a)  and  8(b),  respectively.  These 
spectra  are  characterized  by  a  broad  absorption  band,  ex¬ 
tending  throughout  the  visible  to  the  band  edge,  which  in¬ 
creases  with  increasing  Fe  content.  In  addition,  the  extend¬ 
ed  band  edge  appears  to  shift  to  higher  wavelengths  with 
increasing  Fe  concentration.  The  E±c  polarization  shows 
the  emergence  of  a  well-defined  shoulder  at  620  nm  for  Fe 
levels  greater  than  500  ppm.  The  slight  inconsistencies  in 
the  observed  trends  (i.e.,  the  intersection  of  the  5-  and  250- 
ppm  data  with  the  50-  and  500-ppm  curves)  were  due  to 
scattering  from  residual  90°  domain  walls  in  the  5-  and  250- 
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Fig.  9.  Effects  of  oxidation  (Pi.  =  1  atm.800°('i  and  reduction  tP  .  =  l(r*  atm.  800°Cl  on  the  absorption  spectra  of  BaTiO  :  la)  high  puritv, 
(b)  50-ppm  Fe.  (c)  500-ppm  Fe.  (d)  1000-ppm  Fe  i /.  perpendicular  toe). 


ppm  samples,  which  increased  the  measured  ahsorhance. 
These  inconsistencies  were  removed  when  electrical  poling 
was  performed;  however,  not  all  the  samples  shown  survived 
this  procedure. 

A  second  set  of  samples  from  boules  of  high-purity.  50-. 
500-,  and  1000-ppm  Fe  were  measured  after  annealing  at 
800°C  in  oxygen  (1  atm)  and  in  100  ppm  ol  O.'Ar  ( 10' 1  atm), 
respectively.  Oxygen  partial  pressures  were  controlled  by- 
using  premixed  Ar/Oj  gases  and  were  monitored  by  using  a 
calcia-stabilized  zirconia  electrochemical  oxygen  sensor  in 
series  with  the  sample  furnace. 

Before  annealing,  the  c  faces  of  the  samples  were  painted 
with  platinum  paste,  and  the  crystals  were  mounted  between 
platinum  sheet  electrodes  in  adjacent  spring-loaded  sample 
holders.  The  crystals  were  heated  slowly  through  the  Curie 
temperature  (4-6°C/h)  and  then  were  heated  to  SOO°C  in  1 
h.  Samples  were  annealed  for  equilibration  time1-  based  on 
the  bulk  diffusion  data  of  Wernicke  w  ( — 36  h)  and  subse¬ 
quently  quenched  to  just  above  T,  (where  the  oxygen  diffu- 
sivity  is  negligible)  in  order  to  freeze  in  the  high-tempera¬ 
ture  defect  structure.  Quenching  was  performed  by  trans¬ 
ferring  the  gas-sealed  quartz-glass  tube  from  the  spin  tube 
furnace  used  for  annealing  to  an  adjacent  turnace  main¬ 
tained  at  l?5°C.  A  field  of  1000  V/cm  was  gradually  ap¬ 
plied,  and  the  crystals  were  slowly  cooled  through  the  Curie 
temperature  (4-6°C/h)  and  then  to  room  temperature  be¬ 
fore  removing  the  field.  This  procedure  eliminated  the  need 
for  repoling  after  the  high-temperature  anneal. 

The  effects  of  oxidation  (1  atm)  and  reduction  H0_<  atml 
on  pure,  50-,  500-.  and  1000-ppm  Fe-doped  samples  are 
shown  in  Fig.  9.  Oxidation  increased  the  height  of  the  broad 
visible  absorption  band,  turning  the  more  heavily  doped 
crystals  a  aarker  brown  and  the  lightly  doped  crystals  a 
slightly  deeper  yellow  or  amber.  The  effect  of  reduction  was 
to  reduce  the  height  of  the  620-nm  band  drastically,  chang¬ 
ing  the  heavily  doped  crystals  from  brown  to  yellow  while 
leaving  the  high-purity  crystal  almost  colorless. 

An  interpretation  of  Figs,  8  and  9  can  be  based  on  the 


Log  ( iron  concentration  ( ppm)] 


Fig.  10  Log  versus  log|Fel.  ,.  illustrating  the  trend  toward  the 
theoretical  3  2  dependence  at  higher  Fe  concentrations  where  the 
820- nm  shoulder  i  an  be  resolved 


Fig  11  P, ,  dependence  of  the  620  nm  absorption  coefficient 
i  Fe;"  I  approach! n_  the  predicted  1/4  slope  with  increasing  Fe  con¬ 
tent 
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previous  identification  of  charge-transfer  bands  in  SrTiOa 
and  on  the  observed  Fe  concentration  and  Po.-dependent 
trends.  The  broad  absorption  band  in  the  visible  would 
then  correspond  to  an  Fe4’  charge-transfer  band,  with  the 
620-nm  (2.0-eV)  shoulder  corresponding  to  the  590-nm  fea¬ 
ture  in  SrTiOi:Fe.  The  shoulder  in  BaTiO  i  is  shifted  to 
slightly  longer  wavelength  because  of  the  differences  in  the 
crystal  field.  The  analog  to  the  440-nm  peak  in  the  SrTiO  ■ 
spectrum  occurs  at  absorbance:-  out  of  the  range  of  the  spec¬ 
trometer  for  the  sample  thiiKi.tases  used  (3-4.5  mm)  and 
probably  explains  the  apparent  shifts  in  the  band  edge  with 
increasing  Fe  content. 

Qualitatively,  the  interpretation  above  conforms  well  with 
the  proposed  defect  model.  In  the  as-grown  condition  the 
Fe  ions  are  present  in  both  trivalent  and  tetravalent  states, 
and  the  relative  concentration  of  Fe44  increases  with  in¬ 
creasing  Fe  content.  This  manifests  itself  as  a  nonlinear 
increase  in  the  620-nm  absorption  peak  with  increasing  Fe 
content,  as  observed  in  Fig.  8(b).  When  the  crystals  are 
oxidized,  more  of  the  Fe:,+  is  converted  to  Fe4+,  and  vice 
versa  on  reduction,  resulting  in  the  trends  seen  in  Fig.  9.  It 
is  also  significant  to  note  that  no  absorptions  corresponding 
to  Fe2+  (825  nm  in  SrTiO.i)  were  observed  even  for  the  most 
reducing  conditions  investigated  herein.  This  is  consistent 
with  the  defect  model’s  prediction  concerning  the  low  con¬ 
centration  of  Fe2+  in  these  samples. 

Quantitatively,  the  Fe  concentration  dependence  and  oxy¬ 
gen  partial  pressure  dependence  of  the  620-nm  absorption 
deviate  somewhat  from  the  Fe4+  concentration  behavior 
predicted  by  the  defect  model,  but  the  fit  improves  with 
increasing  Fe  content.  Figure  10  shows  that  although  the 
620-nm  absorption  increases  more  gradually  than  predicted 
for  Fe4+  at  low  Fe  concentrations,  the  3/2  slope  predicted  by 
Eq.  (18)  is  closely  followed  in  samples  in  which  the  620-nm 
shoulder  can  be  resolved  (>500  ppm  Fe).  Similarly,  the  Po, 
dependence  of  the  620-nm  absorption  coefficient  plotted  in 
Fig.  11  is  low  in  lightly  doped  samples  but  approaches  the 
1/4  slope  predicted  for  Fe4+  [Eq.  (15)]  in  the  500-  and  1000- 
ppm  Fe-doped  samples. 

PHOTOREFRACTIVE  MEASUREMENTS 

High-purity  and  Fe-doped  BaTiO.j  crystals  in  the  as-grown, 
oxidized,  and  reduced  condition  were  characterized,  using  a 
variety  of  photorefractive  techniques,  and  the  observed 
trends  were  correlated  with  the  variations  in  Fe  concentra¬ 
tion  and  valence  described  above.  In  all  the  experiments 
described  below,  an  argon-ion  laser  operating  in  the  TEMoo 
mode  was  used  to  write  diffraction  gratings  in  the  samples, 
and  a  He-Ne  laser  (A  =  633  nm)  was  used  to  generate  the 
readout  beam  when  required.  All  measurements  were  taken 
at  room  temperature  ( T  -  295  K). 

Saturation  diffraction  efficiency  measurements  provided 
a  preliminary  estimate  of  the  relative  magnitude  of  the  pho¬ 
torefractive  activity  in  the  various  samples.  These  mea¬ 
surements  were  applied  to  high-purity  50-,  500-,  and  1000- 
ppm  Fe-doped  crystals  in  the  as-grown  condition  and  to 
samples  from  the  same  boules  that  were  oxidized  (Po.  =  1 
atm)  and  subsequently  reduced  (Po.  =  10"4  atm).  In  these 
experiments  gratings  with  a  spacing  of  2  pm  (A  =  514.4  nm, 
20  *  14°;  grating  vector  lying  in  the  plane  containing  the  c 
axis)  were  written  in  the  samples,  and  the  intensity  of  the 
diffracted  readout  beam  was  recorded  as  a  function  of  time 


until  saturation  was  reached.  The  results  for  the  500-ppm 
Fe-doped  crystal  are  shown  in  Fig.  12  and  indicate  that  the 
saturation  diffracted-beam  intensity  dropped  by  almost  a 
factor  of  10  when  the  crystal  was  oxidized  while  increasing 
by  a  lesser  degree  when  it  was  reduced.  Similar  results  were 
obtained  for  the  50-  and  1000-ppm  samples. 

Additional  saturation  diffraction  efficiency  measure¬ 
ments  were  performed  on  the  as-grown  and  reduced  samples 
in  which  the  data  were  normalized  with  respect  to  the  inci¬ 
dent  intensity  (1  pW/cm2)  of  the  633-nm  readout  laser  beam. 
The  writing  beams  in  this  case  were  of  wavelength  488  nm, 
with  intensities  l\  =  17  jiW/cm2  and  I>  =  870  jiW/cm:,  and 
were  combined  at  an  external  angle,  28  =  27°,  by  using  the 
anisotropic  configuration  (grating  vector  /C_Lc).  The  dif¬ 
fracted-beam  intensities  were  measured  by  an  EG&G  photo¬ 
diode.  All  samples  were  copolished  to  identical  thicknesses 
of  4.5  mm  ±  3%  to  permit  direct  comparison  of  the  results, 
which  are  shown  in  Fig.  13  as  a  function  of  Fe  concentration. 

Two-beam-coupling  measurements  were  used  to  deter¬ 
mine  the  sign  of  the  dominant  carrier  and  the  effective  trap 
density  for  the  pure  as-grown  and  reduced  crystals.  In  these 
experiments,  the  pump  and  signal  beam  (A  =  488  nm)  were 
ordinary  polarized  with  intensities  of  Ip  -  150  mW/cm2  and 
h  ~  180  MW/cm2,  respectively,  thus  keeping  the  modulation 
index  m  =  2  (lpIs)U2/(Ip  +  Is)  small  (m  «  1)  for  compatibil¬ 
ity  with  the  standard  photorefractive  models.13  40  41  Under 
these  conditions  the  coupling  gain  T  of  the  two  beams  is 
defined  by  the  equation13 


where  Is?  and  Isnp  are  the  intensities  of  the  signal  beam  with 
and  without  the  pump  beam,  respectively;  l  is  the  crystal 


Fig.  12.  Effects  of  oxidation  and  reduction  on  the  diffracted-beam 
intensity  versus  time  for  BaTiO  i:500-ppm  Fe.  Similar  results  were 
achieved  for  50-  and  1000-ppm  Fe-doped  BaTiOt. 
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Fig  13.  Saturation  diffraction  efficiency  of  as-grown  (filled  cir- 
clesl  and  reduced  (open  circles)  BaTiO,  as  a  function  of  Fe  concen¬ 
tration. 
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Fig.  14.  Beam  coupling  gain  as  a  function  of  grating  vector  K  ot 
high-purity  BaTiO ..  The  filled  circles  r,  present  the  as  grown  crys¬ 
tal.  the  open  circles  represent  the  reduced  crystal,  and  the  solid  line 
is  the  theoretical  fit. 


thickness;  and  d  is  the  angle  of  the  beams  inside  the  crystal. 
The  gain  of  Eq.  (21)  can  be  expressed  by1"" 


4tr  bn 

mb  cos(f) ) 


(22) 


where 


bn  =  - 


(23) 


and  where  n..  is  the  ordinary  refractive  index,  r, ,  is  the 
appropriate  electro-optic  coefficient,  and  £,(  is  the  steady- 
state  space-charge  field,  which  has  the  form41 


[<K)mkT 

e 

1  +  (K/K,,)’ 


(24) 


in  which  k  is  Boltzmann's  constant,  c  is  the  electronic 
charge,  and  j '(K)  represents  the  amount  of  electron-hole 
competition.  The  term  K,r  is  related  to  the  effective  trap 
density  \'t  by  the  equation1 V4" 


K„ 


e-.Y-r 

«„kT  ’ 


(25) 


where  t  is  the  relative  dielectric  constant  and  <ois  the  permit¬ 
tivity  of  free  space. 

The  beam-coupling  data  are  plotted  as  a  function  of  grat¬ 
ing  vector  K  in  Fig  14.  From  the  direction  of  beam  coupling 
it  was  determined  that  the  charge  carriers  were  positive  for 
both  the  as-grown  and  reduced  crystals.  Fitting  these  data 
to  the  above  equations,  we  find  that  f ( K '  is  fairly  constant 
over  the  range  of  K  values  used  (8  =  2-40°).  Using  a  value  of 
150  for  the  dielectric  constant  associated  with  this  crystal 
orientation,19  we  find  that  the  effective  trap  density  is  (2.9  ± 
0.4)  X  10"lfi  cm"'  and  (2.7  ±  0.5)  X  10lfi  cm"  for  the  as- 
grown  and  reduced  crystal,  respectively.  Note  that  within 
the  limits  of  error,  reduction  had  no  effect  on  the  measured 
trap  density  of  the  high-purity  crystals. 

Light-induced  grating  erasure  rate  measurements  were 
also  used  to  confirm  the  effective  trap  densities  in  the  pure 
samples.  This  technique  was  used  previously  to  determine 


materials  properties  of  Bi i  SiO  (Ref.  42)  and  BaTiO ( i  Ref. 
13)  and  involves  writing  a  grating  in  the  crystal  by  using  two 
coherent  plane  waves,  removing  the  writing  beams,  and 
Hooding  the  crystal  with  a  uniform  erase  beam.  The  decay 
of  the  grating  diffraction  efficiency  is  simultaneously  moni¬ 
tored  with  a  low-power  readout  beam.  When  the  modula¬ 
tion  index  of  the  fringe  pattern  and  the  absorption  coeffi¬ 
cient  are  small,  the  decay  rate  1/r  is  generally  in  the  form  of  a 
single  exponential  and  is  related  to  the  grating  vector  K 
through  the  equation1 1J" 


1  =  o  ('T,’h  +  'I,') 

T  ( 


1  +  (tkT'S 
1  +  ( L,,K}: 


(26) 


where  <Tph  is  the  photoconductivity;  nt,  is  the  dark  conductiv¬ 
ity;  NY  is  the  trap  density  as  before;  and  the  diffusion  length 
Li,  is  given  by  L/r  =  nTRkT/e ,  where  m  is  the  carrier  mobility 
and  rg  is  the  recombination  time. 

The  above  experiments  were  performed  for  the  case  of 
anisotropic  diffraction.4 1  The  writing  beams  were  of  wave¬ 
length  488  nm,  with  intensities  7,  =  200  mW/cmJ  and  / j  =  0.1 
mW/cm-’.  The  erase-beam  intensity,  It:  -  14  mW/cm-’,  was 
held  constant  throughout  the  experiment  because  decay 
rates  in  BaTiOi  usually  exhibit  a  nonlinear  intensity  depen¬ 
dence.  The  readout-beam  intensity  was  kept  less  than  1 
/jW/cm-.  This  readout  beam  (X  =  633  nm)  was  chopped  at  a 
frequency  of  1000  Hz,  and  the  diffracted  beam  was  synchro¬ 
nously  detected  by  using  a  photomultiplier  tube  and  a  lock- 
in  amplifier.  The  output  of  the  lock-in  amplifier  was  fed  to 
an  AT&T  6300  computer  equipped  with  a  Data  Translation 
acquisition  board  with  the  sampling  rate  set  at  80  Hz.  The 
grating  decay  data  were  least -squares  fitted  to  a  single  expo¬ 
nential,  and  the  decay  rate,  1/r,  was  calculated  from  the  fit. 
The  measurements  were  repeated  as  a  function  of  grating 
vector  by  varying  the  writing-beam  angles. 

The  decay  rate  as  a  function  of  the  square  of  the  grating 
vector  for  the  range  of  K  values  used  is  shown  in  Fig.  15. 
When  a  value  of  3600  is  used  for  the  dielectric  constant 
appropriate  for  this  orientation19— the  theoretical  fit  gives 
trap  density  values  of  (3.1  ±  0.3)  X  10IB  cm"  1  and  (3.3  ±  0.4) 
X  101”  cm"  1  for  the  high-purity  as-grown  and  reduced  crys¬ 
tals.  respectively  These  values  are  in  good  agreement  with 
those  obtained  from  the  beam-coupling  results. 
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Fig  15  Light-induced  grating  decay  rate  as  a  function  of  the 
grating  vector  squared  for  high-purity  BaTiO  in  the  as-grown 
(filled  circles)  and  reduced  (open  circles)  condition  Because  the 
theoretical  fits  overlapped,  only  one  curve  is  shown  (solid  curvet. 
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Intensity  ( W  /cm? ) 

Fig  16.  The  grating  diffraction  efficiency  deoav  time  as  a  function 
of  the  erase-beam  intensity  for  as-grown  and  reduced  BaTiO,:  (a) 
high  purity,  (b)  doped  with  50  and  500  ppm  of  Fe 

The  intensity  dependence  of  the  light  -induced  decay  time 
was  also  measured  for  the  above  high-purity  crystals  as  well 
as  for  the  50-  and  500-ppm  Fe-doped  crystals  in  the  as- 
grown  and  reduced  condition.  These  measurements  were 
conducted  at  a  grating  spacing  of  A  =  2.8  *<m,  and  the  results 
are  shown  in  Figs.  16(a)  and  16(b)  for  the  high-purity  and 
doped  crystals,  respectively.  These  data  show  that,  like  the 
measured  trap  density,  both  the  intensity  dependence  and 
the  absolute  decay  time  remained  the  same  in  the  as-grown 
and  reduced  high-purity  crystals.  Note  also  that  the  addi¬ 
tion  of  large  amounts  of  Fe  decreases  the  speed  of  the  crys¬ 
tals  and  that  this  speed  is  further  reduced  and  the  intensity 
dependence  altered  by  reduction. 

DISCUSSION  AND  CONCLUSIONS 

On  the  basis  of  the  theoretical  defect  model  for  BaTiO  ,:Fe 
and  the  observed  optical-absorption  spectra,  it  is  apparent 
that  Fe  '*  and  Fe4*  are  the  dominant  valence  states  of  Fe 
present  in  our  samples  in  the  Po  range  investigated.  Thus 
it  is  clear  that  the  popular  model'1  in  which  Fe  *  and  Fe-’*  are 
the  sources  and  traps  of  photorefractive  charge  carriers  does 
not  apply  in  this  case  and  that  trivalent  and  tetravaient  Fe 
would  be  expected  as  the  active  photorefractive  centers  if  Fe 
plays  any  role. 

Despite  the  low  levels  of  transition-metal  impurities  in  our 
high-purity  BaTiOi  crystals,  a  relatively  strong  photorefrac¬ 
tive  effect  was  exhibited  by  these  samples  In  fact,  the 
effective  trap  densities  determined  in  the  as-grown  high 
purity  samples  were  of  the  order  of  3  x  10'"  cm  which  is 
comparable  with  those  values  reported  for  nominally  un 
doped  BaTiO  w  ith  2  to  3  orders  of  magnitude  higher  impu 
nty  levels.  Furthermore,  when  high-purity  BaTiO,  sain 
pies  were  reduced,  the  effective  trap  density,  carrier  sign, 
and  speed  were  virtually  unchanged  Because  the  delect 
model  predicts  that  this  degree  of  reduction  should  yield  a 


nearly  sevenfold  increase  in  the  concentration  of  Fe1*,  the 
observed  photorefractive  behavior  of  the  high-purity  crys¬ 
tals  is  inconsistent  with  a  model  in  which  Fe'*  acts  as  the 
trap  site  for  photorefractive  carriers.  It  appears,  therefore, 
that  Fe  is  not  the  dominant  photorefractive  species  in  these 
high-purity  BaTiO.i  crystals.  Instead,  these  results  suggest 
that  the  photorefractive  traps  are  associated  with  a  defect 
whose  concentration  is  independent  of  oxygen  partial  pres¬ 
sure  in  the  Po,  range  0.21-10  4  atm.  Possibilities  for  such 
Po -independent  centers  would  include  barium  vacancies  or 
oxygen  vacancies  (see  Fig.  7),  but  these  conjectures  require 
further  investigation. 

Preliminary  photorefractive  characterization  of  Fe-doped 
BaTiO,  indicated  little  change  in  the  saturation  diffraction 
efficiency  of  the  as-grown  crystals  as  a  function  of  Fe  con¬ 
tent.  Even  the  results  for  the  reduced  crystals,  which  dis¬ 
played  lower  overall  absorption  and  slightly  higher  efficien¬ 
cy,  showed  no  systematic  dependence  on  Fe  concentration. 
The  large  range  of  Fe  concentrations  spanned  by  the  data  in 
Fig.  14,  along  with  the  associated  variations  in  the  Fe4*/Fe  * 
ratio  given  by  Eq.  (18),  offered  no  visible  improvement  in 
properties. 

Although  the  above  results  seem  to  suggest  that  Fe  centers 
are  not  responsible  for  the  photorefractive  effect  in  the 
doped  crystals,  it  is  important  to  note  that  the  behavior  of 
these  samples  was  different  from  that  of  the  high-purity 
crystals.  Figure  16  indicates  tnat  reduction  affected  both 
the  speed  and  the  intensity  dependence  of  the  grating  era¬ 
sure  process  and  that  these  varied  as  a  function  of  Fe  content 
as  well.  In  addition,  the  observed  decreases  and  increases  in 
saturation  diffraction  efficiency,  which  occurred  on  oxida¬ 
tion  and  reduction,  respectively,  are  not  inconsistent  with  a 
model  involving  Fe4*  and  Fe1*  as  the  sources  and  traps  of 
photorefract  i  ve  charge  carriers.  The  interpretation  of  these 
data  is  complicated,  however,  by  a  large  Po  -dependent  ab¬ 
sorption  coefficient  in  the  Fe-doped  crystals.  For  example, 
the  observed  decrease  in  saturation  diffraction  efficiency  on 
reduction  may  simply  be  due  to  increased  absorption  of  the 
writing  beams  by  nonphotorefractive  processes. 

To  complicate  the  issue  further,  we  have  recently  ob¬ 
served  evidence  of  strong  intensity-dependent  absorption  in 
the  Fe-doped  crystals  as  opposed  to  the  high-puritv  sample, 
which  showed  none.  These  observations  have  inhibited  the 
interpretation  of  beam-coupling  and  grating  decay  measure¬ 
ments  on  the  Fe-doped  crystals,  which  will  be  presented 
elsewhere  in  the  near  future.  In  addition,  work  is  continu¬ 
ing  to  expand  the  characterization  of  the  crystals  used  in  this 
study  to  cover  a  wider  oxygen  partial  pressure  range  in  order 
to  clarify  the  role  of  iron  centers  in  the  photorefractive  pro¬ 
cess  in  BaTiO  further. 
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4  Photoinduced  Optical  Absorption  in  BaTiO'3  : 
Fe 


34 


APPLIED  PHYSICS  LETTERS  1991 


Photoinduced  Optical  Absorption  in  BaTiO^Fe 

Doyle  A.  lemple'1*  and  Cardinal  Ward** 

Departments  of  Physics  and  Electrical  Engineering  and  Computer  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge.  Massachusetts  02139 


Measurements  oi  the  steady  state  photoinduced  absorption  Ao  have  been  performed  on 
single  crystalline  samples  ot  pure  and  iron  doped  BaliOj.  The  results  show  that  Ao  has 
a  sublinear  intensity  dependence,  similar  to  that  observed  in  photorefractive  gratimr  decay 
measurements. 
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Determination  of  the  origin  of  the  centers  responsible  for  the  photorefractive  effect  in 
BaTi03  is  complicated  by  the  presence  of  several  competing  effects  including  si  miniated  pho- 
torefractive  scattering  (commonly  known  as  beam  fanning).  1,2  anisotropic  scattering.  1  1 
and  photoinduced  absorption.  10  To  better  aid  in  modelling  of  these  charge  transport  pro¬ 
cesses,  experimental  methods  must  be  employed  that  effectively  de-couple  these  competing 
effects. 

The  first  measurements  of  the  photoinduced  absorption  in  BaTi03  were  performed  by 
Motes  et  al  6~  8  and  Brost  et  al.  9  More  recently.  Pierce  et  al  10  showed  that  absorption  grat¬ 
ings  could  be  produced  in  BaTi03  using  a  configuration  that  does  not  allow  photorefractive 
beam  coupling.  Using  this  method,  these  autors  were  able  to  determined  the  trap  density 
associated  with  the  absorption  gratings. 

Iron  has  long  been  suspected  as  being  the  active  photorefractive  centers  in  BaTi()3. 
Klein  and  Schwartz  11  found  a  good  correlation  between  the  iron  concentration  and  the 
photorefractive  trap  density  in  various  commercial  samples.  Although  these  samples  had 
several  impurities  that  could  be  photorefractive  centers,  the  strongest  correlation  was  with 
the  iron  concentration.  Godefroy  et  al  12  measured  the  diffraction  efficiency  in  unpoled, 
iron-doped  samples  and  found  a  peak  near  750  ppm.  Previously,  the  diffraction  elficieucy 
and  photorefractive  speed  of  iron  doped  BaTi03  samples  were  reported  by  Schunemann  et 
al.  13  More  recently,  we  reported  measurements  of  the  photorefractive  trap  density  in  the  iron 
doped  samples  using  the  light  induced  grating  erasure  technique.  19  In  this  paper,  we  describe 
a  method  of  measuring  photoinduced  absorption  effects  using  a  configuration  which  is  free 
from  photorefractive  beam  coupling,  allowing  accurate  quantification  of  the  photoinduced 
absorption.  Aq. 

The  crystals  used  in  our  experiments  were  grown  from  a  melt  containing  high-purity 
BaC03  and  Ti02  and  extreme  care  was  taken  to  minimize  melt  contamination  by  the  furnace. 
Spark  mass  spectrographic  analysis  of  the  pure  crystals  showed  less  than  0.3  ppm  iron  (see 
Table  for  a  list  of  the  other  impurities  present).  After  growth  of  two  high  purity  boules. 
several  iron  doped  boules  were  grown  from  the  same  melt  using  the  composition  0.34  Vt 
BaC03  and  0.66  %  [(l-x)Ti02  +  (x/2)Fe203]  where  x  is  the  concentration  in  ppm.  In  these 
experiments  50  ppm,  500  ppm  and  1000  ppm  concentrations  were  used.  Chemical  analysis  of 
these  boules  indicates  a  segregation  coefficient  near  unity  (see  fable  ).  At  lea>t  two  samples 
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from  each  of  these  boules  were  cut  along  (100)  planes  and  then  mechanically  and  electrically 
poled.  One  sample  was  used  as  grown  and  the  other  annealed  at  800  °C  in  a  reducing  oxygen 
paitial  pressure  of  10  4  atmospheres.  For  a  more  complete  description  of  the  growth  and 
preparation  of  these  samples  the  reader  is  referred  to  Ref.  13. 

I  wo  problems  associated  with  photoinduced  absorption  measurements  are  wave-mixing 
effects  and  nonuniform  illumination  along  t fie  path  of  the  probe  beam.  Nonuniform  photoin¬ 
duced  absorption  along  the  path  of  the  probe  beam  results  when  both  the  pump  and  probe 
beams  pass  through  the  center  of  the  sample.  This  causes  the  intensity  of  the  pump  beam 
to  change  along  the  path  ot  the  probe  beam.  Calculation  of  Ao  would  require  a  detailed 
knowledge  of  the  analytical  equations  describing  the  photoinduced  absorption,  which  is  the 
unknown  being  determined.  Notiuniloi m  illumination  effects  can  be  minimized  by  two  simple 
changes.  First,  a  small  diameter  probe  beam,  approximately  1  mm.  is  directed  parallel  to. 
and  just  below  the  crystal  surface  to  be  illuminated.  Second,  expanding  the  diameter  of  the 
pump  beam  so  that  it  is  much  larger  than  the  crystal  to  approximate  uniform  illumination. 

I  he  experimental  setup  shown  in  Fig.  1  minimizes  these  problems  in  several  ways.  Wave 
mixing  effects  are  eliminated  by:  (  1  )  using  two  separate  lasers  for  the  pump  and  probe  beams 
and  rj)  intersecting  the  beams  in  tin-  x-y  plane  of  the  erysi.d.  which  by  symmetry  forbids 
beam  coupling  of  the  writing  beams  in  BaliO,.’*  Using  this  setup.  Ao  is  approximately 
constant  along  the  path  of  the  probe  beam.  Finally,  the  intensity  of  the  probe  beam  is  kept 
below  1/AV/cm2  to  minimize  photoinduced  absorption  effects  due  to  the  probe  beam.  Under 
these  conditions  the  photo-induced  change  in  the  absorption,  measured  by  the  probe  beam, 
is  given  by  ' 

.  r  (w  itliiHii  |>imi|i) 

An  =  -p- 1 :» — — - — — - -.  (1) 

[  »  (  Willi  plLUI|>) 

t 

where  L  is  the  thickness  of  the  crystal  parallel  to  the  probe  beam.  [j",,llout i>l,,,ir>)  js  ^ jie  trans¬ 
mitted  probe  beam  intensity  without  the  pump  beam  present,  and  |((""l'piu"|,)  js  the  intensity 
of  the  probe  beam  with  the  pump  beam  present.  For  these  experiments,  the  pump  source 
was  a  3-watt  argon  ion  laser  operating  at  4SS  nm.  The  pump  beam  was  expanded  and  colli¬ 
mated  using  two  lenses  and  passed  through  a  1/2  wave  plate  and  a  polarizer  before  entering 
the  crystal.  The  half  wave  plate  allowed  the  intensity  of  the  pump  beam  to  be  easily  varied 
by  simply  rotating  plate. 
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A  second  100  mVV  air  cooled  argon  ion  laser  operating  at  4SS  uni  was  used  lot  tin-  probe 
beam.  This  beam  was  expanded  and  collimated  to  a  diameter  of  approximately  tin  and 
directed  through  a  polarizer  positioned  in  a  rotator,  followed  by  another  polarizer  aligned 
along  the  crystal  axis,  and  finally  through  a  1  mm  diameter  aperture.  The  intensity  ot  the 
probe  beam  was  modulated  at  a  frequency  of  1000  hL  using  an  optical  beam  chopper,  and 
synchronously  detected  using  a  photomultiplier  tube  and  a  lock-in  amplifier.  In  all  of  these 
measurements  the  pump  and  probe  beam  polarizations  were  the  same,  that  is,  both  parallel 
to  the  a  axis  (ordinary  polarization)  or  c  axis  (extraordinary  polarization).  The  change 
in  the  absorption  was  measured  at  various  pump  beam  intensities  for  both  ordinary  and 
extraordinary  polarization. 

The  results  of  these  measurements  are  shown  in  Figs. 2  and  3.  1  he  change  in  the  ab¬ 
sorption  for  the  pure  as-grown  crystal  is  not  shown  since  there  was  none  detected  at  the 
intensities  used  in  this  work.  The  rise  time  of  the  induced  absorption  was  faster  than  the 
resolution  of  our  detection  system,  which  was  a  few  milliseconds,  while  the  dark  decay  time 
was  on  the  order  of  10-20  minutes.  Because  of  this,  the  measurements  started  at  the  low 
pump  intensities  and  proceeded  to  the  high  intensities.  Alter  a  run.  the  crystal  was  allowed 
to  return  to  the  dark  equilibrium  conditions.  As  shown  in  the  figures,  the  photoinduced 
absorption  Aa  increases  nonlinearly  with  intensity  for  all  the  samples.  However,  since  the 
experiment  was  designed  to  minimize  intensity  variations  across  the  crystal  surface,  very 
high  intensities  were  not  achievable  and  we  were  unable  to  observe  saturation. 

As  shown  in  table  3,  both  the  absorption  and  the  change  in  the  absorption  were  found  to 
increase  with  iron  content  for  a  given  polarization,  pump  beam  intensity,  and  oxidation  state. 
Also,  for  any  particular  dopant  level,  the  absorption  and  the  photoinduced  absorption  <••  the 
as  grown  sample  was  always  laiger  than  that  of  the  reduced  sample.  1  he  sublinear  intensity 
behavior  of  the  photoinduced  absorption  resembles  that  of  the  sublinear  photoconductivity 
measured  by  Ducharme  and  Feinberg.  15  The  erasure  rate  of  our  as-grown  samples  was  pre¬ 
viously  reported  to  increase  sublinearly  with  intensity  H  and  therefore  seem  to  correspond  to 
the  type  B  crystals  discussed  by  Mahgereftreh.  16  Since  the  intensity  dependence  of  the  pho¬ 
toinduced  absorption  and  the  photoconductivity  are  similar,  it  is  likely  that  shallow  level-., 
as  described  in  the  models  of  Mahgereftreh  16  and  Brost.  1 '  play  a  role  in  the  photoinduced 
absorption.  If  the  shallow  acceptor  levels  were  the  source  of  the  photoinduced  absorption. 
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then  annealing  a  type  B  crystal  in  reducing  atmospheres,  should  result  in  partial  Idling  of 
the  shallow  acceptor  levels  and  thus  a  decrease  in  the  photoinduced  absorption.  As  shown  in 
the  figures  and  table  3.  the  crystals  that  were  annealed  in  reduced  oxygen  partial  pressures 
did  show  a  decrease  in  the  photoinduced  absorption. 

In  conclusion,  we  have  performed  photoinduced  absorption  measurements  on  pure  and 
he  doped  samples  of  BaTi03.  Our  results  show  that:  (1)  there  is  a  photoinduced  absorption 
effect  that  increases  with  the  addition  of  iron,  (2)  the  photoinduced  absorption  increase'- 
sublinearly  with  intensity,  and  (3)  annealing  the  crystals  in  reduced  oxygen  partial  pressures 
lowers  the  photoinduced  absorption.  At  present  it  is  not  clear  if  the  charge  centers  responsible 
for  the  photoinduced  absorption  are  due  to  iron,  or  if  the  charge  compensating  defects,  such  as 
oxygen  vacancies  or  Fe-oxygen  vacancy  complexes  are  the  cause.  Further  experiment >.  such 
as  correlation  of  the  EPR  spectra  with  the  photorefract ive  trap  density  and  photoinduced 
absorption,  would  be  needed  to  clarify  the  exact  nature  of  the  photorefractivc  centers. 
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Tables 


lahle  I:  ‘■'park -Source  Mass  Spectiographic  Analysis  of  Cndoped 
BaTi03.  After  Ref.(  13) 


Element 

Level  (ppm)  present 
in  BaTi03 

Bi+ 

5 

Na+ 

<  2 

Mgi+ 

2 

AC+ 

0.5 

Si4  + 

4 

P5+ 

0.05 

S‘“ 

<  2 

Cl- 

10 

K+ 

<0.1 

Ca2+ 

2 

Cr3;2+ 

<0.02 

.\iiru,+ 

<0.02 

Fe*-3-4+ 

0.3 

Nr  !* 
Cu5-2+ 

0.01 

As5+ 

0.1 

Sr+ 

50 

lable  11:  Chemical  analysis  of  the  Ba  l  iC)3  samples.  Alter  Re|.(  13) 
u  Atomic  absorption  (Northern  Analytical). 

6  Spark-source  mass  spectrometry  (Northern  Analytical). 


Iron  Doping  Level 
added  to  the  Melt  (ppm) 

Iron  Level  in  BaTiO.3 
Crystals  Analysed  (ppm)“ 

0.3* 

50  49 

500 

530 

1000 

yj>o 

Table  III:The  absorption  coefficient  at  4SS  nm  was  measured  using  a  Perkin- 
Lambda  9  double  beam  spectrometer.  The  photoinduced  absorption  data  is  given 
laser  intensity  of  0.5  W/cm2.  All  of  the  data  is  for  ordinary  polarization. 


Fe  Concentration 

Absorption 

(cm  *) 

Ao  ( 

cm  1 ) 

As  Grown 

Reduced 

As  Grown 

Reduced 

Pure 

0.53 

0.47 

- 

0.11 

50  ppm  Fe 

0.76 

0.66 

0.7 

0.31 

500  ppm  Fe 

1.73 

0.63 

1.1 

0.0 

1000  ppm  Fe 

- 

1  .'23 

1.05 

1 .05 

1.  Setup  used  to  study  photo-induced  absorption. 


dimer 
at  a 


42 


Photo-Induced  Absorption  Aa  (cm-1)  Photo-Induced  Absorption  AOt  (cm 


Intensity  (W  /  cm2) 


Intensity  (W/cm2) 


2.  The  photo-induced  absorption  in  BaTiC>3  plotted  as  a  function  of  the  pump  beam 
intensity  for  (a)  pure  and  (b)  50  ppm  Fe. 


Photo-Induced  Absorption  Aa  (cnr^ )  Photo-Induced  Absorption  Aa  (cm^) 


10'3  10"2  10'1  10° 
Intensity  (W  /  cm2) 


3.  The  photo-induced  absorption  in  BaTi03  plotted  as  a  function  of  the  pump  beam 
intensity  for  (a)  500  ppm  Fe  and  (b)  1000  ppm  he. 
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We  have  used  the  light-induced  grating  erasure  technique  to  measure  the  photorefractive 
properties  of  pure  and  iron  doped  BaTi03.  Our  experiments  were  performed  using  an 
anisotropic  configuration  which  forbids  beam  coupling  and  self-diffraction  between  the 
writing  beams,  resulting  in  plane  parallel  intensity  fringes  that  do  not  change  with  time.  This 
allowed  direct  measurement  of  the  charge  transport  processes  without  any  feedback  on 
the  grating  due  to  wave-mixing  processes.  The  results  from  these  experiments  show  that  the 
photorefractive  parameters  known  as  the  trap  density  Npt  and  the  mobility-recombination 
time  product  firn  vary  with  intensity  and  are  not  material  constants. 


I.  INTRODUCTION 

The  photorefractive  properties  of  BaTiOj  and  other 
photorefractive  materials,  have  been  studied  by  a  number 
of  researchers  using  a  variety  of  optical  techniques. 1-6 
Many  of  the  experimental  results  are  found  to  be  described 
by  the  hopping  model1  or  the  band  transport  model.7  One 
method  for  optical  characterization  of  photorefractive  ma¬ 
terials  is  the  light-induced  grating  erasure  technique.’  Us¬ 
ing  this  technique,  Feinberg  found  that  the  grating  spacing 
dependence  of  the  photorefractive  gratings  in  BaTi03  was 
well  described  by  the  models,1,7  and  could  be  used  to  de¬ 
termine  the  photorefractive  trap  density  Npr.  Mullin  and 
Hellwarth  performed  similar  experiments  in  Bi12SiO20 
(BSO)  and  showed  that  the  technique  could  also  be  used 
to  determine  the  mobility  recombination  time  product 
firr.5  However,  Ducharme  and  Feinberg  found  that  the 
photoconductivity  of  BaTi03,  determined  from  the  light- 
induced  grating  decay  rates  of  photorefractive  gratings, 
varies  nonlinearly  with  intensity.8  More  recently,  Brost  et 
al .9  found  that  their  beam  coupling  data  showed  that  the 
photorefractive  trap  density  in  BaTi03  varied  with  inten¬ 
sity. 

In  this  paper  we  report  the  results  of  light-induced 
grating  erasure  measurements  on  pure  and  Fe  doped 
BaTi03  using  anisotropic  gratings  that  do  not  allow  beam 
coupling  or  self-diffraction  between  the  writing  beams.  The 
results  show  the  photorefractive  parameters  known  as  the 
trap  density  Npt  and  the  mobility  recombination  time  prod¬ 
uct  fir,  are  functions  of  intensity. 

II.  THEORY 

A.  Photorefractive  response  time 

Consider  the  band  transport  model  of  Kukhtarev7  in 
which  there  are  donor  and  acceptor  levels  in  the  band  gap 
of  the  crystal.  For  this  example,  we  have  chosen  ( 1 )  the 
acceptor  density  to  be  larger  than  the  donor  density, 
NA  >  No.  and  (2)  the  transport  is  dominated  by  hole 
conduction.8  If  all  the  electrons  from  the  donors  arc  trans¬ 
ferred  to  the  acceptors,  the  density  of  unfilled  acceptors  is 


Na  —  These  unfilled  acceptors  are  assumed  to  be  avail¬ 
able  to  be  redistributed  by  the  light  among  all  the  acceptor 
levels.  If  the  applied  intensity  pattern  is  in  the  form  of  a 
sinusoidal  fringe  pattern 

/(jc)=/0[1  +  m  cos(ATx)],  (1) 

where  m<  1,  and  K  is  the  spatial  frequency  of  the  fringes, 
the  resulting  charge  distribution  may  also  be  sinusoidal. 
For  this  case,  the  charge  transport  solutions  of  Kukhtarev7 
and  Feinberg1  both  predict  an  exponential  response  for  the 
photorefractive  gratings  given  by 

1  gph/l-f  (ekT/SN^Kh 

r  e  \  1  +  {T4ikbT/e)K2  )’  U) 

where 


en[(s/hv)/o  +  0r}(JVA-ND) 


— 


yAd 


(3) 


is  the  photoconductivity,  is  the  mobility,  yr  is  the  recom¬ 
bination  constant,  v  is  the  light  frequency,  s  is  the  photo¬ 
ionization  cross  section,  h  is  Planck’s  constant,  0r  is  the 
thermal  generation  constant,  e  is  the  static  dielectric  con¬ 
stant,  kb  is  the  Boltzmann  constant,  T  is  the  lattice  tem¬ 
perature, 


rr=l/(YAD), 

is  the  recombination  time  and 


(4) 


Np1=ND  (5) 

is  the  photorefractive  trap  density.  The  only  intensity  de¬ 
pendence  in  Eq.  (2)  is  the  photoconductivity,  given  by  Eq. 

(3) ,  which  is  predicted  to  be  linear  in  intensity.  Also,  Eqs. 

(4)  and  (5)  give  the  mobility  recombination  time  product 
and  the  photorefractive  trap  density  as  material  constants. 


B.  Anisotropic  diffraction 

For  the  experiments  reported  here  anisotropic  gratings 
were  used.  That  is  the  grating  vector  lies  in  the  x-y  plane  of 
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FIG  1  Beam  geometry  illustrating  (a)  writing  and  (b)  readout  of  an 
anisotropic  photorefracttve  grating  in  the  x-y  plane  In  tat  the  writing 
beams  A,  and  A:  intersect  In  the  crystal  producing  a  photorefractise 
grating  with  a  wave  vector  K  =  k;  -  k,  In  (b)  an  extraordinary  polar 
ized  readout  beam  R  scatters  into  an  ordinarx  polarized  beam  S.  where 
K  =  ks  -  k„ 


the  crystal.1011  One  of  the  advantages  of  anisotropic  dif¬ 
fraction  is  that  there  are  certain  configurations  in  which 
gratings  can  be  produced  with  no  beam  coupling  or  self- 
diffraction  allowed  between  the  writing  beams  i:  For  two 
waves  intersecting  in  a  photorefractive  crystal  the  amount 
of  coupling  or  self-diffraction  between  the  w  riting  beams  is 
proportional  to  the  effective  susceptibility.1' 

*dT=(?r32)[?N  -  ?0-EJ-  t(, )•?;].  (6) 

where  e  ,  and  e  2  are  the  polarization  vectors  of  the  writ¬ 
ing  beams,  e  R  and  e  s  are  the  polarization  vectors  of  the 
readout  and  scattered  beams  respectively,  ?  0  is  the  second 
rank  optical  dielectric  tensor,  R  is  the  third  rank  electro- 
optic  tensor,  and  E  is  the  space  charge  field.  If  the  two 
waves  interfere  in  the  x-y  plane  of  4  mm  symmetry  BaTiO> 
as  shown  in  Fig.  1(a),  the  effective  susceptibility  is  zero 
when  the  polarizations  of  the  writing  beams  are  both  ordi¬ 
nary  or  both  extraordinary.1'  In  this  case,  the  feedback 
between  the  charge  transport  processes  and  the  wave  dif¬ 
fraction  processes  is  eliminated,  resulting  in  intensity 
fringes  that  are  parallel  planes  and  that  do  not  change  with 
time.  Thus  for  very  small  modulation,  m<l,  the  solutions 
derived  by  Kukhtarev7  and  Feinberg.1  are  very  close  to  the 
actual  experimental  conditions  Since  the  feedback  from 
the  wave-mixing  processes  on  the  charge  transport  pro¬ 
cesses  of  grating  formation  has  been  elinimated.  each  can 
be  studied  independently.  This  is  therefore  an  excellent 
configuration  for  the  measurement  of  photorefractive  pa¬ 
rameters. 

Equation  (6)  does  allow  anisotropic  diffraction  to  oc¬ 
cur  for  readout  of  this  grating.  In  this  instance,  the  effective 
susceptibility  is 

Yrir=(  ?*■?:)[£/»•(  -  e0-E-^-  ?<>)• 

where  e  R  and  e  5  are  the  polarization  vectors  of  the  read¬ 
out  and  scattered  beams  For  a  negative  uniaxial  crystal  the 
phase  matching  angles  for  anisotropic  diffraction  as  shown 
in  Fig.  1(b)  are  given  by  " 


FIG.  2.  Optical  system  used  for  light-induced  grating  erasure  as  a  func¬ 
tion  of  grating  vector  K. 


"o  -  nl 

fiiO 

Ar 

_ _  c:_  /  n\ 

4  sin(0) 

Ar  i 

—  “  sin  v  tf ) 
Aw 

and 

sin(as)  =  [«o  -  n]  +  simfa*)  ]  1/2,  (9) 

where  aR  and  as  are  the  angles  outside  the  crystal  of  the 
extraordinary  polarized  readout  beam  and  ordinary  polar¬ 
ized  scattered  beam  measured  from  the  surface  normal,  n0 
and  ne  are  the  indices  of  refraction  at  the  readout  beam 
wavelength,  6  is  the  angle  of  the  writing  beams  with  the 
surface  normal,  and  A  w  and  kR  are  the  writing  and  readout 
beam  wavelengths. 

III.  EXPERIMENTAL  METHODS 
A.  Light-induced  grating  erasure 

The  light-induced  grating  erasure  technique  involves 
writing  a  photorefractive  grating  in  the  crystal  with  two 
coherent  plane  waves,  then  removing  the  writing  beams 
and  flooding  the  crystal  uniformly  with  an  erase  beam.1  4 
The  grating  decay  rate  is  simultaneously  monitored  with  a 
low  power  (I  /iW/cm;)  readout  beam  incident  at  the 
Bragg  angle 

The  optical  setup  used  for  the  light-induced  decay  rate 
experiment  is  shown  in  Fig.  2.  In  this  setup  an  argon  ion 
laser  was  used  for  the  writing  beams  and  the  erase  beam. 
As  shown  in  Fig  2.  a  beam  splitter  labeled  BS1.  was  placed 
just  after  the  laser  to  produce  two  separate  beams  that  were 
expanded  and  collimated  using  spatial  filters.  The  beam 
reflected  by  BS1,  was  used  as  the  erase  beam  This  beam 
was  directed  onto  mirror  Ml  then  through  the  center  of 
the  crystal  mount  with  the  line  between  Ml  and  the  crystal 
mount  defining  the  system  axis.  The  intensity  of  the  erase 
beam  was  varied  using  a  rotating  half  wave  plate  followed 
by  a  polarizer.  The  beam  transmitted  through  the  beam 
splitter  BS1  was  directed  through  the  second  beam  splitter 
BS  to  produce  two  writing  beams.  These  beams  were  di¬ 
rected  symmetrically  about  the  system  axis.  Optical  table 
rails  were  mounted  parallel  to  the  two  writing  beam  paths 
with  mirrors  mounted  on  each  rail  to  redirect  the  beams 
onto  the  crystal.  The  rails  were  calibrated  to  allow  accurate 
angular  measurements.  The  intensity  of  one  of  the  w  riting 
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beams  was  varied  using  rotating  polarizer  followed  by  a 
fixed  polarizer.  This  was  necessary  to  ensure  that  the  mod¬ 
ulation  index  was  small  m<l. 

For  BaTi03,  phase  matching  for  anisotropic  diffraction 
is  limited  to  an  angular  range  of  3°<a<20°  for  488  nm 
writing  beams  and  633  nm  readout  beams,  where  a  is  the 
incident  angle  of  the  readout  beam  outside  of  the  crystal. 
We  note  here  that  extreme  care  must  be  taken  so  as  not  to 
confuse  higher  diffracted  orders  which  are  easily  observed 
in  the  anisotropic  configuration  in  BaTi03.14  Also,  we 
found  that  small  shifts  in  the  position  of  the  erase  beam  in 
the  crystal  caused  significant  variations  in  the  decay  rates. 
Therefore,  the  crystal  and  the  erase  beam  were  not  moved 
during  the  experiments. 

A  HeNe  laser  was  used  to  readout  the  gratings  (see 
Fig.  2).  This  beam  was  expanded  using  two  lenses  and  the 
intensity  was  varied  using  a  rotating  polarizer  followed  by 
a  fixed  polarizer.  The  beam  was  directed  along  the  optical 
rail  by  two  mirrors  and  modulated  by  a  beam  chopper  at  a 
frequency  of  1000  Hz.  The  readout  beam  angle  could  be 
varied  by  moving  a  mirror  along  the  rail  and  redirecting 
the  beam  to  the  crystal.  The  diffracted  beam  was  synchro¬ 
nously  detected  using  a  photomultiplier  tube  and  a  lock-in 
amplifier  with  the  output  of  the  lock-in  captured  by  an 
AT&T  6300  computer  equipped  with  a  Data  Translation 
data  acquistion  board.  The  sampling  rate  was  set  at  80  Hz 
and  the  grating  decay  rates  were  found  to  be  well  fitted  by 
a  single  exponential.  Also,  a  633  nm  interference  filter  was 
placed  in  front  of  the  photomultiplier  tube  to  eliminate 
bleaching  from  the  argon  laser. 

A  typical  measurement  was  made  as  follows:  First 
shutter  SI  was  opened  and  any  gratings  present  in  the 
crystal  were  erased.  Shutter  SI  was  closed  and  shutter  S2 
opened  allowing  the  writing  beams  to  intersect  in  the  crys¬ 
tal  forming  a  photorefractive  grating.  This  grating  was 
monitored  by  measuring  the  diffraction  of  a  HeNe  beam 
incident  upon  the  crystal  at  the  Bragg  angle.  Eq.  (8). 
When  the  diffraction  efficiency  reached  saturation,  shutter 
S2  was  closed  and  simultaneously  SI  opened,  erasing  the 
grating.  After  three  or  four  decay  times  were  measured,  the 
angle.s  of  the  writing  beams  were  changed  and  the  mea¬ 
surement  repeated  After  one  complete  scan  of  the  angular 
range  was  finished,  the  erase  beam  intensity  was  increased 
and  the  measurement  repeated.  The  writing  beams,  and  the 
erase  beam  were  ordinary  polarized  for  all  the  measure¬ 
ments  reported  here. 

B.  Sample  preparation 

The  crystals  used  in  our  experiments  were  grown  from 
a  melt  containing  high  purity  starting  materials  and  ex¬ 
treme  care  was  taken  to  minimize  melt  contamination  by 
the  furnace.  After  growth  of  two  high  purity  boules,  sev¬ 
eral  iron  doped  boules  were  grown  from  this  melt  at  the 
concentrations  50,  500,  and  1000  ppm  in  the  melt.  Samples 
from  each  of  these  boules  were  cut  along  (100)  planes  and 
then  mechanically  and  electrically  poled  One  sample  from 
each  doping  concentration  was  annealed  at  800  °C  in  a 
reducing  oxygen  partial  pressure  or  10 4  atm.  Another  sam¬ 
ple  was  prepared  from  a  melt  containing  pure  TiO;  and 


FIG  3  Intensity  dependence  of  (A)  the  trap  densitv  and  (B)  the  mo¬ 
bility  recombination  time  product  in  as  grown  BaTiO,  The  dashed  lines 
are  qualitative  guides. 

reagent  grade  BaO  This  sample  will  be  referred  to  as  nom¬ 
inally  pure  since  it  was  grown  from  the  type  of  feed  mate¬ 
rials  commonly  used  in  the  growth  of  BaTiO,.  For  a  more 
complete  description  of  the  growth  and  preparation  of 
these  samples  the  reader  is  referred  to  Ref.  15.  Beam  cou¬ 
pling  measurements  indicate  that  all  the  crystals  exhibited 
p  type  conduction.1'' 

IV.  RESULTS  AND  DISCUSSION 

The  decay  rate  versus  grating  vector  data  was  fitted  to 
Eq.  (2)  for  each  erase  beam  intensity,  and  the  parameters 
known  as  the  photoconductivity  <7ptl,  the  photorefractive 
trap  density  ,Vpr,  and  mobility  recombination  time  product 
fir,  were  calculated  from  these  fits.  The  results  are  plotted 
in  Figs.  3,  4,  and  5.  The  decay  rate  of  the  1000  ppm  crys¬ 
tals  showed  nonexponential  behavior  and  was  not  suitable 
for  this  type  of  analysis. 

The  plots  shown  in  Figs.  3  and  4  indicate  that,  for  a 
given  crystal,  the  trap  density  and  the  inverse  of  the 
mobility-recombination  time  product  have  very  similar  in¬ 
tensity  dependencies.  Since  it  is  unlikely  that  the  mobility 
is  intensity  dependent  this  would  seem  to  indicate  that  the 
trap  density  is  related  to  the  recombination  time  as  in  Eqs. 
(4)  and  (5)  This  behavior  was  found  for  all  the  samples 
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FIG  4.  Intensity  dependence  of  (A)  the  trap  density  and  IB)  the  mo¬ 
bility  recombination  time  product  in  reduced  BaTiO,  The  dashed  lines 
are  qualitative  guides. 


measured  including  the  high  purity  sample  and  the  nomi¬ 
nally  pure  sample.  The  nonlinear  intensity  dependence  of 
the  photoconductivity,  shown  in  Fig.  5,  was  found  to  be 
similar  to  that  observed  by  Ducharme  and  Feinberg.2  8 

Intensity  dependent  trap  densities  can  be  explained  by 
the  normal  models  through  trap  depletion.  However,  it  is 
not  unreasonable  to  assume  that  the  intensity  dependent 
recombination  time  and  trap  density  are  in  some  way  re¬ 
lated  to  the  nonlinear  photoconductivity.  Since  the  trans¬ 
port  seems  to  be  well  fitted  by  a  diffusion  type  photocon¬ 
ductivity,  as  evidenced  by  the  grating  spacing  dependence, 
modification  of  the  model  to  include  the  intensity  depen¬ 
dent  effects  may  be  possible.  One  possibility  is  a  model, 
proposed  by  Rose,  in  which  there  exists  a  distribution  of 
levels  in  the  band  gap.1  However,  direct  of  Rose’s  model 
does  not  seem  appropriate  since  it  disagrees  with  temper¬ 
ature  dependencies  reported  by  Ducharme  and  Feinberg.8 

In  conclusion,  we  have  used  the  light-induced  grating 
erasure  technique  to  measure  the  intensity  dependence  of 
the  photorefractive  properties  of  pure  and  Fe  doped 
BaTi03.  Our  results  show  that  direct  interpretation  of  the 
photorefractive  models  would  lead  to  an  intensity  depen¬ 
dent  trap  density  and  recombination  time  with  a  strong 
correlation  between  the  two  Further  experiments,  such  as 
accurate  measurements  of  the  temperature  dependence 
over  a  wider  range  of  intensities,  are  needed  before  a  sub¬ 
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FIG  5  Intensity  dependence  of  the  photoconductivity  in  (A)  as  grown 
BaTiO,  crystals  and  I  B)  reduced  BaTiO,  crystals.  The  dashed  lines  are 
qualitative  guides 


stantial  modification  of  the  photorefractive  charge  trans¬ 
port  models  can  be  accomplished. 
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A  well-defined  ring  of  anisntropicallv  scattered  light  is  observed  when  a  linearly  polarized  laser  beam  is  incident 
upon  a  birefringent  photorefractive  crystal  such  as  BaTiO^,  LiNbOn.  and  LiTaOa.  An  analysis  is  presented  that 
accurately  predicts  the  cone  angle,  polarization,  and  the  location  of  the  discontinuities  in  these  rings.  Unlike  the 
photovoltaic  model  proposed  recently  bv  others,  the  analysis  presented  herein  is  based  on  the  standard  photorefrac¬ 
tive  theory  and  is  independent  <>l  the  photovoltaic  activity  of  the  crystal.  The  analysis  also  illustrates  the  strict 
phase  matching  and  electro-optic  tensor  symmetry  constraints  that  must  be  satisfied  in  order  to  observe  beam 
coupling  and  normal  and  anisotropic  self  diffraction.  Preliminary  experimental  results  are  presented  for  BaTiO, 
and  LiNbOi  and  are  found  to  he  in  excellent  agreement  with  standard  photorefractive  theory. 


INTRODUCTION 

Anisotropic  scattering  of  a  linear,  ordinary  polarized  laser 
beam  into  a  conical  ring  of  extraordinary  polarized  light  was 
recently  observed  in  LiTaOpCu.1  A  theory  based  on  the 
photovoltaic  effect  has  been  advanced  to  explain  this  phe¬ 
nomenon.  It  is  claimed  that  the  ring  is  a  result  of  scattering 
from  noisy  dynamic  gratings  that  are  formed  by  photogal- 
vanic  currents  in  this  material.1 

In  contrast,  we  show  that  standard  photorefractive  theory 
can  explain  the  cone  angle,  polarization,  and  the  location  of 
the  discontinuities  in  the  anisotropicallv  scattered  diverging 
rings  of  light  that  have  been  observed  in  BaTiO  t  and  doped 
LiNbOi  and  LiTaOa. 

Standard  photorefractive  theory  shows  that  (a)  the  con¬ 
figuration  most  often  used  for  wave  mixing  in  photorefrac¬ 
tive  materials  (which  we  shall  refer  to  as  the  standard  config¬ 
uration)  does  not  give  rise  to  anisotropic  diffraction  and  that 
lb)  the  interaction  geometry  that  gives  rise  to  anisotropic 
gratings  in  birefringent  crystals  is  such  that  the  c  axis  does 
not  lie  in  the  plane  containing  the  grating  vector  K  and  the 
propagation  vector  It/?  of  the  readout  beam.  We  shall  refer 
to  the  specific  configuration  in  which  the  c  axis  is  perpendic¬ 
ular  to  K  and  k«  as  the  anisotropic  configuration. 

In  the  anisotropic  configuration,  we  further  show  ( 1 )  that 
beam  coupling  and  normal  self-diffraction  are  forbidden  in 
some  crystals  for  certain  polarizations  of  the  writing  beams 
and  (2)  that  anisotropic  self-diffraction  of  the  writing  beams 
can  occur  only  for  specific  grating  periods  and  specific  write 
light  polariza*'  >n. 

Finally,  we  report  a  set  of  preliminary  experiments  on 
BaTiO i  and  LiNbO.gFe  that  confirm  that  the  rings  are  ade¬ 
quately  described  by  the  standard  photorefractive  theory. 
The  standard  model  was  previously  employed  to  explain 
anisotropic  diffraction  from  holographic  gratings  in  LiNb- 
O.gFe  (Refs.  2-4)  and  BaTiO,, 7  but  the  conditions  that  per¬ 
mit  grating  formation  without  beam  coupling  and  the  obser¬ 
vation  of  anisotropic  scattering  into  ring  segments  were  not 
reported. 


THEORY 

Consistent  with  standard  photorefractive  theory,  we  consid¬ 
er  the  holographic  recording  system  shown  in  Fig.  1  in  which 
two  plane  waves  <S,(r,  t)  and  fi^tr.  t)  with  wave  vectors  k, 
and  k_>  intersect  in  a  photorefractive  crystal.  The  steady- 
state  plane-wave  fields  are  of  the  form 

£ir.  t)  =  eA(r)  expfitk  •  r  -  w/)|  +  c.c.,  (1) 

where  A.  u.  and  e  are  the  wave  amplitude,  frequency,  and 
polarization  vector,  respectively,  of  the  beams.  The  intensi¬ 
ty  pattern  of  the  interfering  beams  induces  a  space-charge 
electric  field  E«c  by  the  photorefractive  effect  with  the  same 
spatial  variation  as  the  fringe  pattern.  That  is.  ignoring  the 
spatial  dc  bias,  we  have 

A‘  A 

E„  =  Eii',*  •  »■•_.)  j  ’  exp|i(K  •  r  +  0(|  +  c.c.,  (2) 

where  r  is  the  position  coordinate.  E  =  Ek  is  the  dc  compo¬ 
nent  of  the  space-charge  field,'’  7  /,>  =  |.4||-  +  |AJ-,  K  =  k>  — 
k|  is  the  grating  vector  of  the  fringe  pattern,  and  <t>  is  the 
phase  shift  between  the  space-charge  field  and  the  intensity 
fringe  pattern  The  magnitude  of  </>  depends  on  the  specific 
charge-migration  processes  in  the  crystal.67  The  change  in 
the  susceptibility  induced  by  this  space-charge  field  is  given 
by6 

Ay  =  — <  •  r  •  E.(  •  i,  (H) 

where  r  is  the  third-rank  electro-optic  tensor  and  t  is  the 
second-rank  optical  dielectric  tensor.  Using  Eqs.  (2)  and 
t:t).  the  change  in  the  susceptibility  can  he  written  as 

A*  A , 

=  i)\ic,  *  •  r ,)  —  '  exp|t(K  •  r  +  0)]  +  c.c..  14) 

*  U 

where 

n  x  =  —  <  -  r  ■  E  •  i  -  I  ->  I 

Now  consider  readout  ot  tin-  ph.i^e  grating  bv  a  plane 
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Kit;.  I  Setup  tor  recording  phase  holograms  in  a  photorefractive 
crystal. 

wave  of  frequency  u>.  Using  the  slowly  varying  field  ap¬ 
proximation/'  we  can  derive  the  following  coupled-wave 
equations  for  the  incident  and  diffracted  readout  beams  (see 
Fig.  2): 


where  R  and  -S’  are  the  wave  amplitudes  of  the  incident  and 
the  diffracted  readout  beams,  respectively.  Here,  the  cou¬ 
pling  coefficients  s*  and  *s  are  given  by 

’’’Xrff  ,rX,tt 

„  = - - ,  ks.  =  - - - -  (  ‘  ) 

ijwX/(COSo^  rjsAs  COSo  s 

where  A#  =  A.„-  is  the  wavelength  of  the  readout  and  the 
diffracted  beams.  17*.  ijs.  o«,  and  o.s  are  the  indices  of  refrac¬ 
tion  and  the  internal  angles  of  the  incident  and  the  diffract¬ 
ed  readout  beams,  respectively,  and  the  effective  suscepti¬ 
bility  \,t|  is  given  by 

x,u  =  <('i  •  •  Ax  •  fs>-  I81 

where  iq  and  e:  are  the  polarization  vectors  of  the  writing 
beams  and  c«  and  i>s  are  the  polarization  vectors  of  the 
incident  and  diffracted  readout  beams,  respectively.  We 
have  also  assumed  that  the  phase-matching  condition 

K  =  k„-k.  (!" 

is  satisfied. 

Therefore,  in  order  for  diffraction  to  be  observed,  general¬ 
ly  two  conditions  must  he  simultaneously  satisfied:  la) 
phase  matching  |Kq.  I9)|  and  lb)  Xrti  must  be  nonzero  |Eq. 
IM|  For  the  special  case  of  anisotropic  diffraction,  the 
'■vrtimelry  of  the '»\  tensor  must  be  such  that  x.01  is  nonzero 
when  and  i\  are  orthogonal,  and  the  refractive  indices 
experienced  by  k/,-  and  k.  must  permit  the  phase-matching 
condition  to  be  ~atislied  It  should  be  noted  that  beam 
coupling  and  -ell  ditlr,.ction  are  also  described  by  Kq.  (S) 
w  it  It  I  be  -ubstitution  R  -  l.N  =  .4 1  ,  t:  =  —  c)t  and  X;,-  - 

V 11  =  V  lot  In-  cum-  we  get 


X,ff  =  (<',*  ■  i'.Hv-,  ■  ix  •!>,•).  (10) 

Based  on  this  special  form  of  x*rr.  it  must  be  emphasized  that 
with  (C|  *  •  eg)  *  0  it  becomes  possible  to  write  an  anisot  ropic 
grating  that  exhibits  no  isotropic  self-diffraction  or  two- 
beam  coupling  between  the  writing  beams  (cj  •  A x  •  £1*  =  0). 
Specific  examples  are  given  in  the  section  that  follows. 


SYMMETRY  PROPERTIES  OF  THE  COUPLING 
COEFFICIENT 

Standard  Configuration 

In  the  standard  configuration,  the  c  axis  is  coplanar  with 
both  the  grating  vector  and  the  wave  vector  of  the  readout 
beam.  That  is,  the  c  axis  is  parallel  to  the  x  axis  of  the 
coordinate  system,  as  shown  in  Fig.  2(a).  The  nonzero  ele¬ 
ments  of  the  electro-optic  tensor  can  be  found  in  Ref.  9. 
Using  these  in  Eq.  (5),  we  find  that  Ax  for  the  standard 
configuration  has  the  same  form  for  the  4-mm.  2-mm,  and  3- 
m  crystal  classes.  That  is, 

(*Xn  Ax)2  0  \ 

Ax2|  AX22  0  I*  (11) 

0  0  4X.r»/ 

Since  Axi  1  =  Ax:u  =  Ax2.i  =  Ax:i2  =  0,  it  follows  from  Eq.  (8) 
that  ordinary  waves  (z  polarized)  cannot  couple  to  extraordi¬ 
nary  waves  lx-}  polarized)  and  vice  versa  (i.e.,  anisotropic 
diffraction  is  forbidden).  Since  the  diagonal  elements  of 
this  tensor  are  nonzero,  conventional-polarization  isotropic 
beam  coupling  is  allowed. 

Anisotropic  Configuration 

In  the  anisotropic  configuration,  the  c  axis  is  perpendicular 
to  the  plane  containing  the  grating  vector  and  the  wave 
vector  of  the  readout  beam.  That  is,  the  c  axis  is  oriented 
along  the  z  axis  of  the  coordinate  system,  as  in  Fig.  2(b).  In 
this  case  A\  has  t  he  following  forms: 

For  3-m  symnu-t  ry 


Fig.  2.  Oriental  inn  <  if  the  taxis  relative  to  the  coordinate  system  of 
fig.  1  Inr  tilt  the  standard  configuration  and  lb'  the  anisotropic 
1 ,  niliguratmn. 
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Fig.  3.  (a)  Wave- vector  diagram  for  readout  of  an  anisotropic  grating  of  wave  vector  K  in  a  negative  birefringent  crystal  (b)  Wave -vector  dia¬ 

gram  illustrating  anisotropic  self-diffraction  of  extraordinary  polarized  write  beams  k„  into  ordinary  polarized  beams  k,  (of  - 1  and  +  2  orders) 
in  a  negative  birefringent  crystal. 


From  the  phase-matching  condition  |Eq.  (9)|  it  follows 
that  diffraction  of  an  extraordinary  readout  probe  beam  of 
wavelength  A«  into  an  ordinary  beam  is  allowed  in  this  crys- 
(13)  tal  when 


. (  0 

n\  =  — En„  n,.- 1  0 

Vsi  cos  > 

where  n,,  and  n(.  are  the  ordinary  and  extraordinary  indices 
of  refraction,  respectively,  and  y  is  the  angle  that  the  grating 
vector  K  makes  with  the  x  axis.  Since  in  each  of  these  cases 
o x :i i •  and  f>x-u  are  nonzero,  it  follows  that  for  all 
three  crystal  classes  anisotropic  diffraction  is  allowed  when 
the  phase-matching  condition  (Eq.  (9)|  is  satisfied.  That  is, 
diffraction  of  an  ordinary  wave  (x-y-polarized)  into  an  ex¬ 
traordinary  wave  (z-polarized)  and  vice  versa  is  possible. 
However,  since  the  elements  along  the  main  diagonal  of  the 
tensor  are  identically  zero  in  the  4-mm  and  2-mm  crystals, 
coupling  between  the  writing  beams  (ordinary  or  extraordi¬ 
nary  polarizations)  is  forbidden.  Additionally,  since  f>xx\  is 
the  only  diagonal  element  that  is  zero  in  3-m  crystals  jEq. 
( 1 2 ) | ,  only  coupling  between  extraordinary  writing  beams  is 
forbidden  in  this  case. 

Phase  Matching 

To  illustrate  phase  matching  in  an  anisotropic  photorefrac- 
tive  grating,  we  consider  a  grating  formed  in  a  negative 
birefringent  f n,,  >  nf  )  uniaxial  crystal  using  the  anisotropic 
configuration  of  Fig.  2(b).  For  gratings  written  with  two 
beams,  each  of  which  makes  an  external  with  the  y  axis, 
the  grating  vector  is  given  by 


0  r  M  cos  y  \ 

0  r , ,  sin  y  j  >  ( 14) 
r4,  sin  y  0  / 


where  o'«  is  (he  external  angle  of  incidence  of  the  readout 
probe  beam  measured  from  the  y  axis  and  ijn  and  are  the 
indices  of  retraction  of  the  crystal  at  A«.  Additionally,  the 
external  angle  .d*  of  the  diffracted  ordinary  beam  is  given  by 

sin  «'s  =  |r)„2  -  -  sin2  <»'fl|l/2.  (17) 

This  phase-matching  condition  is  also  illustrated  by  the 
readout  wave-vector  diagram  shown  in  Fig.  3(a).  From  reci¬ 
procity.  an  ordinary  polarized  probe  beam  incident  upon  the 
grating  at  an  angle  n\  can  diffract  into  an  extraordinary 
beam  at  an  external  angle  o'«.  To  describe  positive  birefrin¬ 
gent  crystals,  one  simply  interchanges  q„  and  nt:  in  Eqs.  ( 16) 
and  1 17). 

Equation  ( 1(5)  can  also  be  used  to  describe  anisotropic  self¬ 
diffraction  of  an  extraordinary  beam  into  an  ordinary  beam 
and  vice  versa  With  the  substitution  A ^  =  An  and  o«  =  II, 
the  external  Bragg  angle  for  anisotropic  self-dilfraction  into 
a  —  1  order  and  a  +2  order  becomes 


118) 


This  indicates  that  only  extraordinary  beams  can  experience 
anisotropic  sell -diffraction  in  negative  birefringent  crystals 
by  the  standard  photorefract ive  effect  This  fact  is  con¬ 
firmed  by  the  wave-vector  diagram  of  Fig.  3(b)  that  illus¬ 
trates  how  the  phase-matching  condition  is  satisfied  for  an 
isotropic  sell-diltraction  of  extraordinary  polarized  beams 
in  a  negative  uniaxial  crystal. 


where  An  is  the  wavelength  of  the  beams. 
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Fig.  4  Anisotropic  scattering  of  an  extraordinary  polarized  beam 
incident  normally  upon  a  BaTiOi  crystal.  Strong  forward  scatter¬ 
ing  into  a  ring  of  ordinary  polarized  light  with  discontinuities  along 
the  c  axis  is  illustrated  in  the  photograph. 


EXPERIMENTAL  RESULTS 


Anisotropic  Gratings 

Anisotropic  gratings  were  produced  in  single  crystals  of  Ba- 
TiO  i  and  doped  LiNb03  by  using  the  configuration  of  Fig. 
2(b).  Both  writing  beams  of  wavelength  488  nm  had  an 
intensity  of  30  mVV/cm-  and  were  incident  symmetrically 
about  the  a  axis  of  each  crystal  at  an  angle  H. 

Normal  coupling  between  the  writing  beams  was  not  ob¬ 
served  in  BaTiO.i  for  either  ordinary  or  extraordinary  polar¬ 
ization  over  a  wide  range  of  incident  angles.  Examination  of 
Eqs.  (8)  and  (13)  for  the  anisotropic  configuration  shows 
that  for  BaTiO.i  (4-mm  symmetry)  x«.ff  is  zero  when  ep  and  e.s 
are  both  ordinary  or  both  extraordinary,  and  hence  beam 
coupling  is  indeed  forbidden. 

In  negative  birefringent  BaTiO:i,  the  indices  of  refraction 
at  633  nm  are  n,,  =  2.36  and  n„  =  2.42. 10  For  a  grating 
written  with  488-nm  light  at  an  incident  angle  of  20°  (grating 
spacing  A  =  0.7  pm),  Eqs.  ( 17)  and  ( 18)  predict  a  Bragg  angle 
of  17.2°  and  24.7°  for  extraordinary  and  ordinary  633-nm- 
wavelength  readout  beams.  The  corresponding  measured 
Bragg  angles  were  17.5°  and  24.5°  for  extraordinary  and 
ordinary  beams,  respectively.  The  intensity  of  the  633-nm 
readout  laser  beam  was  kept  at  0.7  mW/cm-  for  all  polariza¬ 
tions  used.  Both  the  ordinary  and  the  extraordinary  read¬ 
out  beams  were  observed  to  have  the  same  diffraction  effi¬ 
ciency  (■'■  20rr  I. 

As  shown  by  the  theory,  anisotropic  self-diffraction  of 


extraordinary  beams  is  allowed  in  BaTiOi  when  the  phase¬ 
matching  condition  [Eq.  (9))  is  satisfied  (see  Eqs.  (8),  (13), 
and  (18)  and  Fig.  3(b)].  For  488-nm  writing  beams  (nf  = 
2.45  and  n„  =  2.53), 10  the  Bragg  angle  calculated  by  using  Eq. 

( 18)  for  anisotropic  self-diffraction  into  a  -1  and  a  +2  order 
is  12.5°.  Experimentally,  both  writing  beams  were  observed 
to  undergo  anisotropic  self-diffraction  at  an  input  Bragg 
angle  of  12.5°.  with  a  diffraction  efficiency  of  approximately 
25 *>t  in  the  -1  order  and  less  than  1%  in  the  +2  order.  From 
Eq.  (18),  it  follows  that  there  would  be  no  anisotropic  self¬ 
diffraction  in  BaTi03  for  ordinary  polarized  writing  beams, 
and  none  was  observed.  Similar  results  for  anisotropic  self¬ 
diffraction  into  the  —1  order  were  reported  recently  by 
Kukhtarev  et  at.h 

We  have  also  observed  anisotropic  self-diffraction  when 
one  of  the  writing  beams  was  an  ordinary  wave  and  the  other 
was  an  extraordinary  wave.  In  this  case  only  the  input 
extraordinary  beam  was  observed  to  scatter  weakly  into  a  —  1 
order.  Since  the  Bragg  angle  was  the  same,  12.5°,  and  we 
measured  a  1%  depolarization  in  the  transmitted  beams,  we 
conclude  that  both  writing  beams  experience  depolarization 
in  the  crystal,  allowing  them  to  interfere  and  write  an  aniso¬ 
tropic  photorefractive  grating. 

The  Rings 

When  a  single  extraordinary  (c-polarized)  laser  beam  was 
incident  upon  the  front  surface  of  the  BaTi03  and  LiNb- 
O.iiFe  crystals,  a  well-defined  diverging  ring  of  scattered 
light  polarized  perpendicular  to  the  c  axis  (ordinary  beam) 
was  observed,  in  each  case,  leaving  the  back  face  of  the 
crystal.  The  intensity  in  the  ring  was  found  to  be  maximum 
in  the  direction  perpendicular  to  the  c  axis  and  zero  along 
the  c  axis,  as  shown  in  the  photograph  of  Fig.  4  for  the 
BaTiO.i  crystal.  The  ring  generated  by  the  LiNb03:Fe  crys¬ 
tal  was  much  weaker  than  that  produced  by  the  BaTi03 
crystal.  We  believe  that  the  rings  are  formed  by  diffraction 
of  the  incident  beam  from  noisy  anisotropic  gratings  that 
satisfy  the  phase-matching  condition  for  diffraction  into  — 1- 
order  beams.  These  gratings  are  a  subset  of  those  produced 
through  interference  between  light  scattered  from  optical 
imperfections  in  the  crystal  and  the  incident  beam. 

For  simplicity,  consider  this  process  in  the  a-b  plane  of  the 
crystal  |Fig.  5(a>|.  The  incident  beam  k,  scatters  in  many 


(a)  (b) 


Fig  ’’  W  .ne  t  error  diagrams  1 1  Ins  I  rat  my  how  the  rings  are  lormed  by  scattering  in  BaT  i()t  and  LiNb03.  <  a  1  Shows  how  a  typical  grating 
\<  i  t..r  m  the  n  h  plane  i-  written  b\  interlerence  between  the  incident  extraordinary  polarized  beam  and  a  randomly  scattered  beam.  ( b I 
Sht.w  -  one  "l  i  he  spent  it  grating  sector-  in  the  o-h  plane  that  allows  anisotropic  self-diffraction  of  the  incident  beam  into  a  - 1 -order  beam  on 
i  he  ring 
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directions,  but  only  two  of  the  scattered  beams  k',  lie  at  an 
angle  2t)%d  such  that  the  grating  vector  K  allows  anisotropic 
self-diffraction  of  the  incident  beam  into  a  —  1  -order  beam  ks 
(Fig.  5(b)).  Following  an  analysis  similar  to  that  of  Eq.  (18). 
the  external  apex  angle  of  the  ring  in  the  a-b  plane  is  given 
by 


.  ,,  r  ,  (3ne-  +  "./-'Yl'/ 

L  V  /J 


(19) 


For  BaTiO:t  this  equation  leads  to  an  apex  angle  2t (,  of  50°  for 
normally  incident  488-nm  light.  This  is  in  excellent  agree¬ 
ment  with  the  measured  apex  angle  of  50.2°.  Similarly,  in 
LiNbO  >:Fe  the  calculated  and  measured  angles  of  the  ring  in 
the  a-b  plane  are  52°  and  51°.  respectively.  In  general,  this 
type  of  phase  matching  is  allowed  in  a  ring  around  the 
incident  beam.  The  discontinuity  observed  in  the  ring  cor¬ 
responds  to  the  case  in  which  the  c  axis  lies  in  the  plane 
containing  the  grating  vector  and  the  incident  beam  (stan¬ 
dard  configuration).  This  is  consistent  with  the  photore- 
fractive  theory  that  shows  that  anisotropic  diffraction  is 
forbidden  in  both  LiNbO  i  and  BaTiO;,  for  this  case. 

A  similar  anisotropically  scattered  ring  has  been  reported 
by  Odulov  et  al.1  in  a  LiTaO-gCu  crystal  (positive  birefrin¬ 
gence)  with  an  ordinary  polarized  input  beam.  These  au¬ 
thors  assume  that  the  ring  was  produced  by  noisy  gratings 
formed  by  the  photovoltaic  effect  from  interference  between 
orthogonally  polarized  light  beams.  Their  photovoltaic 
model  predicts  an  apex  angle  of  15°  in  LiTaO,  (n„  =  2.248 
and  nE  -  2.252  at  X  =  440  nm),  which  appears  to  be  in 
relatively  good  agreement  with  their  measured  value  of  13°. 
However,  the  photovoltaic  theory  does  not  predict  the  ob¬ 
served  discontinuity  in  the  ring  along  the  c  axis,  and,  when 
applied  to  BaTiO,,  it  predicts  an  apex  angle  of  75°,  which 
does  not  agree  with  the  observed  value  of  50.2°. 

If,  however,  we  assume  that  the  ring  observed  in  LiTaOa 
was  produced  by  an  anisotropic  grating  through  the  normal 
photorefractive  effect,  Eq.  (19)  would  predict  an  apex  angle 
of  11°.  which  is  also  in  good  agreement  with  the  measured 
value. 


CONCLUDING  REMARKS 

Using  standard  photorefractive  theory,  we  have  described 
the  symmetry  and  phase-matching  conditions  that  must  be 
satisfied  in  order  to  write  anisotropic  gratings  in  birefringent 
photorefractive  crystals.  We  have  identified  the  specific 
configurations  for  writing  anisotropic  gratings  that  do  not 
permit  coupling  between  the  writing  beams.  Such  gratings 
were  written  in  BaTiO,,  and  indeed  no  beam  coupling  was 
observed.  However,  these  gratings  could  be  read  out  bv  a 
laser  beam  of  a  different  wavelength. 


Anisotropic  self-diffraction  was  shown  to  occur  for  only 
one  polarization  of  the  writing  beams  and  at  only  one  grating 
spacing.  For  BaTiO  ,,  LiNbO;, :Fe,  and  LiTaO:i,  anisotropic 
self-diffraction  into  —1  and  +2  orders  is  possible,  and  we 
have  observed  both  orders  in  BaTiO;,. 

A  ring  of  anisotropically  scattered  light  has  been  observed 
when  each  of  these  materials  was  illuminated  with  a  single 
laser  beam.  This  ring  was  shown  to  be  formed  by  diffraction 
from  noisy  anisotropic  gratings  that  satisfy  the  phase¬ 
matching  condition  for  self-diffraction  into  a  -1  order.  A 
previous  model,1  which  was  based  on  gratings  produced  by 
the  photovoltaic  effect,  does  not  consistently  predict  the 
observed  apex  angle  and  the  discontinuities  in  the  ring  of 
light  scattered  from  these  materials.  In  contrast,  the 
photorefractive  analysis  presented  herein  consistently  pre¬ 
dicts  both  the  apex  angle  and  the  discontinuities  in  the  ring 
for  each  of  the  above  materials. 
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the  tensor  photovoltaic  elicit.  anisotropic  scattering  along  the  c  axis  in  am  crystals  is  not  forbidden  by  symmetry 


In  an  earlier  paper.  Odoulov  et  a/.1  reported  on  the  observa¬ 
tion  of  a  single  ring  of  anisotropically  scattered  light  with  an 
apex  half-angle  of  6.5°  fora  LiTaO,  crystal  that  was  pumped 
with  a  single  ordinary  polarized  beam.  A  tensor  photovolta¬ 
ic  model,  which  predicted  an  angle  of  7.5°.  was  advanced  in 
that  paper.  In  our  previous  paper-  we  proposed  a  photore¬ 
fractive  three-wave  mixing  interaction  that  explained  the 
rings  that  we  had  observed  in  BaTiO;  and  LiNhO,.  and  it 
predicted  an  angle  of  5.5°  for  Li'I'aOi.  Because  (1)  the 
measured  cone  angle  in  LiTaO;  was  equally  close  to  the 
calculated  photovoltaic  and  photorefractive  values,  and  (2) 
the  photovoltaic  model  does  not  accurately  predict  the  apex 
angle  of  the  rings  that  we  observed  in  BaTiO  t  and  LiNhO  i  or 
explain  the  nulls  in  the  rings  along  the  c  axis,  we  concluded 
that  the  authors  had  not  sufficiently  proved  that  the  single 
ring  that  they  had  observed  in  LiTaO  t  was  d  ie  to  the  photo¬ 
voltaic  effect. 

In  light  of  their  new  data,  which  now  show  two  rings,  we 
agree  that  the  ring  that  they  originally  ohserved  may  be  due 
to  the  photovoltaic  effect.  However,  our  statement  that  the 
photovoltaic  model,  as  presented  in  their  paper,  does  not 
forbid  anisotropic  scattering  along  the  c  axis  is  correct,  as  is 
shown  below  . 

A  normal  photorefractive  grating  is  produced  by  a  spatial¬ 
ly  varying  intensity  pattern  that  induces  charge  redistribu¬ 
tion  through  diffusion.  Because  of  the  nature  of  this  pro¬ 
cess,  the  space-charge  field  is  directed  along  the  grating 
vector.  Using  this  constraint,  it  was  shown  in  Ref.  2  that 
anisotropic  diffraction  along  the  c  axis  of  3m  crystals  is 
forbidden  by  symmetry  for  normal  photorefractive  gratings. 

On  the  other  hand,  photovoltaic-induced  gratings  are  pro¬ 
duced  by  an  asymmetry  in  the  photoionization  and  recombi¬ 
nation  processes,  which  are  in  turn  related  to  the  spatial 
dependence  of  the  ellipticity  of  the  polarization  of  the  opti¬ 
cal  field  within  the  crystal.  Consider  two  waves  (ordinary 
and  extraordinary ).  intersecting  in  the  plane  containing  the 
<  axis  (see  Fig.  1 ) 

Let  the  ordinary-polarized  pump  wave 

E.  =  .t  .4  ex p [/ 1  k  •  r  —  >)  +  c.c.  ( 1  > 

and  the  extraordinary  scattered  wave 


intersect  in  the  crystal  at  an  angle  »  in  the  plane  containing 
the  r  axis,  as  shown.  According  to  Refs.  1.  4.  and  5.  the 
associated  photovoltaic  currents  are  given  by 

■l  = 

where  d,.*  are  the  photovoltaic  tensor  elements  and  E,  and 
E,,  are  the  electric-field  amplitudes  of  the  light  waves. 

Based  on  the  symmetry  of  the  dpt  tensor  for  3m  crystals 
(independent  symmetric  elements  are  J'vl.  J'LU,  dw.  anc* 
d, „  and  the  antisymmetric  elements  are  ^:!1  and  we 

see  that  spatially  oscillating  currents  along  the  grating  vec¬ 
tor  K  =  k.  -  k  are  forbidden  for  the  geometry  shown  in  Fig. 
(.  However,  current  J)  perpendicular  to  the  grating  vector 
is  not  forbidden,  as  can  be  seen  when  Eqs.  (1)  and  (2)  are 
substituted  into  Eq.  (3).  Since  the  resulting  space-charge 
field  is  given  by  E,"  =  J,/o„  where  n,  is  the  photoconductivi¬ 
ty.  the  photovoltaic  model  predicts  a  space-charge  field  that 
is  directed  along  the  x  axis  and  spatially  modulated  in  the  y- 
c  plane  of  Fig.  1.  For  3m  crystals,  such  a  space-charge  field 
would  induce  a  change  in  susceptibility  through  the  electro¬ 
optic  effect  that  is  given  by 

_i\  =  —  <  lr;mnE,," (41 


and  h,i>  the-  lorm 


E  =  .4  I  r  cos  it  —  \  sin  ftlexp[/ 1  k,  ■  r  -  ^ /  lj  +  i  c.  ( 2 1  Fig.  1  I  In  Mattering  configuration  under  om>iderat i, m 
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Fig  2.  Anisotropic  scattering  ol  a  single  evtraordinary-pularized 
pump  beam  into  two  ordinary  polarized  rings  (with  nullsi  in  a 
heax  tlv  doped  BaTiO  ,:Fe  crystal  1 1  (Min  parts  in  1 1 >’  Fe  in  melt  i 


The  symmetry  of  Eq.  (51  shows  that  coupling  between  or¬ 
thogonally  polarized  beams  in  the  y-g  plane  is  allowed.  The 
strength  of  this  coupling  is.  of  course,  dependent  on  the 
photoconductivity  and  the  magnitude  ot  the  photovoltaic 
and  electro-optic  tensor  elements  involved,  but  it  is  not 
forbidden  by  symmetry.  Therefore,  unlike  in  the  standard 
photorefractive  case,  the  photovoltaic  space-charge  fields 
for  the  configuration  in  Fig.  1  are  perpendicular  to  the  grat¬ 


ing  vector  K  and  can  lead  to  anisot  ropic  dil  tract  ion  along  the 
(  axis. 

In  addition  to  the  single  ring  lapex  halt-angle  ot  2n.l°i 
that  we  reported  earlier  in  Ref.  2  lor  an  nominally  undoped 
BaTiO  .  crystal,  we  recently  observed  a  second  antsotropical- 
Iv  scattered  ring  from  a  heavily  doped  BaTiO  Ke  crystal 
(lotto  parts  in  10''  Fe  in  melt).  This  second  ring  is  much 
weaker  than  the  first  and  contains  two  pairs  of  nulls  located 
along  the  t  and  a  axes,  as  shown  in  Fig.  2.  This  second  ring 
satisfies  the  phase-matching  equation 

2k;,‘  =  k,,"  +  k,  lb) 

with  a  measured  apex  half-angle  of  27. r>°. 

The  authors  are  also  with  the  Center  for  Materials  Science 
and  Engineering,  Massachusetts  Institute  of  Technology. 
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High-order  anisotropic  diffraction  and  anisotropic  self  diffraction  in  barium  titanate  (BaTiOj).  strontium  barium 
niobate  (Sr.Bai-jNb.iOs,  or  SBN),  and  barium  strontium  potassium  sodium  niobate  (Ba;-,Sr,Ki-v 
Na^NbsOis,  or  BSKNN)  have  been  observed:  up  to  fourth  order  in  BaTiO;t  and  SBN  and  up  to  fifth  order  in 
BSKNN.  It  is  shown  that  the  scattered  beams  arising  from  anisotropic  self-diffraction  are  composed  of  multiple 
orders.  To  model  anisotropic  diffraction,  high-order  terms  have  been  included  in  Kukhtarev’s  solution  of  the 
transport  equations  for  diffusion-dominated  transport.  Light-induced  grating  decay  rates  and  diffraction  efficien¬ 
cies  were  measured  as  a  function  of  grating  vector  and  modulation  index  for  the  first  three  orders  in  BaTiC>3  and 
were  found  to  be  in  good  agreement  with  the  model. 


1.  INTRODUCTION 

Anisotropic  photorefractive  gratings'  2  are  characterized  by 
the  property  that,  under  the  proper  phase-matching  condi¬ 
tions,  an  ordinary  polarized  readout  beam  can  diffract  into 
an  extraordinary  polarized  beam  and  vice  versa.  Holo¬ 
graphic  recording  of  anisotropic  photorefractive  gratings  in 
electro-optic  crystals  has  been  studied  extensively  in  recent 
years.1-9  In  the  simplest  case  this  involves  gratings  whose 
wave  vectors  lie  in  the  plane  perpendicular  to  the  c  axis. 

The  diffraction  efficiencies  and  rise  times  of  these  gratings 
were  used  to  determine  the  photoconductivities  and  diffu¬ 
sion  lengths  of  these  materials.  Measurement  of  the  phase¬ 
matching  angles  was  used  to  determine  the  birefringence.3 
The  signal-processing  applications  include  such  novel  tech¬ 
niques  as  dynamic  light  deflection,'1  spatial  light  modula¬ 
tion,5  wave-front  inversion,6"  and,  more  recently,  amplified 
reflection  through  anisotropic  four-wave  mixing.6  Also,  an¬ 
isotropic  scattering  into  conical  rings  of  light  was  shown  to 
be  the  result  of  anisotropic  gratings  written  through  inter¬ 
ference  between  the  incident  pump  beam  and  beams  scat¬ 
tered  in  the  crystal  owing  to  optical  inhomogeneities. 19  The 
angular  dependence  of  this  cone  was  used  to  determine  the 
dispersion  of  the  birefringence  in  KNbOy.3  In  general,  the 
polarization  orthogonality  of  the  readout  and  diffracted 
beams  results  in  an  improved  signal-to-noise  ratio,  as  point¬ 
ed  out  by  Stepanov  et  at.1 

Anisotropic  diffraction  into  a  second  order  was  observed 
previously9  in  BaTi03  by  using  the  configuration  for  aniso¬ 
tropic  self-diffraction.29  In  this  case  extraordinary  polar¬ 
ized  writing  beams  intersect  in  the  x-y  plane  of  the  crystal  to 
form  an  anisotropic  photorefractive  grating.  These  writing 
beams  simultaneously  diffract  into  both  a  first  and  a  second 
order.  Modeling  in  this  instance  is  complicated  by  the  cou¬ 
pling  of  the  two  writing  beams  to  the  four  diffracted  beams. 

In  this  paper  we  investigate  high-order  anisotropic  dif¬ 
fraction  in  barium  titanate  (BaTiOi),  strontium  barium  nio¬ 
bate  (Sri-JBaINb.)0(;,  or  SBN),  and  barium  strontium  po¬ 
tassium  sodium  niobate  <Ba2-,Sr.rK|->NavNb‘,Ois,  or 
BSKNN),  and  we  report  the  first  observation,  to  our  knowl¬ 
edge,  of  anisotropic  self-diffraction  into  orders  higher  than 
the  second  in  these  materials.  We  also  show  that  the  aniso- 


tropically  scattered  beams  arising  from  anisotropic  self-dif¬ 
fraction  are  composed  of  multiple  orders.  To  model  high- 
order  anisotropic  diffraction,  we  use  a  perturbation  ap¬ 
proach,  similar  to  that  used  by  Refregier  et  a/.,10  for 
including  high-order  terms  in  the  model  of  Kukhtarev  et 
a/."  Using  this  model,  we  derive  (1)  approximate  solutions 
for  the  grating  vector  dependence  of  the  second-  and  third- 
order  light-induced  grating  decay  rates  and  (2)  the  modula¬ 
tion  index  dependence  of  the  diffraction  efficiency  of  these 
orders.  We  have  measured  both  of  these  characteristics  for 
the  first  three  orders  in  BaTiOj  and  have  found  reasonable 
fits  to  the  model.  We  also  report  our  measurements  on 
anisotropic  self-diffraction  into  orders  as  high  as  the  fourth 
order  in  BaTiO  and  SBN  and  the  fifth  order  in  BSKNN. 

2.  ANISOTROPIC  DIFFRACTION 

To  review  the  coupling  permitted  for  anisotropic  gratings  in 
4-mm  crystals,  consider  the  grating  formed  by  the  intersec¬ 
tion  of  two  coherent  ordinary  polarized  plane  waves  with 
amplitudes  .4!  and  4.  and  wave  vectors  ki  and  kj  incident 
symmetrically  about  the  y  axis  of  a  4-mm  crystal,  as  shown  in 
Fig.  1(a).  In  this  case  the  grating  vector  K  =  kj  -  k;  is  along 
the  x  axis  of  the  crystal.  The  resulting  interference  fringes 
induce  a  sinusoidal  distribution  in  the  photogenerated  carri¬ 
ers,  which  set  up  a  space-charge  field  with  the  same  spatial 
periodicity  as  the  fringe  pattern.  This  space  charge  then 
modulates  the  refractive  index  through  the  electro-optic  ef¬ 
fect. 

Now  consider  readout  of  this  grating,  as  in  Fig.  1(b),  by  a 
steady-state  plane-wave  field  of  the  form 

A(r,  t)  =  e4(r)exp[ifk  •  r  -  u.f)j  +  c.c.,  (1) 

where  ?  and  uj  are  the  polarization  vector  and  the  frequency, 
respectively.  The  coupled-wave  equations  governing  dif¬ 
fraction  are 


—  ~  IK  rS, 
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Fig.  1.  Schematic  diagrams  illustrating  (a)  writing  and  ( b>  readout 
of  an  anisotropic  photorefractive  grating  in  a  4-mm  crystal.  In  (a) 
the  writing  beams  Ai  and  Aj  intersect  in  the  crystal,  producing  a 
photorefractive  grating  with  a  wave  vector  K  =  k2  -  It!.  In  < hi  an 
ordinary  polarized  readout  beam  R  scatters  into  an  extraordinary 
polarized  beam  S,  where  k,  »  kg  +  K. 


Fig.  2.  Phase-matching  diagrams  showing  (a)  anisotropic  first- 
order  readout  of  a  grating  of  wave  vector  K  and  (b)  anisotropic 
third-order  readout  of  the  same  grating. 


where  R  and  S  are  the  wave  amplitudes  of  the  readout  and 
diffracted  beams,  respectively.  The  coupling  coefficients 
are  given  by 


"X.ft 

- 1 

nr  A,  cos(ar,n) 

(4) 

'X.ft 

- - -  | 

n,  A,  co8(a,‘n) 

(5) 

As  shown  previously,9  for  this  crystal  orientation,  ordinary- 
polarized  writing  beams  cannot  self-diffract  from  such  a 
grating  (£]*  •  5x  •  =  0).  In  this  case,  the  intensity  fringe 

pattern  in  the  crystal  is  uniform,  assuming  small  absorption. 
Without  this  coupling  the  assumptions  in  the  charge-trans¬ 
port  models  of  Kukhtarev  et  al.n  and  Feinberg  et  alu  are 
more  valid  than  for  the  case  of  strong  beam  coupling.  Equa¬ 
tions  ( 2)— (7)  also  show  that  an  ordinary  polarized  readout 
beam  can  diffract  only  into  an  extraordinary  polarized  beam 
for  gratings  produced  as  in  Fig.  1(a). 

The  phase-matching  condition  required  for  observing  an¬ 
isotropic  diffraction  into  high  orders  is  illustrated  in  Fig.  2. 
As  shown  in  Fig.  2(a),  the  readout  beam  can  undergo  first- 
order  diffraction  into  an  extraordinary  polarized  beam  when 
the  phase-matching  condition  K  =  k3  -  k,  is  satisfied.  How¬ 
ever,  as  illustrated  in  Fig.  2(b),  for  a  given  grating  vector  K 
the  phase-matching  condition  may  be  satisfied  for  higher 
orders  of  K.  This  condition  can  be  expressed  as  AfK,  =  K,  - 
K,  or  as 


sin(o) 


M  V  ~  nf2) 
M\r  sin  8 


\,M  sin  d 

A*  " 


(8) 


where  a  is  the  angle  of  the  extraordinary  beam  outside  the 
crystal;  n„  and  n,  are  the  ordinary  and  extraordinary  indices 
of  refraction,  respectively,  at  the  readout-beam  wavelength; 
A w  and  Ar  are  the  writing-  and  readout-beam  wavelengths, 
respectively;  M  is  the  diffracted  order;  and  6  is  the  writing- 
beam  angle  outside  the  crystal. 

Because  ordinary  polarized  writing  beams  in  4-mm  crys¬ 
tals  do  not  experience  beam  coupling  for  grating  vectors  in 
the  x-y  plane,  these  gratings  fit  closely  the  model  of  Kogel- 
nik.13  Therefore  the  diffraction  efficiency  of  Afth  order  is 
given  by 

Sid) 

V  Ri  0) 

where 

IDP  - 

nrnsV  toi-  ar  cos  as 

and 

&X  =  ~i-  (?•  Ew)  •  l  (11) 

where  Ey  is  the  Afth-order  space-charge  field,  d  is  the  crystal 
thickness,  nr  and  n,  are  the  refractive  indices  of  the  readout 
and  scattered  beams,  respectively,  or  and  a,  are  the  angles  of 
propagation  in  the  crystal,  lo~  1 1  +  h-  and  Ar  is  the  readout- 
beam  wavelength. 


=  sin2!f)ld. 


(9) 


where  Ar  *  A,  is  the  wavelength  of  the  readout  and  diffracted 
beams  and  rw,  n„  ar1D,  and  a,1"  are  the  indices  of  refraction 
(n’s)  and  angles  (a’s)  inside  the  crystal  of  the  readout  and 
diffracted  beams,  respectively.  The  effective  susceptibility 
x.ft  is 

X«fr  **•*,).  (6) 

where 

ix  -  Ek  •  «.  (7) 

f  is  the  third-rank  electro-optic  tensor,  i  is  the  second-rank 
optical  dielectric  tensor,  and  EK  is  the  space -charge  field. 


3.  LIGHT-INDUCED  GRATING  DYNAMICS 


The  equations  used  to  describe  the  photorefractive 
of  a  material  are  given  by" 

response 

dn.  dNn+  i 

—  =  — -  +  -  VJ, 

dt  dt  e 

(12) 

dNn*  t*  \ 

if-  =  Uv  +  TND~iV 

(13) 

J  =  en’i  £5C  +  eDfn, 

(14) 
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V(«„£5C)  =  <?(AV  ~  *V4  ~  n),  (15) 

where  nr ,  No,  NA,  and  Np*  are  the  concentrations  of  the  free 
electrons,  donors,  acceptors,  and  ionized  donors,  respective¬ 
ly,  e  is  the  charge  on  an  electron,  J  is  the  current  density,  a  is 
the  photoconductivity,  /  is  the  intensity,  D  is  the  diffusion 
constant,  is  the  recombination  rate,  3  is  the  thermal 
generation  rate,  <  is  the  relative  dielectric  constant,  u  is  the 
mobility,  and  Eit  is  the  space-charge  field.  For  anisotropic 
gratings  in  the  x-y  plane  of  4-ram  crystals,  the  fringes  are 
flat  and  planar  for  ordinary  polarized  writing  beams,  and  the 
fringe  intensity  distribution  is  given  by 

1  =  1,,  1  +  ™  exp hKx)  +  c.c.  .  (16) 

Hence  the  generation  rate  at  any  point  in  the  grating  is 
constant,  assuming  no  charge  depletion.  We  assume  solu¬ 
tions  of  the  form 


Nr/  =  Nm,+ 

+  Vo)  A"  [j  j +  exp(iKx)  +  c.c.). 

(17) 

n  =  n0  + 

‘Afn,  exp(iKx)  +  c.c.], 

(18) 

J  =  Ju  + 

%[J\  expdKx 1  +  c.c.). 

(19) 

E,c  =  £n  + 

'/'_,[£,  exp(iA'x)  +  c.c.). 

(20) 

where  Npi*,  and  E\  are  the  modulated  portion  of  the 

ionized  donor  density,  free-electron  density,  current  densi¬ 
ty,  and  space-charge  field,  respectively,  and  m  =  2fct  ■ 
^2>(/ l/g)1  "/(/ 1  +  h).  Then  the  decay  rate  of  the  first-order 
space-charge  field  is  given  by"-1-’ 


where  NT  =  (A’q/A’.A)(A'o  -  SA). 

To  solve  for  the  high-order  harmonics  of  the  grating  we 
can  use  a  perturbation  approach10  by  assuming  high-order 
solutions  of  the  form 

A '[,*  =  N[„*  +  ‘4IAV  exp(tA’x)  +  •  •  • 

+  A 'pm*  exp(Mihx)  +  c.c.),  (22) 

n  =  nu  +  V2 1 ri |  expUh'x)  +  •  •  •  +  nAf  exp(AfiKx)  +  c.c), 

(23) 

J  =  J„  +  %\.JX  exp(iKx)  +  ••  •  +  expt MiKx)  +  c.c). 

(24) 

£sc  =  £ri  +  V 2 ( E ,  exp(iKx)  +  •  •  •  +  EAf  exp(.\/i Kx)  +  c.c.) 

(25) 

For  a  second-order  perturbation,  we  ( 1 )  neglect  orders  in  the 
harmonic  solutions  higher  than  the  second  order;  (2)  assume 
that  the  second-order  space-charge  field  does  not  signifi¬ 
cantly  at.ect  the  first-order  space-charge  field,  that  is,  the 
first-order  solution  is  still  valid;  and  (3)  neglect  any  gratings 
written  by  interference  of  high-order  diffracted  beams  Us¬ 
ing  these  assumptions,  we  substitute  Eqs.  (22M25)  into 
Eqs.  (12M15)  and  group  terms  with  exp(2iAx)  From  this 
we  find  the  time  dependence  of  the  second-order  space- 
charge  field  to  be 
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d£,  £, 

,  =  -  -  A(K)mE]  -  B(/0£f,  (26) 

dr  To 

where 

1  eMriu  1  +  ( kTt/e2NT)i2K )J 
r,  ~  <  1  +  (nkT/yNA)(2K)2 

A(K )  and  B(K)  are  functions  of  K,  and  Et  is  the  first-order 
space-charge  field  derived  previously.11 12 

For  light-induced  erasure  of  a  previously  written  grating, 
the  general  solution  of  Eq.  (26)  is 


where  we  have  used  the  fact  that  the  first-order  field  decays 
at  a  rate  exp(—  t/r).  This  shows  that  in  general  the  erasure 
rate  of  the  second  order  is  nonexponential.  However,  de¬ 
pending  on  the  relative  sizes  of  r:,  A\,  and  B\,  the  decay  rate 
could  resemble  that  of  a  single  exponential  with  a  decay  time 
of  r2.  Also,  because  the  first-order  space-charge  field  is 
proportional  to  m,  we  see  from  Eq.  (26)  that  the  second- 
order  space-charge  field  is  proportional  to  m2. 

Similarly,  using  the  results  for  the  first-  and  second-order 
space-charge  fields  and  applying  third -order  perturbation 
theory,  we  find  that  the  expression  governing  the  time  de¬ 
pendence  of  the  third-order  space-charge  field  is 

dE  E, 

=  -■  -£(£,*.£,£,).  (29) 

at  r:! 

where  F  is  a  complex  function  of  the  first-  and  second-order 
fields.  This  shows  that  in  general  the  decay  rate  of  the  third 
order  is  also  nonexponential.  However,  when  the  driving 
term  F  is  small  the  grating  can  show  an  exponential  decay 
with  a  single  decay  time  r  (,  where 

1  +  (kl_  (3 K)- 
1  cnn„  e\V7 

7  i  '  1  +  M kJ  (3/0- 

Now,  because  the  complex  function  F  contains  terms  such  as 
Ei3  and  E\E2.  we  see  from  Eq.  (29)  that  the  third-order 
space-charge  field  is  proportional  to  m \  It  is  interesting  to 
note  that  when  the  decay  rates  of  the  three  orders  are  single 
exponentials,  the  decay  times  have  the  same  form:  the  first 
order  is  a  function  of  K,  the  second  order  is  a  function  of  2 K, 
and  the  third  order  is  a  function  of  3 K. 

BaTiO.!,  SBN,  and  BSKNN  crystals  were  cut  along  the 
(100)  axis,  polished  to  optical  quality,  and  poled  to  obtain 
approximately  single-domain  material.  In  particular,  the 
BaTiO  (  crystal  was  grown  from  a  melt  containing  high-puri¬ 
ty  BaO  and  TiO_,  with  extra  precautions  taken  to  minimize 
melt  contamination  from  the  furnace. 

High-order  diffraction  efficiencies  and  decay  rates  of  an¬ 
isotropic  diffraction  were  measured  for  the  high-puritv  sam¬ 
ple  of  BaTiO.  up  to  third  order.  The  writing  and  erase 
beams  all  originated  from  an  argon-ion  laser  (Coherent  In¬ 
nova  70)  operating  in  the  TEMnn  mode  with  X  =  488  nm. 
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The  readout  beam  (A  =  633  nm)  was  chopped  at  a  frequency 
of  1000  Hz,  and  the  diffracted  beam  was  sychronousiy  de¬ 
tected  by  using  a  photomultiplier  tube  and  a  lock-in  amplifi¬ 
er.  The  output  of  the  lock-in  amplifier  was  fed  to  an  AT&T 
6300  computer  equipped  with  a  Data  Translation  data-ac- 
quistion  board  with  the  sampling  rate  set  at  80  Hz. 

The  light-induced  decay  rates  of  the  first  three  orders  in 
BaTi(>3  were  found  to  be  well  fitted  by  a  single  exponential. 
Beam-coupling  measurements  on  this  sample  determined 
the  dominant  conduction  mechanism  to  be  holes  with  practi¬ 
cally  no  hole-electron  competition.  Nonexponential  decays 
were  observed  for  high-order  anisotropic  diffraction  in  the 
SBN  and  BSKNN  crystals;  thus  our  most  detailed  experi- 


Fig  3.  Decay  rates  of  the  first  three  diffraction  orders  as  a  function 
of  grating  vector  for  a  BaTiOi  crystal  (a)  The  first  order  is  plotted 
versus  K2,  (bl  the  second  order  is  plotted  versus  (2K):.  and  (c)  the 
third  order  is  plotted  versus  (3K>2.  The  filled  circles  are  the  data 
points,  and  the  solid  curves  are  the  theoretical  fit 


Fig.  4.  Square  root  of  the  diffraction  efficiency  plotted  versus  the 
modulation  index  m  for  the  three  orders  of  anisotropic  diffraction  in 
BaTiO  ,.  The  first,  second,  and  third  orders  are  plotted  as  filled 
circles,  open  triangles,  and  open  circles,  respectively.  The  solid 
curves  are  the  theoretical  fits  to  the  linear  regions. 

ments  were  conducted  on  the  BaTiO.i  crystal.  The  decay- 
rates  of  the  three  orders  were  measured  as  a  function  of 
grating  vector  and  were  fitted  to  Eqs.  (21),  (27),  and  (30)  by- 
using  a  nonlinear  least-squares-fitting  routine.  The  writ¬ 
ing-beam  intensities  were  /)  =  0.2  mW/cm2  and  h  =  31  mW/ 
cm2  for  the  first  order,  I\  =  7.3  mW/cm2  and  h  =  21  mW/cm2 
for  the  second  order,  and  I\  —  19  mW/cm2  and  h  -  34  mW/ 
cm2  for  the  third  order.  Because  grating  decay  rates  in 
BaTiC>3  usually  exhibit  a  nonlinear  intensity  dependence, 
the  erase- beam  intensity  was  held  constant  throughout  the 
experiments  at  11.5  mW/cm2  for  the  first  order,  4.5  mW/cm2 
for  the  second  order,  and  4  mW/cm2  for  the  third  order.  The 
readout-beam  intensity  Ir  was  kept  below  1  /iW/cm2  in  all 
experiments. 

The  results  of  the  light-induced  grating-erasure-versus- 
grating-vector  experiments  are  shown  in  Fig.  3.  A  nonlinear 
dependence  of  the  decay  rate  on  the  square  of  the  grating 
vector  was  observed  for  anisotropic  diffraction  into  the  first 
and  second  orders.  The  theoretical  fit  of  Eqs.  (21),  (27),  and 
(30)  to  the  data  is  also  plotted  as  solid  lines  in  Fig.  3.  This 
gives  effective  trap  densities,  St,  of  (1.5  ±  0.6)  X  1016/cm\ 
(0.65  ±  0.7)  X  1016/cnv  ,  and  (2.1  -k  0.4)  X  lO'Vcm3  from  the 
first-,  second-,  and  third-order  data,  respectively.  Taking 
the  uncertainties  into  account,  we  find  that  these  values  are 
in  fair  agreement. 

The  diffraction  efficiency  was  also  measured  as  a  function 
of  the  modulation  index  for  each  order.  In  this  experiment 
one  of  the  writing-beam  intensities  was  kept  constant  at  27 
mW/cm2  and  the  intensity  of  the  other  varied.  The  square 
root  of  the  diffraction  efficiency  is  plotted  versus  m  in  Fig.  4 
for  the  three  orders.  In  this  model  the  square  root  of  the 
diffraction  efficiency  is  linearly  proportional  to  the  space- 
charge  Held  in  the  absence  of  beam  coupling  and  for  low- 
diffraction  efficiencies  (Eq.  (9)].  Fitting  our  data  for  the 
three  orders  to  this  equation,  we  find  the  values  of  M  in  mM 
to  be  1,  1.9,  and  3.2  for  the  first,  second,  and  third  orders, 
respectively.  We  see  that  the  theoretical  Fit  (solid  line  in  the 
figures)  is  in  good  agreement  with  the  measured  efficiencies 
for  the  first  two  orders  with  a  reasonable  fit  for  the  third 
order  at  small  m. 

This  result  is  in  contrast  to  the  results  of  Huignard  and 
LeduM  for  Bi^SiOso-  They  found  that  the  diffraction  effi¬ 
ciencies  of  the  different  orders  followed  the  equation 
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(l-m3)1'1 

m 


(31) 


We  find  that  our  data  do  not  fit  this  equation  for  large  m. 


4.  ANISOTROPIC  SELF-DIFFRACTION 

Anisotropic  self-diffraction  occurs  when  the  phase-match¬ 
ing  conditions  are  satisfied  for  the  writing  beams.  The  con¬ 
figuration  that  we  used  to  study  anisotropic  self-diffraction 
is  shown  in  Fig.  5.  Here,  two  extraordinary  polarized  writ¬ 
ing  beams,  A(  and  A>,  of  wave  vector  ki  and  k2  intersect  in 
the  x-y  plane  of  a  4-mm  crystal  to  form  an  anisotropic 
grating  of  wave  vector  K  =  kj  -  k2.  For  plane-wave  light 
fields  of  the  form 

A(r,  t)  =  eA(r)  exp[i(k  •  r  -  wt) j  +  c.c.,  (32) 

using  the  slowly  varying  field  approximation  we  can  derive 
the  following  coupled-wave  equations  governing  self-diffrac¬ 
tion: 


Fig.  6.  Phase-matching  diagrams  showing  examples  of  anisotropic 
self-diffraction  of  extraordinary  polarized  beams  into  ordinary  po¬ 
larized  (a)  first  and  second  orders  and  (bi  third  and  fourth  orders. 


8A, 

— 2tA4  cos(0J  ^  =x1M3  +  »4A4, 

(33) 

dA, 

-2**2  C0S^,n>  6v'  =  *-3A4  + 

(34) 

6A3 

-21*3  cos(i£,n)  *  <fjA,  +  K2A,, 

dy 

(35) 

dA. 

—2ikA  cos(^in)  --  =  *3*Ao  +  x4*.44. 

d>- 

(36) 

Here  8m  and  ^,n  are  the  angles  of  the  writing  and  diffracted 
beams,  respectively,  inside  the  crystal.  The  coupling  con¬ 
stants  are  given  by 


Fig.  5.  Setup  for  observing  anisotropic  self-diffraction  in  photore- 
tractive  crystals  The  angle  4  is  the  phase  matching  angle  defined 
in  Eq  (43) 


Table  1.  Phase-Matching  Angles  of  the  Writing 
Beams  Outside  the  Crystal  for  High-Order 
Anisotropic  Self-Diffraction 


Orders 

BaTi03 

SBN 

BSKNN 

First  and  second0 

12 

8.5 

10.2 

Second  and  third 

6.9 

4.9 

6.2 

Third  and  fourth 

4.9 

3.5 

4.2 

Fourth  and  fifth 

— 

— 

3.4 

0  Measured  values 


X,  =  -<  •'IR  •  Em  •  ?,  (41) 

X2«  -i  (42) 

where  Ey  and  £.vr+i  are  the  Af-  and  Af  +  1 -order  space- 
charge  fields,  is  the  free-space  permittivity,  and  c  is  the 
free-space  speed  of  light.  We  also  see  that  although  the 
writing  beams  A,  and  A_>  do  not  experience  self-diffraction 
into  each  other,  they  are  indirectly  coupled  because  £,vr  is 
proportional  to  ( / 1 /2)w  -. 

The  anisotropic  self-diffraction  phase-matching  diagram 
is  shown  in  Fig.  6(a).  Close  examination  of  this  figure  clear¬ 
ly  shows  that  each  diffracted  beam  is  composed  of  two  or¬ 
ders.  For  instance,  beam  A  t  is  the  result  of  self-diffraction 
of  beam  A!  into  first  order  (k:)  =  k)  +  K)  and  beam  A2  into 
second  order  (k  .  =  k2  +  2K).  This  can  be  observed  experi¬ 
mentally  simply  by  blocking  beam  Ai  or  k  >  and  noting  the 
presence  of  the  diffracted  beam  A  i.  Thus  a  measure  of  the 
rise  or  decay  rate  of  these  diffracted  beams  involves  more 
than  one  time  constant.  Because  the  diffraction  intensity  of 
the  orders  goes  as  m2M,  where  m  is  the  modulation  index  and 
M  is  the  order  of  the  diffracted  beam,  the  departure  from  a 
single  time  constant  should  be  significant  only  for  large  m 
gratings. 

Figure  6(b)  shows  typical  phase-matching  conditions  re¬ 
quired  to  observe  anisotropic  self-diffraction  into  orders 
higher  than  the  second  order.  In  the  particular  case  shown 
the  diffracted  beams  are  composed  of  both  third  and  fourth 
orders  (k-t  =  k*  +  3 K  =  k2  +  4K).  In  general,  anisotropic 
63  self-diffraction  occurs  in  pairs  of  order —Af  and  Af  +  1.  The 
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phase-matching  conditions  can  be  expressed  mathematical¬ 
ly  by 


sin(ft)  = 


n  -  -  nk: 
(2M  +  llJ  -  1 


(43) 


sin(^)  =  (2A/  +  1 ) 


n  -  -  n 
<2Af  +  1): 


(44) 


where  ft  and  ^  are  the  angles  outside  the  crystal  of  the  writing 
and  diffracted  beams,  respectively.  If  the  high-order  space- 
charge  fields  are  neglected,  Eqs.  (33)-(36)  and  (43)  and  (44) 
reduce  to  those  of  Kukhtarev  et  al.- 

High-order  self-diffraction  experiments  were  performed 
on  the  crystals  using  a  Coherent  Innova  70  argon-ion  laser, 
with  A  =  488  A,  operating  in  the  TEMi»,  mode.  The  beams 
were  roughly  collimated  to  approximately  5  mm  in  diameter 
with  an  intensity  ratio  I\/l>  =  3. 

First,  the  angles  for  anisotropic  self-diffraction  into  first 
and  second  orders  (A/  =  1 ),  as  in  Fig.  6(a),  were  measured  for 
each  sample.  These  were  ft  =  12°,  8.5°,  and  10.8°  for  the 
BaTiO),  SBN\  and  BSKNN  crystals,  respectively.  From 
Eq.  (43)  these  values  of  ft  imply  that  n, 2  —  n,,-  =  0.346  for 
BaTiO),  0.1734  for  SBN,  and  0.281  for  BSKNN.  I'sing 
these  values  of  n,,2  -  n,,2  in  Eq.  (43),  we  calculated  the  angles 
of  the  writing  beams  needed  to  observe  the  higher  orders. 
These  are  shown  in  Table  1.  Note  that  up  to  the  fifth  order 
was  observed  in  BSKNN.  The  measured  angles  were  all 
within  the  experimental  error,  ±0.1°  of  the  calculated  val¬ 
ues. 


5.  SUMMARY 

In  conclusion,  we  have  observed  high-order  anisotropic  self¬ 
diffraction  in  BaTiO),  SBN,  and  BSKNN.  The  phase¬ 
matching  angles  were  measured  and  found  to  be  in  excellent 
agreement  with  the  theory.  Also,  it  was  shown  that  the 
beams  that  are  diffracted  through  anisotropic  self-diffrac¬ 
tion  are  composed  of  two  distinct  orders.  To  model  high- 
order  anisotropic  diffraction,  we  included  high-order  terms 
in  Kukhtarev's  solution  of  the  transport  equations.  The 
light-induced  decay  rates  versus  grating  vector  squared  and 
the  diffraction  efficiencies  versus  modulation  index  were 
measured  for  the  first  three  orders  in  BaTiO-)  and  were 
found  to  be  in  fair  agreement  with  the  model. 


ACKNOWLEDGMENTS 

This  research  was  supported  in  part  by  the  Defense  Ad¬ 
vanced  Research  Projects  Agency  under  U.S.  Air  Force  con¬ 
tract  F19628-85-K-0050.  The  SBN  and  BSKNN  crystals 
were  supplied  by  R.  R.  Neurgaonkar  of  Rockwell  Interna¬ 
tional.  The  high-purity  BaTiOt  crystal  was  grown  in  the 
MIT  Crystal  Physics  and  Optical  Electronics  Laboratory. 

The  authors  are  also  with  the  Departments  of  Electrical 
Engineering  and  Computer  Science,  Massachusetts  Insti¬ 
tute  of  Technology. 

REFERENCES 

1.  S.  I.  Stepanov,  M.  P.  Petrov,  and  A.  A.  Kashilm.  "Optical  dif¬ 
fraction  with  polarization-plane  rotation  in  a  volume  hologram 
in  an  electrooptic  crystal,"  Sov.  Tech.  Phys.  Lett.  3,  345  1 1977) 

2.  N.  V’.  Kukhtarev,  E.  Kratzig,  H  C.  Kulich,  and  R.  A  Rupp. 
"Anisotropic  self-diffraction  in  BaTiO,,"  Appl.  Phys.  B  35,  17 
(1984). 

3.  R.  A.  Rupp  and  F.  W.  Drees,  “Light-induced  scattering  in  pho- 
torefractive  crystals,"  Appl.  Phys.  B  39,  223  (1986). 

4.  E.  Voit,  C.  Zaido,  and  P.  Gunter.  “Optically  induced  variable 
light  deflection  by  anisotropic  diffraction  in  photorefractive 
KN’bOi,”  Opt.  Lett.  11,  309  (1986). 

5.  E.  Voit  and  P.  Gunter,  "Photorefractive  spatial  light  modula¬ 
tion  by  anisotropic  self-diffraction  in  KNbNOa  Crystals."  Opt. 
Lett.  12,  10(1987). 

6.  N.  V.  Kuktharev  and  S.  G.  Odulov,  “Wave  front  inversion  in 
anisotropic  self-diffraction  of  laser  beams,"  Sov.  Tech.  Phys. 
Lett.  6,  503  (1980). 

7  N.  V.  Kuktharev,  "Wavefront  reversal  of  optical  beams  in  an 
anisotropic  media,”  Sov.  J  Quantum  Electron.  11, 878  (1981 ). 

8.  D.  A  Temple  and  C.  Warde,  “Amplified  reflection  through 
anisotropic  four-wave  mixing  in  BaTiO), "  submitted  to  Opt. 
Lett. 

9.  D.  A.  Temple  and  C  Warde.  "Anisotropic  scattering  in  photore¬ 
fractive  crystals."  J.  Opt.  Soc.  Am.  B  3,  337  (1986) 

10.  P.  Refregier,  L.  Solymar.  H.  Rajbenbach,  and  J.  Huignard, 
“Beam  coupling  in  photorefractive  Bii.-SiO-jo  crystals  with  mov¬ 
ing  gratings.  Theory  and  experiments,"  J.  Appl  Phys.  58,  949 
(1985). 

11.  N.  V.  Kukhtarev.  V  B.  Markov.  S.  G.  Odulov.  M.  S.  Soskin.  and 
V.  L.  Vineteskii,  “Holographic  storage  in  electrooptic  crystals  I. 
Steady  state,"  Ferroelectrics  22, 949  (1979). 

12.  J.  Feinberg,  D  Heiman.  A.  R.  Tanguay,  Jr.,  and  R  W.  Hell- 
warth,  "Photorefractive  effects  and  light  induced  charge  migra¬ 
tion  in  barium  titanate."  -)  Appl  Phys.  51,  1297  (1980). 

13.  H.  Kugelnik.  "Coupled  wa\e  theory  for  thick  hologram  gra¬ 
tings."  Bell  Syst  Tech.  J.  48.  2909  ( 1969). 

14.  J  P,  Huignard  and  B  Ledu.  "Collinear  Bragg  diffraction  in 
photorefractive  Bii.SiOj,.."  Opt.  Lett.  7,  310  (1982). 


64 


8  Amplified  Reflection  Through  Anisotropic  Four- 
Wave  Mixing  in  BaTiO^ 


6  i 


Reprinted  from  Optics  Letters,  Voi.  13.  page  764,  September  1988. 

Copyright  ©  1988  by  the  Optica)  Society  of  America  and  reprinted  by  permission  of  the  copyright  owner. 


Amplified  reflection  through  anisotropic  four-wave  mixing  in 

BaTi03 


D.  A.  Temple 

Depart/m-nt  o*  Physics.  Massachusetts  institute  of  Technology.  Cambridge.  Massachusetts  i’LT.W 

].  Chang 

Department  of  Materia/ Science.  Massachusetts  /nstitute  of  Tn  hn>Jog\.  Cambridge.  Massachusetts  02 l.W 

C.  Warde 

Department  of  Electrical  Engineering  and  Eomputrr  Si  leni.e.  Massachusetts  Institute  of  Technology.  Cambridge.  Massachusetts  021  JO 


Ref rived  |,inikiry  27.  1988.  accepted  )une  20  1988 

We  have  demonstrated  cw  non  -phase-conjugate  amplified  reflect  ion  in  single-crystal  BaTiOi  using  anisotropic 
four- wave  mixing  We  show  that  the  reflected  beam  is  polarized  orthogonally  to  the  readout  beam  and  that  it  is  a 
true  reflection  of  the  image  beam,  not  a  phase  conjugate.  The  maximum  reflectivity  measured  was  approximately 

200rr. 


Amplified  reflection  has  been  observed  previously  in 
BaTi03  using  isotropic  four-wave  mixing.1  1  In  the 
normal  configuration  shown  in  Fig.  1(a).  the  object  A i 
and  pump  A 2  beams  intersect  in  the  crystal  to  form  a 
photorefractive  grating.  A  second  pump  beam  A 1 
counterpropagating  with  the  first  pump  beam  scatters 
from  this  grating  into  a  fourth  beam  A4,  which  has  the 
same  polarization  as  the  pump  beam  A;  and  is  the 
phase  conjugate  vT  the  object  beam  A,.  In  BaTiO,  the 
phase-conjugate  reflectivity  has  been  found  to  exceed 
10(Ft.1 

The  configuration  for  observing  amplified  reflection 
through  anisotropic  four-wave  mixing  in  BaTiO,  is 
shown  in  Fig.  1(b).  As  seen  from  the  phase-matching 
diagram  (Fig.  2),  it  is  a  simple  variation  of  the  configu¬ 
ration  used  for  observing  anisotropic  self-diffrac- 
tion.4  6  In  this  system,  two  coherent  ordinary-polar¬ 
ized  plane  waves  Ai  and  A->  with  wave  vectors  k|  and  kj 
intersect  in  the  x-y  plane  of  the  crystal,  producing  an 
anisotropic  grating  through  the  photorefractive  effect. 
An  ordinary-polarized  plane  wave  A  1  of  wave  vector  k  1 
is  used  to  read  out  the  grating,  as  in  Fig  2.  producing 
an  extraordinary-polarized  wave  A4  of  wave  vector  k4. 
Unlike  the  normal  phase-conjugate  reflection  shown 
in  Fig.  1(a),  the  unique  features  of  this  system  are  1 1 ) 
the  readout  beam  is  always  ordinary  polarized  and  the 
reflected  beam  is  always  extraordinary  polarized,  (2) 
ordinary-polarized  writing  beams  cannot  couple  to  the 
grating,  which  eliminates  beam-coupling  effects,*’  (3) 
the  reflected  beam  is  a  real  reflection  of  the  image 
beam  A|,  not  a  phase  conjugate,  and  (4)  anisotropic 
diffraction  in  BaTiO;  makes  use  of  the  r,_  electro¬ 
optic  coefficient  exclusively.  The  major  disadvan- 
tages  of  this  mirror  are  that  strict  phase-matching 
conditions  must  be  met  in  order  to  achieve  maximum 
efficiency  and  phase-conjugate  reflection  i'  not  possi¬ 
ble  in  this  configuration. 


To  model  this  mirror,  consider  the  three  ordinary- 
polarized  plane  waves  A,.  A-j,  and  A:(  intersecting  in 
the  x-y  plane  of  a  photorefractive  crystal  as  in  Fig. 
Kb).  Although  there  are  three  gratings  formed  by 
interference  between  beams  A,  and  Aj,  A.  and  A:i,  and 
A_>  and  A  t,  the  ordinary-polarized  writing  beams  Ai 
and  A-i  are  not  properly  phase  matched  to  couple  to 
any  of  these  gratings. *  Hence,  the  only  allowed  cou¬ 
pling  is  diffraction  of  the  readout  beam  A into  beam 
A4.  Therefore  we  only  need  to  consider  the  grating 
produced  by  interference  between  beams  Ai  and  A_>, 
where  K  =  k, .  —  k|.  Now  .  by  using  the  slowly  varying 
field  approximation,  H  we  can  derive  the  coupled- 
wave  equations  for  this  process.  For  plane  waves  of 
the  form 


E(r,  f  1  =  Airiexp[i(k  •  r  -  u.f )]  +  c.c..  ( 1 ) 


the  coupled  equations  are 

dA 

\ 

d\  "H 

(  v. . 

1A,  expuo). 

(2) 

dA, 

/A  *A , ' 

\ 

d> 

A ,  exp(io). 

(3) 

The  coupling  coefficients  are  given  by 


r'te  •*  •  e.  He  ,  ■ 
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where  A  and  are  the  wave  vector  and  frequency  of 
the  beams,  respectively.  /,,  =  I\  +  /_.  ty.  e_.  e  .  and  e, 
are  the  polariz.u ion  vectors.  X  is  the  wavelength.  «  and 
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beam?  .4  and  ,4_  write  a  photoretractive  gr  itmg.  1  hi?  grat¬ 
ing  i?  read  nut  hv  beam  A  ,  which  scatter-  i:.i.-  n-am 


*  K, 


Fig  2.  Phase-matching  diagram  for  observing  anisotropic 
tour  wave  mixing  in  a  negative  birefringenr  photoretractive 
i rystal 


■y  are  the  angles  of  the  writing  and  readmit  beams,  o  is 
the  phase  between  the  intensity  fringes  and  the  index¬ 
grating  planes,  and  n  and  'i.  are  the  ordinary  and 
extraordinary  indices  of  refraction.  I  he  effective  sus¬ 
ceptibility  is  given  by 

Xe!l  =  •'r-E.).?,  '*>' 


where  r  is  the  second-rank  optical  dielectric  constant. 
r  is  the  third-rank  electro-optic  coefficient,  and  E.  is 
the  space-charge  field.  From  Eqs.  t2>  and  id>  we  see 
that  the  reflected  beam  is  a  standard  reflection  of  the 
image  beam  .4,  ^  .4;,  propagating  back  along  the  di¬ 
rection  of  the  image  beam.  kt  =  ~k;.  Thus,  the  crys¬ 
tal  acts  as  if  it  were  a  standard  mirror  with  its  surface 
normal  parallel  to  k(.  This  is  similar  to  the  phase- 
conjugation  case  except  that  the  reflection  is  a  stan¬ 
dard  reflection.  When  the  beam  ratios  are  such  that 
/.,  / ,  »  /.,  the  reflectivity  R  can  be  approximated  by 


-  V 

sin  t  v  * c < i*‘ * 2" ) 

=  n  rt, 

1  V  A  / 

cost  H  l(  ,  >?(  y  1 
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L  J 
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i  ,  i 


where  /  i?  the  interaction  length  o*  the  readout  beam 


/,.  Since  the  angles  are  fixed  owing  to  the  phase¬ 
matching  conditions,  this  relation  shows  that  the  re¬ 
flectivity  can  only  be  maximized  by  optimizing  the 
beam  overlap,  a  severe  limitation  for  this  system.  For 
small  modulation  (m  «  1)  the  space-charge  field  £„c 
can  be  approximated  by-  “ 

{kT/e)K  <8) 
1  +  (ikT/e~.\r )K~ 

where  e  is  the  charge  on  an  electron,  e  is  the  dielectric 
constant,  k  is  the  Boltzmann  constant.  T  is  the  tem¬ 
perature.  K  =  Uir/Xlsin  (*  is  the  magnitude  of  the 
grating  vector,  and  X r  is  the  trap  density.  I'sing  a 
trap  density  of  1.5  X  10ih  cm~ ',  which  was  determined 
previously,  values  for  the  refractive  index  and  dielec  ¬ 
tric  constant  determined  by  Wemple  et  a/..1"  r4J  =  820 
m  V  i  Ref.  11),  and  an  interaction  length  of  approxi¬ 
mately  4  mm  in  Eq.  (7).  we  estimate  the  maximum 
reflectivity  to  be  approximately  4 GOT.  This  is  in  rea¬ 
sonable  agreement  with  the  maximum  reflectivity  that 
we  were  able  to  attain,  which  was  2009.  The  differ¬ 
ence  is  probably  due  to  our  estimate  of  the  interaction 
length  l  or  to  the  fact  that  Eq.  (6)  neglects  reflections 
at  the  crystal  surface. 

Extraordinary-polarized  writing  beams  can  also  be 
used  to  write  the  gratings.  However,  in  this  case  the 
process  is  greatly  complicated  by  self-diffraction  of 
each  writing  beam  into  both  a  first  and  second  order. s 
Qualitatively  we  expect  this  to  cause  the  grating 
planes  to  vary  in  the  crystal  owing  to  the  large  varia¬ 
tion  of  the  writing-beam  phases  and  intensities,  as 
each  is  diffracted  into  the  two  orders.  This  should 
have  the  effect  of  reducing  the  grating  strength  and 
therefore  the  reflectivity.  For  this  case,  we  measured 
a  reflectivity  of  approximately  50rc  for  the  same  writ¬ 
ing-  and  readout-beam  intensities. 


Fig.  3  Optical  -etups  used  to  studs  anisotropic  tour  wave 
mixing  In  setup  on.  ail  the  waves  are  collimated  and  the 
reflet livitv  vu-  measured  a?  a  function  of  the  pump  beam 
ratio.  Setup  'to  was  ust-d  to  demonstrate  mirror  action 
using  a  resolution  c  hart  as  the  transparence  T  The  angle? 
were  «  =  «  =  12.0°  and  v  =  4'te.  Lenses  shown  have 

different  lot  al  a  ncth? 
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Fig.  4.  Graph  of  the  mirror  reflectivity  as  a  function  of  the 
pump-beam  ratio  I  JI  i.  The  intensities  of  the  writing  beams 
7i  and  /?  were  kept  constant,  with  the  filled  circles  represent¬ 
ing  1\  =  850  *iW  and  7  =  .3.4  m\V  and  the  open  circles 
representing  7i  =  100  fiW  and  I<  =  20  m\V. 


Fig.  5.  Photograph  of  the  object  after  anisotropic  reflec¬ 
tion. 


The  systems  used  for  studying  amplified  reflection 
through  anisotropic  four-wave  mixing  are  shown  in 
Fig.  3.  All  the  beams  originated  from  an  argon-ion 
laser  operating  at  488  nra  in  the  TEM,*,  mode.  For  the 
configuration  in  Fig.  3ta),  the  reflectivity  IJI\  was 
measured  as  a  function  of  the  readout -pump-beam 
intensity  h.  The  writing-beam  intensities  were  held 
constant,  and  all  the  beams  were  collimated  in  an 
effort  to  approximate  the  theory.1  The  results  of 
these  measurements  are  plotted  in  Fig.  4,  with  the 
filled  circles  representing  7t  =  850  p W  and  7j  =  3.4  m\V 
and  the  open  circles  representing  7i  =  100  pVV  and  I>  = 
20  mW.  The  lines  are  the  theoretical  fits  to  Eq.  (7). 
As  shown,  the  reflectivity  was  found  to  be  a  fairly  good 
fit  to  the  model  in  the  region  of  undepleted  pumps. 

The  reflecting  properties  of  this  mirror  were  investi¬ 
gated  by  using  the  system  shown  in  Fig.  3ibi.  Here, 


the  object  beam  I\  was  expanded  and  recollimated, 
and  a  resolution  chart  placed  in  the  path  was  imaged 
with  a  lens  L  (focal  length  38  cm)  through  a  beam 
splitter  BS  onto  the  crystal.  The  crystal  was  slightly 
rotated  to  avoid  overlap  of  the  surface  reflection  of  A  ■> 
with  the  anisotropic  reflected  image  A4.  A  photo¬ 
graph  of  an  amplified  reflected  image  is  shown  in  Fig. 
5.  The  lens  in  front  of  the  film  was  used  for  reimaging 
and  magnification;  its  focal  length  was  much  different 
from  that  in  front  of  transparency  T.  Three  separate 
observations  confirmed  the  reflecting  properties  of 
this  mirror:  (1)  when  an  aberrator  was  placed  in  the 
path  of  beam  1  between  the  beam  splitter  and  the 
crystal  the  reflected  image  at  the  screen  was  greatly- 
distorted;  (2)  when  transparency  T  was  removed  and 
the  lens  of  focal  length  /  =  38  cm  placed  22  cm  in  front 
of  the  crystal,  the  reflected  beam  was  found  to  focus  at 
a  distance  16  cm  in  front  of  the  crystal  and  the  trans¬ 
mitted  beam  was  found  to  focus  at  a  distance  16  cm 
behind  the  crystal;  and  (3)  when  a  conventional  mirror 
was  used  in  place  of  the  crystal  the  position  of  the  focal 
point  was  unchanged. 

In  summary,  we  have  demonstrated  non-phase-con¬ 
jugate  amplified  reflection  using  anisotropic  four- 
wave  mixing  in  photorefractive  BaTiO;i.  The  reflect¬ 
ed  beam  was  shown  to  be  a  true  reflection  of  the  image 
beam,  with  a  measured  amplification  of  approximate¬ 
ly  200%.  This  has  potential  for  improving  the  signal- 
to-noise  ratio  of  image  amplification  without  produc¬ 
ing  a  phase-conjugate  image12  or  contrast  reversal.12 

This  research  was  supported  in  part  by  the  Defense 
Advanced  Research  Projects  Agency  under  U.S.  Air 
Force  contract  F19628-85-K-0050. 
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Abstract— Unusual  beam  fanning  has  been  observed  in  heavil)  doped 
BaTiOjiCr  along  with  an  increased  photovoltaic  effect.  The  nature  and 
symmetry  of  this  beam  fanning  are  described,  and  the  results  of  the 
photovoltaic  measurements  are  presented. 

Introduction 

BARIUM  T1TANATE  (BaTiCM  is  a  photorefractive 
material  with  large  electrooptic  coefficients.  Beam 
fanning  has  been  observed  in  nominally  "undoped"  sam¬ 
ples  of  the  material,  and  this  phenomena  is  believed  to  be 
important  in  the  process  that  is  responsible  for  the  excel¬ 
lent  self-pumped  phase  conjugating  property  of  the  cry  stal 

m-[3). 

Although  models  of  the  photorefractive  effect  have  ex¬ 
isted  since  the  early  work  of  Chen  (4]-(6|.  the  details  of 
the  effect  in  BaTiO,  (i.e..  the  excitation  and  trapping  cen¬ 
ters),  and  the  exact  nature  of  beam  fanning  have  not  yet 
been  firmly  established.  Recent  papers  have  pointed  to 
iron  as  the  cause  of  the  photorefractive  effect  in  BaTiO, 
[7J,  but  as  we  shall  see  herein  it  is  not  the  only  dopant 
that  can  alter  the  photorefractive  properties  of  the  mate¬ 
rial  In  this  paper,  we  discuss  some  of  the  properties  of 
BaTiO,  doped  with  chromium 

Sample  Preparation 

A  BaTiO,  boule  was  grown  by  the  top  seeded  solution 
method  (81  with  20  ppm  chromium  in  the  melt.  A  sample 
was  cut  from  the  boule,  polished,  and  electrically  poled 
to  yield  a  crystal  with  final  dimensions  of  6.7  x  5.6  x 
3  6  mm.  An  undoped  sample  grown  from  a  nominally  pure 
melt  with  dimensions  4.1  x  4.9  x  4  1  mm  was  used  for 
comparison 

Diffraction  Efficiency 

Attempts  were  made  to  measure  the  rise  time  and  sat¬ 
uration  level  of  the  diffraction  efficiency  of  the  Ba¬ 
TiO,  :  Cr  sample  by  writing  a  grating  with  an  Ar-ion  laser 
(wavelength  488  nm)  and  reading  the  grating  with  a  HeNe 
laser  (wavelength  632.8  nm).  These  attempts  were  frus¬ 
trated  by  the  observation  that  the  diffracted  intensity  did 
not  nse  to  a  steady-state  value  as  was  the  case  for  the 
nominally  pure  sample.  Instead,  the  diffracted  beam  in¬ 
tensity  initially  rose,  but  then  decayed  toward  zero,  as 
shown  in  Fig  1 . 
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Fig  1  Diffraction  efficiency  as  a  function  of  time  for  a  BaTiO,  Cr  cry  s¬ 
tal  The  writing  beams  were  on  continuously  from  i  =  0 

Observation  of  the  crystal  during  the  writing  process 
revealed  that  light  from  the  writing  beams  was  being  scat¬ 
tered  out  of  the  crystal  along  well-defined  paths  (beam 
fanning).  This  scattered  light  became  noticeable  at  ap¬ 
proximately  the  same  time  that  the  intensity  of  the  dif¬ 
fracted  HeNe  laser  beam  began  decaying,  and  the  rise  in 
intensity  of  the  scattered  light  continued  to  correlate  with 
the  funher  decay  in  the  intensity  of  the  diffracted  beam. 
Clearly,  while  some  beam  fanning  is  essential  for  self- 
pumped  phase  conjugation,  strong  beam  fanning  can  be 
detrimental  since  it  may  affect  or  limit  the  diffraction  ef¬ 
ficiency  as  illustrated  in  Fig.  1. 

Beam  Fanning 

Beam  fanning  in  the  Cr  doped  BaTi03  crystal  was  stud¬ 
ied  using  a  collimated  (2  mm  diameter)  Ar-ion  laser 
beam,  with  a  wavelength  of  488  nm.  Three  cases  of  beam 
fanning  for  a  single  pump  beam  initially  incident  along 
one  of  the  high  symmetry  axes  of  the  crystal  are  discussed 
below:  1)  initial  propagation  antiparallel  to  the  c  axis,  2) 
initial  propagation  parallel  to  the  c  axis,  and  3)  initial 
propagation  parallel  to  an  a  axis.  The  c  axis  was  found 
optically  using  the  crossed  polarizer  technique,  and  the 
+  c  axis  was  determined  by  observing  the  voltage  devel¬ 
oped  across  the  c  axis  when  the  crystal  was  heated.  The 
face  which  becomes  negative  upon  heating  the  crystal  in 
the  dark  is  the  tc  face. 

The  most  sinking  beam  fanning  was  observed  for  the 
case  of  the  incident  beam  propagating  antiparallel  to  the 
c  axis  of  the  cry  stal  Initially  the  beam  propagated  through 
the  crystal,  bu.  after  a  short  time  two  distinct  beam  paths 
were  observed  within  the  crystal  diverging  from  the  in¬ 
cident  beam.  Fig  2(a)  is  a  photograph  of  the  crystal  il¬ 
lustrating  this  phenomenon.  The  fanning  beams  are  ob¬ 
servable  through  light  scattered  from  crystal 
imperfections.  At  steady  state,  most  of  the  incident  power 
was  diverted  into  these  two  paths.  The  two  paths  were 
always  observed  to  lie  in  the  plane  of  polarization  of  the 
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Fig  2  B«am  fanning  in  a  BaTiO,  Cr  crystal  (a)  Beam  propagating  an- 
tiparallel  to  the  c  axis  (bl  Beam  propagating  parallel  to  the  <  axis  (c) 
Beam  propagating  parallel  to  the  a  axis 

incident  light,  and  diverged  from  the  path  of  the  incident 
beam  close  to  the  entrant  surface.  The  two  beams  did  not 
remain  collimated  within  the  crystal.  The  fanning  beams 
diverged  within  a  range  of  angles  between  approximately 
55  and  80°  from  the  propagation  direction  of  the  incident 
beam. 

For  the  case  of  the  incident  pump  beam  propagating  in 
the  direction  of  the  +  c  axis,  no  fanning  beams  were  ob¬ 
served  inside  the  crystal  The  incident  beam,  however, 
was  strongly  attenuated  as  it  propagated  into  the  crystal 


2123 

(see  Fig.  2(b)]  and  two  strong  backscattered  beams  were 
observed  outside  the  crystal.  These  backscattered  beams 
had  properties  similar  to  the  fanning  beams  discussed 
above:  1)  the  build-up  time  is  approximately  the  same,  2) 
the  two  backscattered  beams  lie  in  the  plane  of  polariza¬ 
tion  of  the  incident  beam,  and  3)  the  backscattered  beams 
do  not  remain  collimated.  Fig.  2(b)  shows  the  beam  path 
within  the  crystal  for  this  configuration. 

The  final  case  is  that  of  propagation  along  an  a  axis. 
Since  the  a  axis  is  not  polar,  there  is  no  distinction  be¬ 
tween  antiparallel  and  parallel.  When  the  polarization  of 
the  incident  beam  is  along  the  c  axis,  as  before,  the  inci¬ 
dent  beam  initially  propagates  through  the  crystal.  Then, 
as  time  progresses,  a  single  fanning  beam,  which  is  di¬ 
rected  toward  the  +c  axis  and  lies  in  the  plane  formed  by 
the  c  axis  and  the  propagation  direction,  is  observed.  In 
steady  state,  most  of  the  power  is  transferred  to  the  fan¬ 
ning  beam.  Fig.  2(c)  shows  the  trajectory  of  this  single 
fanning  beam  inside  the  crystal.  This  is  the  same  scatter¬ 
ing  geometry  that  gives  rise  to  a  pair  of  anisotropically 
scattered  cones  of  transmitted  light  with  nulls  along  the  c 
axis  in  earlier  work  on  undoped  BaTiO,  and  BaTiO, :  Fe. 
The  apex  half  angles  of  the  cones  were  25.1°  and  37.5° 
[9],  [10],  When  the  polarization  of  the  input  beam  is  par¬ 
allel  to  the  other  a  axis  the  light  is  scattered  at  large  angles 
(approaching  90°)  to  the  incident  beam  direction.  The 
scattering  is  out  of  the  plane  formed  by  the  polarization 
and  propagation  vectors  of  the  incident  beam,  and  is  to¬ 
ward  the  c  axis.  This  case  is  more  closely  related  to  the 
configuration  used  in  self-pumped  phase-conjugate  mir¬ 
rors,  and  yields  results  similar  to  those  previously  ob¬ 
served  by  others  in  “undoped”  crystals  [  1  ]— ]3] . 

The  general  properties  of  the  fanning  beams  did  not  vary 
with  wavelength  for  the  visible  lines  (514.5,  496.5, 
488.0,  and  476.5  nm)  of  the  Ar-ion  laser.  Also  the  beam 
fanning  pattern  was  found  to  be  independent  of  the  inten¬ 
sity  of  the  incident  beam  except  that  at  lower  intensities, 
longer  times  were  required  in  order  to  observe  the  fanning 
beams. 

In  order  to  determine  whether  the  Gaussian  intensity 
distribution  of  the  laser  beam  influenced  the  beam  fan¬ 
ning,  several  beamstop  experiments  were  performed.  Ap- 
erturing  the  beam  to  a  smaller  diameter  decreased  the 
width  of  the  fanning  beam,  as  expected.  Blocking  one- 
half  of  the  beam  in  order  to  obtain  a  sharp  intensity  gra¬ 
dient  caused  a  decrease  in  the  intensity  of  the  fanning 
beam  on  the  blocked  side,  but  the  general  shape  and  di¬ 
rection  of  the  fanning  beam  remained  the  same  as  in  the 
unblocked  case.  Finally,  a  beam  stop  was  placed  in  the 
center  of  the  Gaussian  beam.  Except  for  the  reduced  in¬ 
tensity,  no  change  in  the  beam  fanning  was  observed 

No  fanning  was  observed  for  the  red  line  of  a  HeNe 
laser  It  is  not  clear  whether  this  was  due  to  the  low  power 
of  the  HeNe  laser,  the  relatively  low  absorption  of  the 
crystal  in  the  red,  or  a  combination  of  both  The  intensity 
of  the  HeNe  laser  beam  was  approximately  15  mW /cm\ 
which  is  small  compared  with  the  intensity  of  the  Ar-ion 
laser  beam  of  approximately  500  mW/cmv 
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A  photorefractive  beam-coupling  model  in  which  the 
directional  properties  of  the  fanning  beams  are  deter¬ 
mined  by  the  symmetry  of  the  crystal  ( 4  mm  ).  the  polar¬ 
ization  of  the  incoming  laser  beam  and  the  symmetry  of 
the  electrooptic  tensor  can  be  used  to  explain  the  case  for 
input  light  propagation  along  the  c  axis  We  assume  that 
the  beam  fanning  along  preferred  directions  initiates  from 
scattering  off  gratings  established  by  interference  between 
the  input  beam  and  randomly-scattered  beams.  The  initial 
fanning  beam  itself  then  also  suffers  beam  fanning,  lead¬ 
ing  to  beam  trajectories  inside  the  crystal  which  are  de¬ 
termined  by  the  polarization  of  the  incident  light,  the  sign 
of  the  dominant  charge  carriers,  the  symmetry  and  mag¬ 
nitude  of  the  electrooptic  tensor  elements,  and  the  refrac¬ 
tive  indexes  of  the  crystal  through  a  set  of  phase  matching 
and  beam  coupling  conditions 

In  cases  1)  and  2)  (propagation  along  the  <  axis),  the 
optical  properties  of  the  cry  stal  are  isotropic  about  the  in¬ 
cident  beam.  Thus  beam  fanning  in  these  two  cases  must 
be  invariant  under  a  1 80°  rotation  ot  the  crystal  about  the 
incident  beam.  Since  the  fanning  beams  are  started  by 
random  scattering  and  reaches  a  steady  state,  some  sort 
of  feedback  mechanism  is  active  It  beam  coupling  is  the 
feedback  mechanism,  then  the  scattered  beams  of  interest 
must  have  a  polarization  component  parallel  to  that  of  the 
incident  beam  (this  is  required  in  order  for  the  two  beams 
to  interfere)  Assuming  diffusion  is  the  transport  mecha¬ 
nism  (no  external  fields  were  applied),  the  photorefractive 
grating  vector  will  lie  in  the  plane  defined  by  the  propa¬ 
gation  vectors  of  the  incident  and  scattered  beams. 

In  order  to  show  that  the  strongest  fanning  beams  in  the 
crystal  must  lie  in  the  plane  of  polarization  of  the  incident 
beam,  we  must  examine  the  coupling  due  to  the  electroop¬ 
tic  effect.  Let  us  define  a  coordinate  system  with  the  +z 
axis  in  the  direction  of  the  c  axis  and  the  x  and  y  axes 
along  the  equivalent  a  and  b  axes  of  the  crystal  We  will 
assume  that  the  incident  pump  beam  is  propagating  along 
the  ;  axis  and.  without  loss  of  generality  .  polarized  along 
the  x  axis  The  interference  between  a  typical  scattered 
beam  and  the  incident  beam  will  create  a  sinusoidal  elec- 
tnc  field  with  a  grating  vector  which  is  the  difference  be¬ 
tween  the  two  propagation  vectors  Through  the  elec¬ 
trooptic  effect  this  creates  a  grating  in  the  impemuttivity 
However,  only  those  gratings  which  can  couple  energy 
from  the  pump  to  the  scattered  signal  beam  w  ill  be  rein¬ 
forced  and  build  up  to  a  steady-state  strength. 

It  can  be  shown  that  the  impermittivuy  tensor  for  this 
geometry  has  the  form 


0 

0 

AA  H 

0 

A  A  i  i 

0 

*0 

0 

+ 

0 

A  A;: 

A  Am 

0 

0 

K, 

■4  A 1 1 

A  Am 

AA,,_ 

(  1 ) 

From  this  form,  it  is  seen  that  the  r-polarized  pump  can 
transfer  energy  by  beam  coupling  only  to  beams  polarized 
in  the  re  plane  This  means  that  only  randomly -scattered 
beams  polarized  in  the  re  plane  (extraordinary  beams)  will 


be  reinforced  by  energy  diffracted  out  of  the  pump  beam 
By  considering  the  propagation  of  light  in  anisotropic  me¬ 
dia.  we  see  that  such  beams  must  either  have  their  prop¬ 
agation  vector  in  the  .re  plane  or  exactly  along  the  y  axis. 
Thus  a  pump  beam  incident  along  the  c  axis  can  generate 
a  fanning  beam  which  lies  in  the  plane  formed  by  this  axis 
and  the  polarization  vector  of  the  pump  beam,  and  the 
fanning  beam  will  be  polarized  in  this  plane.  This  is  pre¬ 
cisely  what  was  observed  experimentally  in  Fig.  2(a)  A 
second  beam  polarized  parallel  the  polarization  of  the 
pump  and  propagating  perpendi<  ular  to  the  plane  defined 
by  the  polarization  and  propagation  vectors  of  the  pump 
is  also  permitted.  This  may  partly  explain  the  intense  scat¬ 
tering  of  the  pump  beam  observed  just  inside  the  crystal 
in  Fig.  2(a)  and  Co). 

Absorption 

Fig  3  shows  the  absorption  spectra  of  the  Cr  doped  and 
undoped  crystals.  As  can  be  seen,  the  addition  of  Cr  has 
two  major  effects  on  the  absorption  1)  the  band  edge  is 
effectively  moved  toward  lower  energy,  which  gives  the 
crystal  its  reddish  color,  and  2)  the  absorption  becomes 
more  anisotropic.  There  is  little  structure  in  the  absorp¬ 
tion  curves  except  for  a  broad  peak  at  about  600  nm  in 
the  c  axis  absorption.  Without  further  data,  it  is  not  pos¬ 
sible  to  assign  the  peak  to  Cr  transitions. 

Photovoltaic  Effect 

Another  significant  characteristic  of  this  heavily  doped 
BaTiO, :  Cr  cry  stal  is  the  magnitude  of  its  bulk  photovol¬ 
taic  effect  Bulk  photovoltaic  currents  in  noncentrosym- 
metnc  crystals,  including  BaTiOj,  have  been  observed  by 
many  researchers  [1 1  ] -t  13].  Phenomenological  models  of 
the  photovoltaic  effect  have  been  presented  by  Glass  [13] 
and  Belinicher  and  Sturman  |14],  Following  the  notation 
of  Belinicher  and  Sturman.  the  current  density  j\  '  is  given 
by 

jT  =  0ljkE';'E\--"  (2) 

where  0,)k  is  the  bulk  photovoltaic  tensor,  and  E'~  and 
are  components  of  the  optical  field  of  frequency  m 
In  BaTiOi.  the  nonzero  photovoltaic  tensor  elements  are 

0in  =  0ni  —  0221  ~  02i2  =  011’  0iii  =  0ii ■  and  0 in 
=  0 122  ~  0H 

Fig.  4(a)  schematically  shows  the  arrangement  used  for 
measuring  the  bulk  photovoltaic  currents  as  a  function  of 
the  angle  0  between  the  optical  polarization  and  the  c  axis. 
From  (2),  the  current  density  along  the  c  axis  is  given  by 

jT  =  |  E  |‘  [  0)n  sin"  8  +  0 m  cos*  8]  (3  ) 

Although  allowed  by  (2),  no  measurable  current  was 
observed  along  the  a  axis  for  any  value  of  8.  The  mea¬ 
sured  c  axis  currents  normalized  by  the  light  intensity  are 
shown  in  Fig  4(b).  The  solid  line  is  a  fit  to  the  above 
functional  torm.  and  the  inset  shows  that  the  current  var¬ 
ies  linearly  with  intensity. 

Fig  3  shows  the  absorption  curves  for  both  the  Cr 
doped  and  undoped  crystals.  Because  the  absorption  of 
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Fig  3  Polarization  dependence  of  Ihe  absorption  coefficient  of  nominally 
“undoped"  BaTiO,  and  BaTiO,  Cr  cryslals 
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Fig  4  (a)  Schematic  diagram  of  the  apparatus  used  to  measure  photovol¬ 
taic  current  (b)  The  c  axis  photovoltaic  current  normalized  to  intensity 
as  a  function  of  the  angle  between  the  polarization  vector  of  the  light 
and  the  c  axis  for  a  BaTiO,  Cr  crystal  The  solid  line  is  a  numerical  fit. 
and  the  inset  shows  the  mammum  current  versus  intensity 

the  doped  crystal  is  not  negligible,  it  must  be  included  in 
the  photovoltaic  model.  Since  the  current  is  linear  with 
intensity,  the  equation  for  the  current  becomes 

i]  =  |  £  |*  W  j  dl  {Sjn  sin"  6  exp  (  -a,  / ) 

+  B333  cos2  6  exp  (  -ary/)  }  (4) 

where  W  and  L  are  the  dimensions  of  the  electrodes  per¬ 
pendicular  to,  and  along  the  direction  of  propagation  of 
the  light  beam,  respectively,  and  a!  and  a,  are  the  a  or  c 
axis  absorption  coefficients,  respectively. 

Using  the  parameters  obtained  by  curve  fitting  the  data, 
and  the  absorption  of  the  crystal,  the  photovoltaic  tensor 
coefficients  at  the  488  nm  wavelength  were  found  to  be 
/Sjjj  =  1.1  x  10'7  and  03lt  =  6.7  x  10“8  A/W.  These 
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values  are  an  order  of  magnitude  larger  than  those  mea¬ 
sured  on  the  undoped  crystal.  Dividing  03,3  by  the  ab¬ 
sorption  at  488  nm  gives  a  Glass  constant  i,  of  1.2  x 
10~*  A-cm/W  which  is  an  order  of  magnitude  larger  than 
that  reported  in  LiNbO,  [13]. 

Summary 

The  increased  absorption  of  the  Cr  doped  BaTi03  crys¬ 
tal  (whether  it  is  due  to  Cr  or  some  compensating  defect) 
is  probably  the  source  of  the  increased  photovoltaic  activ¬ 
ity.  Unlike  the  case  of  LiNbO,  [13],  dividing  the  photo¬ 
voltaic  tensor  elements  for  the  Cr  and  undoped  samples 
by  the  respective  absorption  coefficient  of  each  sample  did 
not  yield  the  same  Glass  constant;  however  further  re¬ 
search  needs  to  be  done  on  samples  doped  with  varying 
amounts  of  Cr.  Since  the  behavior  of  the  beam  fanning 
can  be  explained  by  a  photorefractive  beam  coupling 
model,  it  seems  likely  that  the  increase  in  beam  fanning 
is  simply  due  to  an  increase  in  the  two  beam  coupling 
coefficient  which,  in  turn,  is  probably  due  to  the  high  ab¬ 
sorption  of  the  Cr  doped  crystal.  Unfortunately,  standard 
beam  coupling  gain  measurements  could  not  be  per¬ 
formed  reliably  due  to  strong  beam  fanning. 
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Czochralski  grown  SBN:Fe  and  SBN:Ce  crystals  have  been  studied  by 
photoconductivity  and  temperature  dependent  conductivity  in  order  to 
understand  the  nature  of  trapping  centers  in  these  materials. 


INTRODUCTION 

The  tungsten  bronze  crystal  SrQ  6  Ba0  4  Nb,  Oe  (SBN-60)  is  a  photorefractive 
material  that  has  recently  been  grown  in  large  single  crystals.  SBN-60  has 
only  one  phase  transition  at  75  0  C,  the  paraelectric-ferroelectric 
transition,  and  therefore  is  not  damaged  by  cooling.1  Also,  since  its 
tetragonal  distortion  is  relatively  large,  compared  with  BaTiOs ,  90°  domains 
are  uncommon.  However,  like  BaTiOj,  the  origin  of  the  photorefractive  effect 
and  the  role  of  defects  and  impurities  in  the  effect  arc  unknown. 

As  part  of  a  program  to  study  the  nature  of  dopants  in  SBN-60,  Neurgaonkar 
et  al.J  have  grown  Ce  and  Fe  doped  SBN  crystals.  The  sensitivity  of  Ce  doped 
SBN  is  two  orders  of  magnitude  larger  than  the  undoped  material  while  the 
response  time  is  reduced  approximately  one  order  of  magnitude.  In  order  to 
better  understand  the  nature  of  the  charge  transport  and  trapping,  these 
materials  were  investigated  using  photoconductivity  and  temperature- 
dependent  conductivity  techniques. 


EXPERIMENTAL  RESULTS 

The  experimental  apparatus  consists  of  a  sample  chamber,  a  Kcithly  617 

electrometer,  an  Ar-ion  laser,  and  an  AT&T  personal  computer.  The  sample 

chamber  is  a  Faraday  cage  with  an  optical  window,  and  it  houses  a  teflon 

crystal  holder.  The  Faraday  cage  is  nessecary  in  order  to  reduce  electrical 

noise,  and  teflon  is  used  for  the  crystal  holder  as  its  resistivity  is  much 

greater  than  that  of  the  crystal.  The  computer  controls  a  shutter  which 
exposes  the  crystal  to  the  laser  beam,  and  it  records  the  current  measured 
by  the  electrometer.  The  dark  conductivity  was  measured  as  a  function  of 

temperature  by  heating  the  Faraday  cage  and  monitoring  the  crystal 

temperature  with  a  thermistor  probe  near  the  crystal. 

Photoconductivity  rise  and  decay  curves  were  recorded  over  a  range  of 
intensities  and  for  several  lines  of  the  Ar-ion  laser.  The  rise  and  decay 
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times  of  the  SBN:Ce  were  too  short  to  be  resolved  by  the  electrometer  so  a 
storage  oscilloscope  was  used  to  record  the  response  of  this  crystal.  In 
order  to  achieve  sufficient  signal,  the  intensity  was  set  to  the  maximum 
availible  from  the  laser.  The  photocurent  response  time  (1/e  point)  was 
found  to  have  a  decay  time  of  approximately  22  ms  and  a  rise  time  of  18  ms. 
These  response  times  are  of  the  same  order  of  magnitude  as  those  observed  in 
photorefractive  studies  of  this  material.1 

The  response  of  the  SBN.Fe  crystal  was  much  slower.  Figure  1.  shows 
typical  rise  and  decay  curves  for  the  SBN:Fe  crystal  for  two  different 
intensities  at  a  wavelength  of  488nm.  These  curves  were  fitted  to  a  sum  of 
exponentials.3  The  decay  curves  are  best  fit  to  a  sum  of  three  exponentials 
plus  a  constant  dark  current,  while  the  rise  curves  are  composed  of  the 
steady  state  current  minus  the  sum  of  two  exponentials.  The  steady  state 
photocurrent  as  a  function  of  intensity  is  shown  in  Fig  2.  From  this  plot  we 
find  that  the  photocurrent  is  given  by: 


where  ipj,  is  in  nA  and  I  is  in  mW/cm* .  The  three  time  constants  for  the 
decay  curves  are: 


Tj  -  0.5  10.2  sec 
Tj  -  2.6  10.4  sec 
t3  -  18  19  sec 

The  largest  time  constant  may  be  due  to  thermal  decay.  Since  the  decay  time 
constants  are  expected  to  change  with  temperature  and  no  attempt  was  made  to 
control  the  temperature  in  the  photoconductivity  experiments,  the 
uncertainty  in  the  above  time  constant  values  was  expected. 

The  time  constants  for  the  rise  curves  vary  with  intensity.  Further 
experiments  are  being  conducted  in  order  to  measure  the  time  constants  as  a 
function  of  temperature,  and  such  measurements  should  yield  the  depths  of 
the  traps.  An  Arrhenius  plot  [ln(i)  vs  I/T]  of  the  dark  current  for  SBN:Fe 
is  shown  in  Fig  3.  This  shows  a  trap  depth  of  approximately  0.8  eV. 

Photoconductivity  studies  indicate  that  SBN:Ce  is  significantly  faster  than 
SBN.Fe.  Based  on  the  observed  multiple  time  constant  decay  data,  the 
response  time  of  SBN:Fe  appears  to  be  associated  with  several  trapping 
levels  or  distributions  of  trapping  levels. 
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Fig.  1.  Photocurrent  rise  and  decay  in  SBN:Fe  for  laser 
intensities:  (a)  35.8  mW/cm  ,  (b)  445  mW/cmJ 
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Fig.  3.  Arrhenius  plot  of  the  dark  current  yielding  an 
activation  energy  of  0.8  eV 
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11  Collaboration  with  Related  Research  Group 

Our  collaboration  with  related  research  groups  included  : 

1.  An  exchange  of  crystals  (BaTiOa  for  SBN)  with  Dr.  R.  R.  Neurgoankar 
of  Rockwell.  We  jointly  participated  in  the  measurement  of  the  photo¬ 
conductivity  in  the  SBN  sample  (see  Section  10) 

2.  An  exchange  of  crystals  (BaTiOa  for  BSO)  with  Professor  Armand  Tan- 
guay  of  the  University  of  Southern  California.  We  measured  photolumi¬ 
nescence  in  the  BSO  samples. 

3.  Providing  a  BaTiOa  crystal  called  ”  Doyle”  to  Professor  Jack  Feinberg  of 
the  University  of  Southern  California. 

4.  Collaborating  with  and  providing  a  BaTiOs  crystal  to  and  Dr.  Char- 
maine  Gilbreath  of  the  Naval  Research  Laboratoiy  for  use  in  wavefront 
shaping  of  a  diode  laser  for  space  communications. 

5.  Technical  exchanges  with  personnel  from  Sanders  Associates  on  the  crys¬ 
tal  growth  of  barium  titanate. 
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